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The  U.S.  Navy  mobile  blood  bank  is  responsible  for  providing  blood  to  the  National  Naval 
Medical  Center  located  in  Bethesda,  MD  as  well  as  to  other  military  and  civilian  hospitals  in  the 
surrounding  area  (e.g.  VA  hospital).  The  mobile  blood  bank  travels  to  various  locations  throughout  the 
year  with  all  equipment  and  personnel  support  required  to  draw  blood  from  donors.  The  Navy  wanted  our 
team  to  study  their  system  to  determine  if  it  could  be  improved.  We  then  determined  that  our  project  for 
the  Navy  mobile  blood  bank  was  to  provide  three  products:  Recommendations  to  improve  the  actual  blood 
donating  process  as  determined  through  computer  simulation  modeling  of  their  system  (e.g.  reduce  time  in 
system):  Information  through  forecasting  where  to  go  to  get  blood  products  and  forecasting  likely 
quantities  at  each  location;  lastly,  from  the  computer  simulation,  a  chart  to  determine  more  accurately  the 
number  of  servers  required  for  a  mission  based  on  anticipated  donors. 

The  Navy  mobile  blood  bank  visits,  annually  over  44  different  organizations  throughout 
Washington,  DC,  Maryland,  Virginia  and  Pennsylvania.  The  Navy  averages  191  blood  drawings  each 
year,  or  16  drawing  per  month,  to  meet  the  hospital's  needs. 


TIME  PERIOD 

TOTAL  #  VISITS 

TOTAL  BAGS  OF  BLOOD 

Jan  1994  -  Dec  1994 

199 

6340 

Jan  1995  -  Dec  1995 

182 

6111 

Jan  1996-  Feb  1996 

27 

892 

Jan  1994  -  Feb  1996 

408 

13343 

Figure  1.  Yearly  sums  for  all  data 


The  mobile  blood  bank  has  an  assigned  team  of  approximately  eight  to  twelve  civilian  and 
military  personnel;  the  team  can  be  augmented  by  hospital  personnel  temporarily  beyond  twelve  for  large 
drawings.  The  blood  drawings  at  the  various  locations  are  scheduled  one  year  in  advance  by  a  member  of 
the  Navy’s  blood  bank.  The  blood  bank  does  not  require  the  organizations  they  visit  to  schedule  donors. 
The  locations  visited  by  the  mobile  blood  bank  have  personnel  that  are  not  always  in  fixed  locations  due  to 
their  jobs  and  schedules.  As  noted  in  the  study  by  Jennifer  Michaels  et  al.  In  "A  Simulation  Study  of 
Donor  Scheduling  Systems  for  the  American  Red  Cross",  1992,  the  Navy's  clients  are  not  conducive  to  a 
scheduling  system.  They  stated  in  their  study  ..."any  company  that  does  not  have  most  of  its  employees  on 
site  during  the  course  of  the  day,  will  benefit  from  a  more  flexible  scheduling  system."1  Therefore,  our 
study  team  did  not  examine  the  impact  of  scheduling  donors  to  improve  the  system’s  effectiveness. 

1  Michaels,  Jennifer,  John  Brennan,  Bruce  Golden  and  Michael  Fu,  "A  Simulation  Study  of 
Donor  Scheduling  Systems  for  the  American  Red  Cross",  Computer  Ops  Res,  Vol.  20,  No  2, 
1992,  pg  212. 
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The  Navy,  as  does  the  American  Red  Cross,  relies  heavily  on  repeat  donors.  Just  as  with  the  Red 
Cross,  the  Navy’s  blood  program  is  voluntary.  Additionally,  the  Navy  has  a  smaller  population  to  draw 
from  for  donors,  but  they  represent  a  healthy  segment  of  the  population.  The  donors  are  from  the  military 
community  (active  duty,  retired,  DoD  civilians  and  family  members)  account  for  less  than  six  percent  of 
the  population.  To  keep  this  small  donor  population  happy  and  returning  the  Navy  was  interested  in  ways 
to  improve  the  overall  blood  donating  system.  They  wanted  to  improve  the  time  required  to  go  through  the 
system,  which  includes  the  total  time  as  well  as  time  in  queues. 

To  simulate  the  process  we  first  had  to  have  an  understanding  of  the  operation  at  the  different 
locations  the  mobile  blood  bank  visited.  Depending  on  the  location,  space  allocated  by  the  organization, 
and  projected  number  of  donors,  these  factors  effected  the  size  of  the  blood  team  and  the  set-up  of  the 
process.  Our  team  went  with  the  blood  bank  to  four  locations  to  observe  the  various  factors  and  collect 
data  on  the  following  dates: 

1.  U.S.  Naval  Academy,  Annapolis,  MD  Nov  27,  1995 

2.  U.S.  Naval  Academy,  Annapolis,  MD  Feb  27,  1996 

3.  National  Naval  Medical  Center,  Bethesda,  MD  Mar  5,  1996 

4.  U.S.  Naval  Academy,  Annapolis,  MD  Mar  26,  1996 

The  set-ups  were  very  similar  at  each  of  the  locations.  The  system  has  approximately  eight 
stations  with  the  number  of  servers  varying  depending  on  the  anticipated  number  of  donors.  Additionally, 
the  first  four  stations  were  sometimes  combined: 

1.  Registration/personal  history. 

2.  Vital  signs. 

3.  Hemoglobin  check. 

4.  Deferral  check  (computer  based  or  hard  copy). 

5.  Interview. 

6.  Bag  issue. 

7.  Phlebotomy  (blood  donating). 

8.  Recovery  (food  and  drinks). 

Station  1 :  Personal  History. 

The  first  station  had  infinite  capacity  because  it  did  not  require  any  servers.  The  amount  of  time 
spent  at  this  first  station  was  dependent  upon  how  fast  the  individual  could  fill  out  the  personnel  history 
and  answer  the  questions.  When  the  potential  donor  has  completed  the  paperwork  he  or  she  moves  to  the 
next  station. 
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Station  2  and  3:  Vital  signs  and  hemoglobin. 

The  vital  signs  check  includes:  temperature,  blood  pressure  and  pulse  check.  The  hemoglobin 
check  is  a  simple  prick  of  the  individuals  finger  followed  by  a  simple  test.  At  both  of  these  stations  a 
potential  donor  could  be  deferred  or  sent  to  station  four.  The  second  and  third  stations  were  difficult  to 
measure  since  the  donors  could  go  to  station  three  before  station  two  or  visa  versa  depending  upon  if  there 
was  a  line  at  one  of  these  two  stations.  Also,  there  were  times  when  a  server  would  combine  the  tasks  of 
station  two  and  three  at  one  location. 

Station  4:  Deferment  check. 

Station  four  consisted  of  one  server  with  a  laptop.  The  server  would  check  the  data  base  to 
determine  if  the  donor  had  been  deferred  from  donating  blood.  Reasons  for  deferral  included  overseas 
assignments,  immunization,  etc.  If  the  donor  is  not  deferred  he  or  she  move  on  to  the  next  station.  This 
station  has  on  occasion  been  combined  with  stations  two,  three  or  five. 

Station  5:  Interview. 

Station  five  is  the  interview  station.  The  sever  or  severs  at  this  station  trained  and  tested  on  the 
military  and  federal  regulations  for  blood  donating.  The  interviewers  re-ask  the  questions  potential  donors 
answered  on  the  personal  history  form  to  include  emphasis  on  sexual  preference,  use  of  drugs,  tattoos, 
travel  outside  the  United  States,  etc.  The  questions  are  asked  in  a  screened  off  area  to  preserve  the 
individuals  privacy.  Once  the  donor  has  successfully  completed  the  interview  he  or  she  goes  to  the  bag 
issue  table,  station  six.  If  donors  fails  this  station  they  depart  the  system. 

Station  6:  Bag  issue. 

At  station  six  the  individual  selects  one  of  two  bar-coded  stickers  which  tells  the  lab  to  use  or  not 
use  their  blood.  From  this  station  the  donor  goes  to  station  seven  to  donate  blood. 

Station  7:  Phlebotomy. 

The  beds  are  set  up  in  groups  of  threes  with  one  phlebotomist  assigned  to  each  group.  The 
number  of  beds  set  up  is  dependent  upon  the  anticipated  number  of  donors.  The  donor  is  sent  to  the  first 
available  bed  to  actually  donate  blood.  The  process  of  filling  the  bag  takes  four  to  seven  minutes  on 
average.  If  the  donor  takes  longer  than  10  minutes  to  fill  the  bag,  the  blood  cannot  be  used  and  that  fact  is 
noted  on  the  donor’s  bag. 

Station  8:  Canteen. 

Upon  completion,  the  donor  then  goes  to  the  canteen  for  snacks  and  beverages.  We  did  not  measure 
station  eight  because  donors  can  control  their  departure  time  once  they  are  feeling  “OK”  to  leave. 
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For  the  blood  donating  process  our  team  assumed  the  donor’s  system  time  started  when  he  or  she 
walked  up  to  the  registration  table  to  fill  out  the  personal  history  form.  The  process  ended  when  the  donor 
left  the  system  because  of  deferment  at  one  of  the  stations  which  we  noted  in  our  time  logs,  or  left  the  bed 
after  donating  blood.  For  this  study  our  team  measured  donor’s  times  at  each  station  in  minutes.  We  did 
not  record  the  seconds  due  to  limited  personnel  to  monitor  the  system.  Our  study  was  similar  to  other 
studies  done  on  the  American  Red  Cross’  mobile  blood  banks  in  terms  of  the  constraints  on  the  system.* 1 2 
Our  problem  had  the  following  constraints:  The  arrival  rate  of  our  donors  was  random  because  the  Navy 
does  not  use  scheduling,  with  arrivals  constrained  to  a  four  hour  window  on  average;  the  system  has  a 
limited  number  of  servers  at  each  station,  due  to  personnel  constraints  and  resource  constraints,  such  as 
equipment  to  take  blood  pressure;  and  because  tests  and/or  questions  completed  at  each  station  can  result  in 
moving  on  to  another  station  or  deferment. 

For  our  base  case  model  we  used  the  following  system  configuration: 


history 

vital  signs 

finger  prick 

deferral 

check 

4  servers 


interview 


2  servers 


bag  issue 


1  server 


9  beds 


3 

servers 


Figure  2.  Base  Case  System 

2 

Brennan,  Hohn,  Bruce  Golden,  and  Harold  Rapport,  "Go  with  the  Flow:  Improving  Read 
Cross  Bloodmobiles  Using  Simulation  Analysis.",  Interfaces  Vol  22,  Sep-Oct  1992,  pp  1-13. 
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II.  MODEL  VALIDATION  AND  VERIFICATION. 


A.  Verification. 

We  modeled  the  Navy’s  system  using  the  computer  simulation  package  GPSS/H  (version  3)  to 
simulate  the  process  as  well  as  Proof  Animation  to  simulate  the  process  independently.  To  verify  our 
programs  we  changed  all  the  service  time  distributions  in  our  model  to  exponential  distributions  with  the 
same  mean  as  our  sample  data.  We  computed  our  theoretical  results  using  Queuing  Analyses  with  TK 
solver  for  windows  to  verify  the  results  generated  by  the  computer  simulation  model  (see  Appendix  A, 
QTK  output).  Queuing  Analyses  with  TK  solver  gave  the  following  results  (where  E(IAT)  is  the  Expected 
inter-arrival  time,  E(ST)  is  the  Expected  service  time,  W  is  Expected  waiting  time  in  the  system,  Wq  is  the 
Expected  waiting  time  in  the  queue,  L  is  the  Expected  system  size,  Lq  it  the  Expected  queue  size,  Pi  is  the 
probability  of  a  server  idle).: 
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Station  3 

2.67 

4.183 

3 

10.6078 

7.6078 

2.5359 

1.8187 

0.7172 

Station  4 

21.1698 

4.183 

21.1698 

21.6302 

0.4604 

5.171 

0.1101 

0.5627 

Figure  3.  QTK  Results. 

To  compare  our  results  with  our  theoretical  results  required  steady  state  of  the  system.  To 
approximate  steady  state  we  ran  the  exponential  computer  simulation  for  forty-eight  hours  and  500 
replications.  The  long  run  length  successfully  overwhelmed  the  system  and  reduce  the  standard  deviation 
to  achieve  steady  state  in  order  get  closer  to  theoretical  solution.  The  absolute  error  between  our  program’s 
results  (GPSS/H)  and  QTK’s  was  minimal,  with  ±  .06  minutes  for  all  point  estimates;  thus  verifying  our 
open  Jackson  network  computer  simulation,  (see  figure  4,  GPSS/H  and  QTK  absolute  error  comparison). 
The  difference  from  GPSS/H  results  and  the  theoretical  is  due  to  the  rounding  error  in  our  data. 
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Station  2 
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Station  3 
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Figure  4.  Absolute  Error  Comparison 


B.  Validation. 

We  were  able  to  validate  the  base  case  model  by  comparing  the  simulated  model  with  the  actual 
system  on  which  we  had  previously  collected  data.  We  compared  the  base  case  simulation  with  the  actual 
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data  from  our  visit  to  the  USNA  on  March  26,  1996.  In  the  base  case,  running  the  model  for  500 
replications  for  the  four  hour  drawing,  the  model  had  71  donors  go  through  the  system.  For  the  actual 
collection  of  data,  there  were  71  donors  go  through  the  system,  three  which  were  discarded  for  lost  data. 
Validating  our  system  allowed  us  to  use  the  base  case  to  study  changes  to  the  number  of  servers  and  inter¬ 
arrival  rates. 

In  addition,  we  had  the  previous  26  months  of  historical  information  on  the  mobile  blood  bank. 
From  this  historical  data  we  knew  the  date  of  the  drawing,  the  number  of  donors  at  each  location,  and  the 
actual  number  deferred.  This  information  was  then  compared  to  the  hospital’s  flat  logs  that  are  maintained 
by  the  blood  lab.  The  flat  logs  register  donors  by  bag  ID  number.  The  flat  logs  record  the  actual  number 
of  bags  processed  from  the  blood  drawings.  They  also  denote  which  bags  of  blood  were  usable  and  which 
ones  were  not.  For  privacy  reasons  we  did  not  record  the  details  of  why  the  blood  was  not  useable.  We 
took  a  random  sample  from  the  flat  logs,  109  sample  dates  all  together.  The  overall  difference  between  the 
flat  log  data  and  the  data  collected  by  the  mobile  blood  bank  personnel  was  less  than  a  five  percent  error 
which  we  considered  acceptable  (se  Appendix  B,  Flatlog  Comparison). 

III.  GENERATION  OF  INPUT. 

We  collected  donor  inter-arrival  times,  service  times,  and  deferments  at  the  different  blood 
drawings  manually  using  synchronized  watches,  (see  Appendix  C  for  the  times).  We  visited  the  Naval 
Academy  a  total  of  three  times.  At  each,  the  drawings  were  similar  in  size  and  configuration.  The  blood 
drawing  at  the  National  Naval  Medical  Center  was  smaller  in  size,  but  of  similar  system  configuration.  All 
four  drawings  were  four  hours  in  duration.  From  the  data  we  were  able  to  determine  the  distributions  of 
the  arrival  rates  and  service  times.  We  used  Unifit  II  to  fit  the  data  to  the  best  distribution  by  matching  the 
first  four  moments  as  closely  as  possible,  (see  appendix  D  for  detailed  Unifit  II  printouts).  Similar  to  the 
study  by  John  E.  Brennan  at  al.,  “Go  with  the  Flow:  Improving  Red  Cross  Bloodmobiles  Using  Simulation 
Analysis”,  1992  we  also  examined  the  arrival  pattern  of  the  donors.  The  drawings  at  the  USNA  were 
scheduled  from  1400  hours  until  1800  hours.  We  determined  that  the  donors  were  following  a  bimodal 
function  as  suspected,  (see  Figure  5  for  Donor  Pattern  charts).  This  is  due  to  the  student’s  schedules,  who 
predominately  get  out  of  class  at  1530  hours  and  leave  activities  an  hour  or  so  later.  The  figure  below 
shows  the  donor  patterns  for  all  three  visits  to  the  USNA  with  the  times  in  minutes. 
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Figure  5.  Donor  Arrival  Patterns. 

In  collecting  the  service  times  we  attempted  to  time  all  the  different  stations.  At  station  one  we 
were  able  to  collect  the  service  time  for  potential  donors  to  complete  the  personnel  history  form  in  minutes. 
However,  we  were  unable  to  collect  the  times  separately  of  the  vitals,  hemoglobin  and  deferral  checks.  At 
these  three  station  donors  were  able  to  go  to  stations  out  of  order  as  they  became  available.  If  we  had  had 
an  automated  system  for  tracking  service  times  then  this  data  would  have  been  available.  Therefore,  we 
dealt  with  only  the  total  service  time  to  complete  stations  one  through  four  and  combined  them  into  one 
station  for  our  model.  At  the  remaining  stations,  (interview,  bag  issue  and  phlebotomy)  we  were  able  to 
collect  data  and  determine  their  distributions  for  the  computer  simulation  model,  (see  Figure  6,  Station 
statistics  for  Base  Case  Model). 
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Figure  6.  Station  Statistics  for  Base  Case  Model. 

For  the  interviewers,  station  two  in  our  model,  we  collected  the  service  times  for  each  of  the 
interviewers  to  include  break  times.  When  an  interviewer  went  on  break  this  caused  the  queue  to  build  up 
in  front  of  one  server.  We  were  able  to  model  a  server  on  break  in  GPSS/H  to  note  the  overall  impact  on 
the  system.  For  this  station,  the  service  times  of  the  two  interviewers  was  very  similar,  the  differences 
were  negligible.  Therefore,  for  simplicity  within  our  model  we  assumed  the  same  service  time  for  each 
interviewer.  However,  for  station  three  in  the  model,  bag  issue,  the  service  times  were  short.  The  service 
times  were  one  of  three  times:  two,  three  or  four  minutes  because  we  had  truncated  the  times  to  minutes, 
not  seconds.  For  this  station  we  modeled  it  as  a  discrete  distribution. 
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At  the  phlebotomy,  station  four  in  the  model,  we  did  not  look  at  the  service  times  of  the  individual 
blood  bank  employees.  Our  data  at  this  station  combines  into  one  service  time  the  various  stages  the  donor 
goes  through:  Preparation  of  the  arm,  actual  filling  of  the  bag,  and  removal  of  the  needle.  Because  we 
calculated  the  service  time  as  the  minutes  between  time  the  donor  went  to  a  bed  to  the  time  he  or  she  left 
the  bed,  we  assumed  each  bed  was  a  server,  not  the  attendant/phlebotomist  of  three  beds. 

IV.  THE  EXPERIMENTAL  DESIGN: 

A.  Computer  Simulation  and  Model  of  the  Navy’s  system. 

The  Navy’s  mobile  blood  bank  donor  process  is  an  FCFS  open  Jackson  network  with  donors 
arriving  randomly  into  station  one.  In  our  base  case  computer  simulation  model  donors  arrived  into  station 
one  (personal  history  paperwork,  vital  signs,  hemoglobin  and  deferral  check)  with  98  percent  move  on  to 
station  two  (interview).  At  station  two  17  percent  are  deferred  with  83  percent  moving  to  station  three  (bag 
issue).  All  donors  in  our  system  moved  from  three  to  station  four  (phlebotomy)  and  then  out  of  the  system. 

We  were  interested  in  the  effects  on  the  system  if  the  number  of  servers  were  varied.  By  varying 
the  system  to  improve  one  area,  such  as  the  interview  station,  we  did  not  want  to  create  queue  build  ups  at 
other  stations,  like  the  beds.  We  ran  four  variations  on  our  simulation  model  of  the  Navy’s  process.  They 
included: 

1 .  Exponential  Case.  Used  to  determine  steady  state  to  verify  our  model. 

2.  Base  Case.  Modeled  the  real  system  we  observed.  Interviewers  at  Station  two  were  on  a 
clock  to  put  one  interviewer  on  break  twenty  minutes  every  other  hour,  (90  to  1 10  minutes 
and  210  to  240  minutes).  Changes  were  made  to  this  model  to  note  improvements  to  the 
overall  system. 

3.  Two  Interviewers.  Used  base  case  but  with  two  interviews  in  place  at  all  times. 

4.  Decreased  the  number  of  beds.  Changed  the  number  of  beds  from  nine  to  six. 

For  each  case  the  number  of  servers  at  each  station  was  as  follows: 


CASE 

STATION  1 

STATION  2 

STATION  3 

STATION  4  | 

Expon  Case 

4 

2 

1 

9 

Base  Case 
(interviewers  on 
breaks) 

4 

1-2 

1 

9 

Two  interviewers 
always 

4 

2 

1 

9 

Decrease  beds:  6 

4 

2 

1 

6 

Figure  7.  Number  of  Servers  at  Each  Station. 
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To  show  the  difference  from  the  actual  observed  system  and  our  simulation  model  we 
used  animation.  We  first  animated  our  visit  to  the  Navy  Academy  on  the  27th  of  February.  This 
system  has  the  same  number  of  servers  at  each  of  the  stations  as  the  Base  case.  However,  the 
animation,  unlike  the  Base  case  model  shows  the  true  arrivals  of  the  donors  to  include  batch 
arrivals  and  the  exact  times  the  interviewers  went  on  break,  reducing  that  station  to  one  server. 
The  second  animation  shows  the  base  case  with  two  servers  dedicated  to  station  two  and  the 
distributions  for  the  inter-arrival  times  and  service  times  as  determined  by  Unifit  II.  The  third 
animation  is  the  same  as  the  second  but  with  the  number  of  beds  reduced  to  six.  All  three 
animation  models  visually  show  how  the  changes  impact  the  system.  Most  notable  was  the 
development  of  queues  at  various  stations.  The  animation  also  show  deferments  being  rejected  at 
the  various  decision  points  in  the  system.  The  deferments  collected  at  the  bottom  of  the  model  to 
show  the  cumulative  sum  of  deferrals  in  a  four  hour  drawing. 

B.  Trend  Analysis. 

The  Navy  mobile  blood  bank  has  maintained  logs  for  the  last  2  years  of  its  operation  (we  have  the 
data  from  January  1994  -  February  1996).  These  logs  include  information  regarding:  (1)  Visit  location,  (2) 
day,  date,  year  of  visit,  (3)  anticipated  donor  turnout,  (4)  actual  donor  turnout,  (5)  amount  of  bags  obtained 
for  the  given  day,  and  (6)  number  of  people  deferred.  (See  Appendix  E )  We  used  this  data  to  look  for 
possible  trends  in  donor  turnout.  We  quickly  realized  that  two  years  of  data  is  definitely  a  minimum 
amount  necessary  for  any  trend  analysis.  However,  we  were  able  to  use  the  results  from  different  statistical 
tests  as  indicators  of  potential  trends.  These  indicators  provide  insight  to  areas  of  potential  interest  and 
those  which  should  be  tracked  in  coming  years.  We  analyzed  the  available  data  using  the  following: 

1.  Runs  tests 

-  Runs  Up/Down  test 

-  Run  length 

2.  Linear  trends  test. 

3.  Additive  time  series  model. 

4.  Auto  regressive/moving  averages. 

C.  Runs  Tests. 

We  implemented  two  types  of  Runs  test  as  a  diagnostic  procedure  in  part  to  check  the 
reasonableness  of  the  assumption  that  our  data  is  a  sequence  of  binary  outcomes  from  independent  and 
identically  distributed  (iid)  Bernoulli  random  variables.  The  first  test  is  a  Runs  Up  /  Down  test,  and  the 
second  is  a  test  taking  into  account  the  length  of  each  run.  The  question  we  were  interested  in  answering 
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was:  ‘is  our  sequence  of  time  series  data  occurring  by  random  chance,  or  is  there  evidence  indicating  a  lack 
of  randomness  in  the  ordering  of  the  data?’ 

1.  Runs  Up  /  Down  Test. 

Initially  we  considered  the  entire  sequence  of  bi-monthly  sums  using  a  Runs  Up  /  Down  test  to 
identity  patterns  in  our  time  series  data  which  are  unlikely  to  occur  if  the  iid  Bernoulli  random  variable 
assumption  is  valid.  To  assess  whether  or  not  our  observed  sequence  of  outcomes  is  incompatible  with  an 
assumption  of  randomness  we  compared  the  observed  number  of  up  /  down  runs  with  the  number  which  is 
expected  if  all  possible  orderings  of  n0  ups  and  downs  are  equally  likely. 

A  test  of  this  nature  will  give  us  some  indication  of  whether  changes  in  our  sequence  is  a 
departure  from  randomness  and  indicative  of  a  persistence  in  its  direction  of  movement  or  where  our 
sequence  contains  a  trend  (e.g.  a  cyclical  pattern).  Our  hypothesis  and  test  statistics  were: 

Hypothesis: 

H0:  Sequence  generated  by  a  random  process 

Hj:  Sequence  generated  by  a  process  containing  either  persistence  or  frequent  changes  in 

direction. 


Test  Statistics: 

E(R)  =  (2n  -  1  )/3  a2(R)  =  ( 1 6n  -  29)  /  90  z*  =  (R  -  E(R))  /  s(R) 


Bi-monthly  Sums: 


Sum(+/-)* 

Run# 

Sum(+/-) 

Run  # 

Sum(47-) 

Run  # 

252 

228- 

273+ 

18 

256- 

27 

173- 

l 

216- 

11 

241- 

19 

277+ 

311+ 

2 

259+ 

338+ 

321  + 

28 

302- 

3 

278+ 

12 

528+ 

20 

234- 

29 

378+ 

4 

265- 

13 

212- 

293+ 

30 

163- 

5 

286+ 

209- 

21 

127- 

31 

251+ 

6 

288+ 

14 

243+ 

22 

461  + 

32 

7 

276- 

199- 

23 

145- 

8 

239- 

15 

286+ 

24 

104- 

181- 

9 

!  300+ 

16 

34- 

33 

227+ 

216- 

156- 

25 

395+ 

34 

373+ 

10 

178- 

17 

235+ 

232- 

353- 

261+ 

291+ 

26 

I  231- 

35 

*+/-  indicates  up  /  down  run 


Our  Results: 

E(R)=  34.33  a2(R)  =  8.922  z^*  =  0.2243 

At  an  a  =  0.05;  -1 .96  >=  |0.2243|  <=  1 .96,  p-value  =  .956 
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We  also  ran  a  Runs  Up  /  Down  test  on  the  correlation  coefficient  obtained  when  comparing  the 
same  point  in  time  for  the  year  1994  and  1994  (e.g.  week  2  of  February  1994  vs.  week  2  of  February 
1995).  We  used  the  same  hypothesis  to  test  our  results. 


Correlation  Coefficients: 


Sum(+A) 

Run  # 

Sum(+/«) 

Run# 

Sum(+/-) 

Run  # 

Sum(+/~) 

Run  # 

-0.065728 

-0.066126 

7 

-0.264036 

11 

-0.124466 

0.035185 

0.021972 

8 

-0.008851 

12 

-0.188667 

l 

0.036203 

4 

-0.061436 

9 

-0.026086 

0.197515 

2 

5 

■mm— 

13 

-0.166030 

3 

0.095876 

6 

0.096705 

14 

0.075071 

-0.050343 

-0.024776 

15 

Our  Results: 

E(R)=  15.667  ct(R)=  1.986 

At  an  a  =  0.05;  -1.96  >=  |0.3359|  <=  1.96 

zM*=  -0.3359 

These  test  results  indicate  the  sequence  is  random;  it  is  appropriate  to  treat  the  observations  as  a 
random  sample  from  an  infinite  population.  However,  according  to  resident  GMU  statistics  expert,  Dr. 
Sutton,  Runs  can  be  quite  ineffective  for  detecting  inconsistency  in  variation  if  the  variation  is,  for 
example,  cyclical  and  the  period  is  not  very  long.  For  our  time  series  data,  we  only  have  two  years  worth 
of  information;  enough  to  make  some  initial  observations,  but  not  enough  information  to  rule  out  possible 
trends.  With  this  in  mind,  we  continued  with  other  methods  of  time  series  analysis. 


2.  Run  Length. 

Performing  a  test  to  see  if  the  longest  run  warrants  anything  other  than  iid  also  yielded 
insignificant  results.  The  longest  run  is  three  and  the  p-value  associated  with  that  amount  is  0.9466.  This 
indicates  there  is  no  reason  to  reject  the  assumption  that  the  data  comes  from  anything  other  than  iid. 

D.  Linear  Trends. 
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Analyzing  the  linear  trends  yielded  very  interesting  results.  We  first  looked  at  the  linear  trend  of  the  data 
across  the  26  months  (See  Figure  8).  The  linear  trend  for  the  26  month  span  is:  Yt  =  274,923  -  2,85162*t , 
indicating  a  decrease  of  approximately  74  bags  over  the  last  two  years  (-2.85162*26).  Initially,  this 
decrease  did  not  appear  to  be  extremely  significant.  Yet,  there  are  additional  studies3  that  claim  an  overall 
decline  in  donor  participation  is  occurring.  We  felt  that  our  declining  trend  together  with  the  studies 
indicating  decline  warranted  a  deeper  look  into  the  data. 


Monthly  Sums  for  All  Data 

Linear  Trend  Model 
Yt  =  55 1.689- 2.85 162*t 


Time 


•  Actual 
A  Fits 

—  Actual 

-  -  Fits 


Figure  8.  Monthly  Sums. 


Looking  at  1994  and  1995  separately  revealed  greater  insight  to  the  source  of  declining 
participation.  The  linear  equation  describing  1994  is:  Yt=  514,651  +  2.10490*t  (See  Figure  9).  This 
equation  indicates  an  increase  of  approximately  25  bags  over  the  year.  Still  not  very  significant  since  we 
are  only  dealing  with  twelve  data  points. 


3  Roberts,  Russell  and  Michael  Wolkoff,  "Improving  the  Quality  and  Quantity  of  Whole  Blood 
Supply:  Limits  to  Voluntary  Arrangements",  Journal  of  Health  Politics,  Policy  and  Law, 
1988,  Vol.  13,  No.  1,  ppl  67-1  77. 
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1994  Monthly  Sums  for  All  Data 

Linear  Trend  Model 
Yt  =  514.651  +  2.10490*t 


Actual 

Fits 

Actual 

Fits 


Figure  9.  1994  Monthly  Sums 


The  linear  equation  for  1995  is:  Yt  =  57L409  -  9.56294H  (See  Figure  10).  This  equation  suggests 
a  loss  of  approximately  115  bags  over  the  year. 


1995  Monthly  Sums  for  All  Data 

Unear  Trend  Model 
Yt  =  571.409  -9.56294*1 

900  _ 

800  _ 

700  _ 

600  - 
500  - 
400  - 
300  _ 

200  - 

1  I  i 

0  5  10 

Time 

Figure  10.  1995  Monthly  Sums. 

The  linear  model  representing  1995  lends  more  support  to  the  notion  of  declining  donor 
participation.  It  is  interesting  to  note  the  average  blood  drawn  for  1995  is  only  slightly  less  than  the 
average  for  1994  (509  vs.  528),  yet  there  is  a  much  more  drastic  negative  trend  in  1995.  Again,  our  data  is 
limited  to  only  two  years,  but  we  can  conclude  areas  which  may  be  potential  indicators  of  actual  trends. 

To  determine  if  the  Navy  Mobile  Blood  Bank  should  be  concerned  with  the  1995  trend  and  focus 
on  methods  to  minimize  the  affects  of  a  declining  donor  population  we  subjected  the  data  to  a  number  of 
additional  statistical  tests,  (see  Appendix  F  for  all  additional  test  results).  Parametric  and  nonparametric 
tests  for  rejecting  iid  based  on  differences  in  the  years  yielded  extremely  insignificant  p-values  (t-test,  sign 
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test,  Wilcoxon  test,  Mann-Whitney).  These  results  imply  that  the  data  is  not  really  following  any  type  of 
trend  even  though  visually  there  appears  to  be  something  occurring. 

E.  Additive  Time  Series  Model. 

Running  an  Additive  Time  Series  Model  in  Minitab  enabled  us  to  look  at  our  data  sequence  with  the 
Seasonal  Component  isolated  (See  Figure  1 1).  This  type  of  model  is  looking  at  seasonal  trends  together 
with  some  type  of  trend  component  (linear  or  exponential  obtained  using  a  least  squares  calculation)  and 
cyclical  component  (deviations  from  the  trend). 

The  trend  obtained  from  this  model  exactly  matched  the  linear  trend  model  discussed  in  the 
previous  section.  For  seasonal  indicators,  we  see  evidence  of  extreme  variation  over  the  seasonal  periods 
(approximately  10  to  90  bags  of  variation).  The  high  variation  and  seasonal  indices  for  March  can  be 
explained  by  the  outlier  data  point  that  is  due  to  a  rare  day  at  the  Naval  Academy  when  253  bags  of  blood 
were  drawn.  Furthermore,  there  does  not  appear  to  be  any  cyclical  trend. 


Additive  Model  for  All  Data  (Monthly  Sums) 


Seasonal  Indices 


r 

1  2  3  4  5  6  7  8  9  10  11  12 

Percent  Variation,  by  Seasonal  Period 

ii 

...lull 

1  2  3  4  5  6  7  8  9  10  11  12 


900 
000 
700 
600 
500 
400 
300 
200 

1  2  3  4  5  6  7  8  9  10  11  12 

150 
50 
-50 
-150 
-250 
4350 

1  2  3  4  5  6  7  8  9  10  11  12 


Original  Data,  by  Seasonal  Period 


Figure  11.  Additive  Model. 


There  is  not  a  lot  to  obtain  from  an  additive  model  with  only  two  years  of  data.  With  two  years  of 
data,  the  model  really  only  has  two  data  points  to  compare  (one  from  1 994,  and  the  other  from  the  same 
period  in  1995).  It  would  be  interesting  to  see  how  data  over  the  next  two  to  three  years  supports  the 
indicators  represented  in  the  figure  above. 

Simply  eyeing  the  raw  data  seems  to  support  the  idea  that  there  are  seasonal  factors  affecting  the 
amount  drawn.  January’s  average  is  consistently  below  the  norm  (43 1  vs.  513),  as  well  as  the  months  of 
April  and  May  (460  and  450  respectively).  September  through  November  tend  to  be  above  average 
months  (570,  550,  and  551  respectively).  These  observations  are  supported  by  the  Seasonal  Indices  chart, 
but  until  more  data  can  be  obtained,  not  a  lot  of  confidence  can  be  offered  from  the  additive  model.  In 
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fact,  a  Durbin-Watson  test  statistic  of  2.15  indicates  there  is  no  support  of  positive/negative  correlation  of 
error  terms.  This  suggests  there  is  no  need  to  even  look  at  time  series  analysis. 

F.  Autoregressive  /  Moving  Averages. 

Autoregressive  /  Moving  Averages  (ARMA)  is  a  form  of  analysis  which  generates  a  model  using  white 
noise  as  the  forcing  terms  in  a  set  of  linear  differences  equations.  This  is  an  iterative  process;  we  had  to  try 
different  combinations  of  AR  and  MA  types  (e.g.  AR(2)  and  MA(2))  to  obtain  a  model  best  fitting  our 
data.  Using  Minitab  to  calculate  the  ARMA  statistic,  we  found  that  an  AR(1)  MA(1)  best  describes  our 
data.  The  Minitab  results  are  as  follows: 


Minitab  Output 

Final  Estimates  of  Parameters 

Type 

Estimate 

AR  1 

0.8449 

MA  1 

1.0052 

Constant 

40.0585 

Mean 

258.225 

To  measure  how  well  the  model  fits  the  data,  we  used  the  Minitab  output  for  the  modified  Box- 
Pierce  chi-square  statistic.  We  computed  the  p-value  to  check  the  significance  of  the  value  for  each  of  the 
lags  (12,  24,36). 

H0:  The  specified  ARMA  model  fits  our  data. 

Hp  The  specified  ARMA  model  does  not  fit  our  data. 


Modified  Box-Pierce  (Ljung-Box)  chi-square  statistic 


Lag 

12 

24 

36 

Chi-square 

6.8(DF=10) 

17.5(DF  =  22) 

25.7(DF=34) 

p-value: 

0.744159 

0..  735185 

0.846087 

We  accept  the  null  hypothesis  based  on  the  observed  p-value  and  are  unable  to  conclude  that  the 
model  obtained  from  Minitab  does  not  fits  our  data.  The  model  is: 


Drawn  at  time  t  =  40.0585  +  0.8449(Drawnt4)  +  Zt  +  1.0052(2^) 

(where  Zt  is  the  “ white  noise  error  ”  ~N(0,o2)). 


Obtaining  a  model  with  such  an  insignificant  p-value  for  rejecting  the  Minitab  model,  suggests  the 
ARMA  process  is  potentially  a  good  method  of  forecasting  coming  months.  However,  we  must  keep  in 
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mind  we  have  only  two  years  of  data  and  the  previous  Durbin-Watson  test  (together  with  several  additional 
tests)  indicates  an  iid  process.  Obtaining  two  to  three  more  years  of  data  would  yield  much  more 
significant  results.  Hence,  we  interpreted  the  information  obtained  from  the  ARMA  process  as  offering  an 
interesting  suggestion  of  the  potential  for  using  the  above  model  for  forecasting  the  data. 

V.  RUN  SUMMARIES. 


The  summary  of  our  runs  output  is  in  Appendix  G.  We  ran  the  various  cases  for  20  runs  for  four 
hours,  500  runs  for  four  hours. 

VI.  ANALYSIS  OF  OUTPUT. 


The  following  adjustments  to  the  system  were  tested  at  500  replications  for  four  hours: 


SYSTEM 

AVG  TIME  IN  SYSTEM 

TOTAL  DONORS 

a.  Exponential  case 

52.09 

68.30 

b.  Base  case  (2  inter-  j 

viewer-  take  breaks) 

38.76 

70.60 

1  c.  Two  interviewers 

(always) 

36.99 

70.60 

d.  Decrease  beds  6 

37.04 

69.80 

Figure  12.  Model  Output  Run  Summary. 


By  changing  the  number  of  servers  at  station  two  (interviewer)  and  at  station  four  (blood  drawing) 
we  noted  slight  improvement  in  the  total  time  through  the  system.  Decreasing  the  number  of  beds  did  not 
change  the  overall  time  in  the  system.  A  similar  decrease  in  total  system  time  was  noted,  as  compared  to 
the  base  case,  when  there  were  two  dedicated  interviewers.  In  both  changes  the  queues  actually  decreased 
(see  Figure  13  ,  Summary  of  Model  Variations).  The  reduction  in  queue  build  up  at  the  interview  station 
was  most  noticeable  in  the  animation  simulation  in  comparing  the  27  February  actual  system  to  the 
modeled  system.  The  February  27  model  has  queues  form  when  one  interviewer  goes  on  break  which  is 
visually  demonstrated  in  the  animation.  The  minimal  difference  between  the  simulated  models  is  due  to 
the  lack  of  any  large  queue  forming  as  noted  in  our  output.  The  Lq  never  builds  up  in  the  Base  case  as  was 
actually  observed  at  the  various  blood  drawing  because  we  were  not  able  to  program  for  batch  arrivals  into 
our  model. 
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CASE 

L  (system  size) 

W  (waiting  time 
in  system) 

Lq  (queue  size) 

Wq  (waiting  time 
in  queue) 

a.  Expon  Case 

13.03 

52.09 

2.27 

9.51 

b.  Base  Case 

8.59 

37.02 

0.44 

1.92 

c.  2  interviewers 

8.55 

36.99 

0.44 

1.95 

d.  Decrease  beds 

8.55 

37.04 

0.44 

1.95 

Figure  13.  Summary  of  Model  Variations. 


VII.  CONCLUSIONS. 

From  the  beginning  of  our  project  we  attempted  to  produce  three  products  to  assist  the  Navy  in 
improving  their  number  of  donors  and  overall  donor  satisfaction.  We  examined  the  Navy’s  larger  blood 
drawings  of  70  to  100  donors.  From  our  study  we  noted  that  the  system  could  be  improved  by  reducing 
the  time  donors  spend  the  system  as  a  whole.  The  total  system  time  in  our  model  was  reduced  by  having 
two  dedicated  interviewers  at  station  two. 

There  is  not  a  significant  reduction  in  system  time  between  the  two  models  because  our  model  did 
not  demonstrate  the  queue  build  up  we  observed  at  the  USNA  blood  drawings.  We  were  unable  to 
program  batched  arrivals  into  the  system  as  we  observed  at  the  drawings.  Batched  arrivals  would  have 
reflected  the  true  system  in  terms  of  queue  build  up  at  the  various  stations.  In  the  animation  model  of 
February  27th  all  of  the  arrivals  times  were  entered  as  discrete  times  which  is  a  more  accurate  reflection  of 
the  system  we  observed  because  it  allows  manual  batching  of  the  donors.  However,  from  our  survey  of 
donors  we  noted  that  station  two  was  a  point  of  dissatisfaction. 

In  the  survey  we  randomly  surveyed  donors  at  the  blood  drawing  at  the  Naval  Academy  on  March 
26th.  We  wanted  to  determine  the  satisfaction  or  dissatisfaction,  of  the  donors  and  their  time  in  system. 

The  survey  included  the  start  and  end  time  in  the  system  (not  including  the  canteen),  and  the  following 
questions: 

-  Was  the  process  too  long? 

-  Is  there  anything  in  the  process  you  would  change? 

-  Would  you  donate  again? 

The  survey  confirmed  our  initial  thoughts.  Out  of  the  33  randomly  selected,  two  were  deferred 
and  3 1  completed  the  system.  Of  the  thirty-one,  seventy-four  percent  of  those  surveyed  concluded  the 
system  time  was  not  too  long,  forty-eight  percent  recommended  changes,  and  one-hundred  percent  said 
that  they  would  donate  again  in  the  future.  On  average,  it  took  those  surveyed  fifty-eight  minutes  to 
complete  the  system.  Those  surveyed  who  did  offer  comments  focused  on  the  interview  station  and 
recommended  more  servers. 

The  students  surveyed  support  our  conclusions  on  the  need  to  maintain  two  interviewers  at  station 
two  at  all  times.  We  recommend  the  supervisor  act  as  a  relief  person  when  assigned  interviewers  need  a 
break.  If  the  supervisor  serves  also  as  an  interviewer  then  the  system’s  queues  at  that  station  develop  as 
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noted  in  our  base  case  model.  At  previous  drawings  we  had  observed  there  were  two  interviewers, 
however,  one  of  the  interviewers  fulfilled  two  jobs:  interviewer  and  overall  supervisor. 

In  addition  to  the  number  of  servers  at  station  two,  we  also  concluded  that  the  number  of  beds 
could  be  decreased  from  nine  to  six  for  drives  with  donor  populations  less  than  or  equal  to  70.  Previously 
the  Navy  had  listed  the  following  table  in  its  operations  manual  as  a  guideline: 


Number  of  Donors 

Number  of  Beds 

Donors  per  Hour 

<30 

4  to  5 

16  to  20 

30  to  45 

6  to  8 

24  to  32 

46  to  60 

9  to  10 

36  to  40 

61  to  80 

1 1  to  12 

44  to  48 

>80 

>  12 

>48 

Figure  14.  Donors  to  beds. 

To  draw  any  more  detailed  conclusions  on  the  number  of  servers  we  would  have  to  implement  our 
proposed  changes  and  observe  their  impact  on  the  system.  However,  the  Navy,  to  conserve  its  resources 
would  benefit  by  further  studying  the  ratio  of  servers  to  donors. 

An  additional  area  of  study  should  be  reducing  the  number  of  decision  points.  Within  the  current 
system  there  is  redundancy  in  the  questions  the  donors  are  asked  from  station  one  to  station  two.  Station 
two  re-asks  the  donor  the  questions  on  the  form.  A  reexamination  of  what  questions  really  need  to  be 
asked  at  each  station  of  the  donors  could  also  reduce  time  within  the  system. 

Finally,  we  noted  in  our  trend  analysis  that  donor  participation  is  decreasing.  But  with  only  26 
months  of  data  we  cannot  draw  any  strong  conclusions.  All  the  tests  performed  show  the  data  to  be  iid. 
Some  of  the  tests  offered  possibilities  of  indicating  potential  trends  but  we  cannot  extracted  anything 
significant  from  the  trend  analysis  without  more  data.  The  data  will  provide  insights  and  should  be  further 
collected  and  monitored  in  the  coming  years. 

Overall  the  Navy’s  system  is  adequate  in  its  current  state.  But  we  know  from  our 
observations,  donor  surveys  and  prior  studies,  that  the  Navy  should  do  what  it  can  to  continuously  improve 
donor  satisfaction  and  examine  other  ways  to  obtain  new  donors. 
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'Input 

Name 

Output 

Unit 

Comment 

M/M/c:Multiple  Servers/Unlimited  Queue 

iat 

3.4843 

min 

Mean  interarrival  time 

St 

14.1443 

min 

Mean  time  to  complete  service 

.287 

lambda 

1/min 

Arrival  rate  (arrivals/unit  of  time) 

.0707 

mu 

1/min 

Service  rate  per  channel  (#/time) 

r 

4.0594 

Avg  #  arrivals  during  avg  service  time 

12 

c 

#  of  servers  in  the  system  (c  >  1) 

rho 

.3383 

Fraction  of  time  each  server  is  busy 

pO 

.0173 

Probability  of  0  in  the  system 

6 

n 

Target  #  of  customers  in  the  system 

pn 

.1073 

Probability  of  n  in  the  system 

Lq 

.0006 

Expected  queue  size 

L 

4.06 

Expected  system  size 

Wq 

.0019 

min 

Expected  waiting  time  in  the  queue 

W 

14.1462 

min 

Expected  waiting  time  in  the  system 

10 

t 

min 

Specific  time  in  the  queue 

Ptq 

0 

Prob.  of  waiting  >=  t  in  the  queue 

PqO 

.9989 

Probability  of  no  wait  in  the  queue 

10 

K 

Max  variable  value  whose  prob  wanted 

pK 

.0058 

Probability  of  K  in  system  (K>=c) 

PK 

.9989 

Probability  of  <=  K  in  system 

1 

d 

min 

Size  of  time  interval  for  plot 

60 

T 

min 

Total  time  horizon  for  prob  plotting 

TWq 

1 

Probability  that  queue  delay  <=  T 

(should  be  1  if  full  plot  is  needed) 

See  Tables  and  Plots 


Variable  Sheet 
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Input 

Name 

Output 

Unit 

4.0992 

iat 

min 

4.4769 

St 

min 

lambda 

.244 

1/min 

mu 

.2234 

1/min 

r 

1.0921 

2 

c 

rho 

.5461 

pO 

.2936 

6 

n 

pn 

.0156 

Lq 

.464 

L 

1.5562 

Wq 

1.9022 

min 

W 

6.3791 

min 

10 

t 

min 

Ptq 

.0508 

PqO 

.6143 

10 

K 

pK 

.0014 

PK 

.9983 

1 

d 

min 

60 

T 

min 

TWq 

1 

Comment 

M/M/c:Multiple  Servers/Unlimited  Queue 

Mean  interarrival  time 

Mean  time  to  complete  service 

Arrival  rate  (arrivals/unit  of  time) 

Service  rate  per  channel  (#/time) 

Avg  #  arrivals  during  avg  service  time 
#  of  servers  in  the  system  (c  >  1) 

Fraction  of  time  each  server  is  busy 

Probability  of  0  in  the  system 
Target  #  of  customers  in  the  system 
Probability  of  n  in  the  system 

Expected  queue  size 
Expected  system  size 
Expected  waiting  time  in  the  queue 
Expected  waiting  time  in  the  system 

Specific  time  in  the  queue 
Prob.  of  waiting  >=  t  in  the  queue 
Probability  of  no  wait  in  the  queue 

Max  variable  value  whose  prob  wanted 
Probability  of  K  in  system  (K>=c) 
Probability  of  <=  K  in  system 

Size  of  time  interval  for  plot 
Total  time  horizon  for  prob  plotting 
Probability  that  queue  delay  <=  T 
(should  be  1  if  full  plot  is  needed) 

See  Tables  and  Plots 


Variable  Sheet 


Untitled 
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Input 

Name 

Output 

Unit 

4.183 

iat 

min 

3 

St 

min 

lambda 

.2391 

1/min 

mu 

.3333 

1/min 

r 

.7172 

1 

c 

rho 

.7172 

pO 

.2828 

6 

n 

pn 

.0385 

Lq 

1.8187 

L 

2.5359 

Wq 

7.6078 

min 

W 

10.6078 

min 

10 

t 

min 

Ptq 

.2794 

PqO 

.2828 

10 

K 

pK 

.0102 

PK 

.9742 

1 

d 

min 

60 

T 

min 

TWq 

.9975 

Comment 

M/M/c:Multiple  Servers/Unlimited  Queue 

Mean  interarrival  time 

Mean  time  to  complete  service 

Arrival  rate  (arrivals/unit  of  time) 

Service  rate  per  channel  (#/time) 

Avg  #  arrivals  during  avg  service  time 
#  of  servers  in  the  system  (c  >  1) 

Fraction  of  time  each  server  is  busy 

Probability  of  0  in  the  system 
Target  #  of  customers  in  the  system 
Probability  of  n  in  the  system 

Expected  queue  size 
Expected  system  size 
Expected  waiting  time  in  the  queue 
Expected  waiting  time  in  the  system 

Specific  time  in  the  queue 
Prob.  of  waiting  >=  t  in  the  queue 
Probability  of  no  wait  in  the  queue 

Max  variable  value  whose  prob  wanted 
Probability  of  K  in  system  (K>=c) 
Probability  of  <=  K  in  system 

Size  of  time  interval  for  plot 
Total  time  horizon  for  prob  plotting 
Probability  that  queue  delay  <=  T 
(should  be  1  if  full  plot  is  needed) 

See  Tables  and  Plots 


Variable  Sheet 


Untitled 


Page 


Input 

Name 

Output 

Unit 

4.183 

iat 

min 

21.1698 

St 

min 

lambda 

.2391 

1/min 

mu 

.0472 

1/min 

r 

5.0609 

9 

c 

rho 

.5623 

pO 

.0062 

6 

n 

pn 

.1458 

Lq 

.1101 

L 

5.171 

Wq 

.4604 

min 

W 

21.6302 

min 

10 

t 

min 

Ptq 

.0133 

PqO 

.9143 

10 

K 

pK 

.0211 

PK 

.9729 

1 

d 

min 

60 

T 

min 

TWq 

1 

Comment 

M/M/c:Multiple  Servers/Unlimited  Queue 

Mean  interarrival  time 

Mean  time  to  complete  service 

Arrival  rate  (arrivals/unit  of  time) 

Service  rate  per  channel  (#/time) 

Avg  #  arrivals  during  avg  service  time 
#  of  servers  in  the  system  (c  >  1) 

Fraction  of  time  each  server  is  busy 

Probability  of  0  in  the  system 
Target  #  of  customers  in  the  system 
Probability  of  n  in  the  system 

Expected  queue  size 
Expected  system  size 
Expected  waiting  time  in  the  queue 
Expected  waiting  time  in  the  system 

Specific  time  in  the  queue 
Prob.  of  waiting  >=  t  in  the  queue 
Probability  of  no  wait  in  the  queue 

Max  variable  value  whose  prob  wanted 
Probability  of  K  in  system  (K>=c) 
Probability  of  <=  K  in  system 

Size  of  time  interval  for  plot 
Total  time  horizon  for  prob  plotting 
Probability  that  queue  delay  <=  T 
(should  be  1  if  full  plot  is  needed) 

See  Tables  and  Plots 
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_ ^ _  Interview _ _ _ _ Bag _ _ _ Take  Blood _ _ Totals 

lnt#|  In  Out  I  Time!  Mini  Def.  In  I  Out  1  Tlmel  Min  Bed#!  In  Out!  Time!  Min  Served  |lni 

1  14:19  14:22  0:03  3  14:22  14:26  0:04  4  2  14:26  14:40  0:14  14  35 
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_ _ i _  Interview _ 

lnt#|  In  |  Out  Tlmei  Min 

1  17:29  17:33  0:04  4 
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Specified  Models  and  Their  Parameters  Sample:  Data  From  A_ARRST1.TXT 

eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee 
Model  1:  Gamma  Distribution 

Location  Parameter  .15789 

Scale  Parameter  6.31422 

Shape  Parameter  .53458 


Quantile  Estimate 
M.L.  Estimate 
M.L.  Estimate 


Model  2:  Weibull  Distribution 
Location  Parameter 
Scale  Parameter 
Shape  Parameter 


0. 

3.19219 

.82883 


Default 
M.L.  Estimate 
M.L.  Estimate 


Model  3 :  Weibull  Distribution 

Location  Parameter  .15789 

Scale  Parameter  2.61592 

Shape  Parameter  .66020 

Model  4:  Gamma  Distribution 

Location  Parameter  0. 

Scale  Parameter  4.64850 

Shape  Parameter  .76010 


Quantile  Estimate 
M.L.  Estimate 
M.L.  Estimate 


Default 
M.L.  Estimate 
M.L.  Estimate 


eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee 
Screen  1/3  -  Press  Fl-4  for  help  or  another  allowed  key:  Logfile:  Open,  On 


UniFit  Manual  Model  Selection  Sample:  Data  From  A_ARRST1.TXT  • 

eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee 

Oaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa<£ 

0  Functional  Groups  (Phase)  0 

LaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaQ 
0  Sample  maintenance  ° 

0  Descriptive  sample  summaries  (I)  ° 

0  Model  specification  (II)  0 

0  Goodness-of-fit  assessment  (III)  0 
0  Inferences  about  model  and  sample  0 
0  Change  to  guided  selection  mode  0 
°  exit  manual  selection  mode  ° 

aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaai 


eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee 
Fl-4  =  Help,  F8-10  =  Logfile,  ESC  =  Exit  Logfile:  Open,  On 


Model  Moment  Comparison  Sample:  Data  From  A_ARRST1.TXT 


eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee 

Model 

Mean 

Variance 

Skewness 

Kurtosis 

aMaaaaaaaaaaaa 

aaaaaaaaaaaaaaa 

aaaaaaaaaaaaaaa 

aaaaaaaaaaaaaaa 

aaaaaaaaaaaaaaa 

S  ,ple  values 

3.53333 

16 . 0410 

1.43210 

4.37573 

1 -Gamma  (E) 

3.53333 

21.3133 

2.73543 

14.2238 

2-Weibull 

3.52972 

18 . 3476 

2.66245 

14.2769 

3-Weibull  (E) 

3.67161 

30.2599 

3.86722 

28.6174 

4 -Gamma 

_ 3 . 53J3.3.1 _ 

16.4247 

2.29400 

10.8937 

5 -Rand.  Walk(E) 

3.53333 

21.1666 

2.77606 

14.5954 

6 -Random  Walk 

3.53333 

18.5726 

2.61623 

13.3952 

7 -Lognormal 

4.43168 

129.208 

24.5693 

4339.30 

8-Lognormal  (E) 

7.99936 

3425.23 

438.151 

1 . 07136E+07 

9-Inv.  Gaussian 

3.53333 

56.1879 

6 . 36441 

70.5095 

A-Exponential 

3.53333 

12.4844 

2.00000 

9.00000 

B-Expo.  (E) 

3.53333 

11.3936 

2.00000 

9.00000 

C-In.  Gaus.  (E) 

3.53333 

225.057 

13.3333 

299.295 

eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee 


Press  Fl-4  for  help  or  ENTER  to  continue:  Logfile:  Not  Open 

Model  Moment  Comparison  Sample:  Data  From  A_ARRST1.TXT  ; 

eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee 


Model 

Mean 

Variance 

Skewness 

Kurtosis 

a  aaaaaaaaaaaa 

aaaaaaaaaaaaaaa 

aaaaaaaaaaaaaaa 

aaaaaaaaaaaaaaa 

aaaaaaaaaaaaaaa 

Sample  Values 

3.53333 

16 . 0410 

1.43210 

4.37573 

1-Gamma  (E) 

3.53333 

21.3133 

2.73543 

14.2238 

2-Weibull 

3.52972 

18 . 3476 

2.66245 

14.2769 

3-Weibull  (E) 

3.67161 

30.2599 

3.86722 

28.6174 

4 -Gamma 

3.53333 

16.4247 

2.29400 

10.8937 

5-Rand.  Walk(E) 

3.53333 

21 . 1666 

2.77606 

14.5954 

6 -Random  Walk 

3.53333 

18.5726 

2.61623 

13 . 3952 

7 -Lognormal 

4.43168 

129.208 

24.5693 

4339.30 

8 -Lognormal  (E) 

7.99936 

3425 . 23 

438.151 

1 . 07136E+07 

9-Inv.  Gaussian 

3.53333 

56.1879 

6 . 36441 

70.5095 

A-Exponential 

3.53333 

12.4844 

2.00000 

9.00000 

B-Expo.  (E) 

3.53333 

11.3936 

2.00000 

9.00000 

C-In.  Gaus.  (E) 

3.53333 

225.057 

13.3333 

299.295 

eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee 

Press  Fl-4  for  help  or  ENTER  to  continue:  Logfile:  Not  Open 


Summary  of  Sample:  Data  From  A_SVRST1.TXT 
eSeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee 


Sample  Characteristic 

f  '  aaaaaaaaaaaaaaaaaaaaaa 

l  .servation  Type 

Number  of  Observations 

Minimum  Observation 

Maximum  Observation 

Mean 

Median 

Variance 

Skewness 


Value 

aaaaaaaaaaaaaaa 
Real  Valued 
67 

1.00000 

26.0000 

14.1493 

14.0000 

21.9774 

.07682 


Guided  Selection  Model  Rankings  For  Sample:  Data  From  A_SVRST1.TXT 
eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee 

P  ge  of  Random  Variable 

During  the  fitting  process  UniFit  considers  distributions  having  any 
reasonable  range  (not  just  the  specified  range),  provided  they  produce 
values  in  the  specified  range  at  least  99.99%  of  the  time. 

Specified  random  variable  range  At  least  0. 

Relative  Evaluation  of  Candidate  Models 

Relative 
Score 

Models  (0-100) 

aaaaaaaaaaaaaaaaaaaaaa  aaaaaaaa 

1- Weibull  95.0 

2- Weibull  (E)  92.5 

3- Extreme  Value  Type  B  85.0 

4- Gamma  82.5 

5- Log-logistic  81.3 

In  addition,  16  other  models  wer 
Current  Primary  Model 
1-Weibull 

i  jolute  Evaluation  of  the  Primary  Model 

Based  on  a  heuristic  evaluation,  there  is  no  current  evidence  for  not  using 
the  primary  model.  If  you  are  doing  simulation,  then  the  primary  model  will 
probably  provide  a  good  representation  for  your  data.  However,  we  recommend 
further  confirmation  of  the  primary  model.  Press  F3  for  more  information. 

Additional  Information  About  the  Primary  Model 

Result  of  an  Anderson-Darling 

goodness-of-fit  test  at  level  0.1  Do  not  reject 


=  .  39% 


Random  Variable  Range 

(if  different  from  that  specified) 

aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 

At  least  .05882 
Unrestricted 


e  considered  having  scores  from 


.0  to  81 


"Error"  in  the  model  mean 
relative  to  the  sample  mean 


.  05582 


M'odel  Moment  Comparisons  With  Sample:  Data  From  A  SVRSTl.TXT 
eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee 


M  del 

c.  aaaaaaaaaaaaa 
Sample  Values 

1-Weibull _ 

^2-Weibull  (E ) 

3- Ext.  Value  B 

4 - Gamma 

5 - Log-logistic 

6- Log-logis. (E) 

7 - Gamma  (E) 

8- Log-Laplace 

9- Log-Lap.  (E) 

A-Lognormal 
B-Lognormal  (E) 
C-Random  Walk 
D-Rand.  Walk ( E ) 
E-Inv.  Gaussian 
F-Pearson  6 
G-In .  Gaus .  (E) 

H-Pearson  6  (E) 
I-Pearson  5 
J-Pearson  5  (E) 
K-Expo.  (E) 
L-Exponential 


Mean 

aaaaaaaaaaaaaaa 

14.1493 

14.0934 
14.0909 
14.5305 

14.1493 
14.6845 
14.6925 

14.1493 
15.3953 
15.4101 
14 . 5498 
14 .5696 

14 . 1493 

14.1493 

14.1493 
15 .5893 

14.1493 
15 .6816 
16 .5942 
16 .8266 

14 . 1493 

14 .1493 


Variance 
aaaaaaaaaaaaaaa 
21 .9774 

22.4297 
22.4513 
35 . 5956 
28 . 9137 
15.4080 
15 .6614 
29 .1362 
70.4480 
71 .5865 
47.8327 
48.9761 
54 . 5156 
56 . 3534 
57 . 3606 
132.546 
59.5196 
142 .769 
285 . 262 
339 . 170 
198.540 
200 . 201 


Skewness 

aaaaaaaaaaaaaaa 

.07682 

. 08521 

709018 - 

1 .13955 
.76006 
2 . 76657 
2 . 80616 
. 76616 
13.5037 
14 . 4007 
1 .53343 
1 . 55903 
1 .43581 
1 . 46102 
1 . 60581 
5 . 87716 
1 . 64258 
6 . 59030 

Does  Not  Exist 
Does  Not  Exist 
2.00000 
2.00000 


Kurtosis 
aaaaaaaaaaaaaat 
. 2.874  34 

2.71145  } 

2.71189 
5.40000 
3 .86654 
40 .6371 
42.5799 
3 .88051 

Does  Not  Exist 
Does  Not  Exist 
7.45283 
7 .61150 
6 . 32929 
6 .44350 
7 . 29772 

Does  Not  Exist 
7.49679 

Does  Not  Exist 
Does  Not  Exist 
Does  Not  Exist 
9 .00000 
9 . 00000 


Models  Available  For  Sample:  Data  From  A_SVRST1.TXT 
eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee 

M'del  1:  Weibull  Distribution  X' 


Location  Parameter 

/  o.  ■fjr' 

Default 

Scale  Parameter 

/  15.717  7^X4-- 

M.L.  Estimate 

Shape  Parameter 

(  3.27388  lb\ 

M.L.  Estimate 

Model  2:  Weibull  Distribution  _ 

Location  Parameter 

.05882 

Quantile  Estimate 

Scale  Parameter 

15.6535 

M.L.  Estimate 

Shape  Parameter 

3.25628 

M.L.  Estimate 

Model  3 :  Extreme  Value 

Type  B  Distribution 

Location  Parameter 

11.8454 

M.L.  Estimate 

Scale  Parameter 

4.65183 

M.L.  Estimate 

Model  4:  Gamma  Distribution 

Location  Parameter 

0. 

Default 

Scale  Parameter 

2.04348 

M.L.  Estimate 

Shape  Parameter 

6.92411 

M.L.  Estimate 

Model  5:  Log-logistic 

Distribution 

Location  Parameter 

0. 

Default 

Scale  Parameter 

13.6332 

M.L.  Estimate 

Shape  Parameter 

4.74092 

M.L.  Estimate 

Model  6 :  Log-logistic 

Distribution 

Location  Parameter 

.05882 

Quantile  Estimate 

Scale  Parameter 

13.5722 

M.L.  Estimate 

Shape  Parameter 

4 .70930 

M.L.  Estimate 

Model  7 :  Gamma  Distribution 

Location  Parameter 

.05882 

Quantile  Estimate 

Scale  Parameter 

2.06780 

M.L.  Estimate 

Shape  Parameter 

6 . 81422 

M.L.  Estimate 

Model  8:  Log-Laplace  Distribution 

Location  Parameter 

0  . 

Default 

Scale  Parameter 

14.0000 

M.L.  Estimate 

Shape  Parameter 

3 . 32170 

M.L.  Estimate 

Model  9:  Log-Laplace  Distribution 

Location  Parameter 

.05882 

Quantile  Estimate 

Scale  Parameter 

13.9412 

M.L.  Estimate 

Shape  Parameter 

3 . 29945 

M.L.  Estimate 

Model  A:  Lognormal  Distribution 

Location  Parameter 

0. 

Default 

Scale  Parameter 

2.57572 

M.L.  Estimate 

Shape  Parameter 

.45135 

M.L.  Estimate 

Model  B:  Lognormal  Distribution 

Location  Parameter 

.05882 

Quantile  Estimate 

Scale  Parameter 

2.57033 

M.L.  Estimate 

Shape  Parameter 

.45730 

M.L.  Estimate 

Model  C:  Random  Walk  Distribution 

Location  Parameter 

0  . 

Default 

*  Scale  Parameter 

.09092 

M.L.  Estimate 

Shape  Parameter 

.  31735 

M.L.  Estimate 

Model  D:  Random  Walk  : 

Distribution 

Location  Parameter 

.  05882 

Quantile  Estimate 

Scale  Parameter 

.  09225 

M.L.  Estimate 

Shape  Parameter 

.  30771 

M.L.  Estimate 

Model  E:  Inverse  Gaussian  Distribution 

Location  Parameter 

0. 

Default 

Scale  Parameter 

14.1493 

M.L.  Estimate 

Shape  Parameter 

49 . 3841 

M.L.  Estimate 

Model  F:  Pearson  Type 

6 

Distribution 

Location  Parameter 

0. 

Default 

Scale  Parameter 

1.00000 

Default 

Shape  1  Parameter 

46 . 0063 

M.L.  Estimate 

Shape  2  Parameter 

3 .95115 

M.L.  Estimate 

Model  G:  Inverse  Gaussian  Distribution 

Location  Parameter 

.05882 

Quantile  Estimate 

Scale  Parameter 

14.0904 

M.L.  Estimate 

Shape  Parameter 

47.0016 

M.L.  Estimate 

Model  H:  Pearson  Type 

6 

Distribution 

Location  Parameter 

. 05882 

Quantile  Estimate 

Scale  Parameter 

1 . 00000 

Default 

Shape  1  Parameter 

44 . 0406 

M.L.  Estimate 

Shape  2  Parameter 

3 .81899 

M.L.  Estimate 

3el  I :  Pearson  Type 

5 

Distribution 

Location  Parameter 

0. 

Default 

Scale  Parameter 

32.6130 

M.L.  Estimate 

Shape  Parameter 

2.96532 

M.L.  Estimate 

Model  J :  Pearson  Type 

5 

Distribution 

Location  Parameter 

.05882 

Quantile  Estimate 

Scale  Parameter 

30 .6675 

M.L.  Estimate 

Shape  Parameter 

2 . 82896 

M.L.  Estimate 

Model  K:  Exponential  Distribution 

Location  Parameter 

.  05882 

Quantile  Estimate 

Scale  Parameter 

14 . 0904 

M.L.  Estimate 

Model  L:  Exponential  Distribution 

Location  Parameter 

0. 

Default 

Scale  Parameter 

14.1493 

M.L.  Estimate 

Summary  of  Sample:  Data  From  A_SVRINT.TXT 
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Sample  Characteristic 

a ' ' aaaaaaaaaaaaaaaaaaaaa 

Observation  Type 

Number  of  Observations 

Minimum  Observation 

Maximum  Observation 

Mean 

Median 

Variance 

Skewness 


Value 

aaaaaaaaaaaaaaa 
Real  Valued 
6  5 

2.00000 
12.0000 
4.47692 
4.00000 
3 . 87837 
1.52038 


GUided  Selection  Model  Rankings  For  Sample:  Data  From  A_SVRINT.TXT 
eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee 

P  nge  of  Random  Variable 

During  the  fitting  process  UniFit  considers  distributions  having  any 
reasonable  range  (not  just  the  specified  range),  provided  they  produce 
values  in  the  specified  range  at  least  99.99%  of  the  time. 

Specified  random  variable  range  At  least  0. 

Relative  Evaluation  of  Candidate  Models 


Relative 

Score 


Models  (0-100) 

aaaaaaaaaaaaaaaaaaaaaa  aaaaaaaa 

1- Pearson  Type  6  80.6 

2- Log-logistic  77.8 

3- Pearson  Type  5  69.4 

4- Lognormal  68.1 

5- Log-Laplace  68.1 


Random  Variable  Range 

(if  different  from  that  specified) 

aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 


In  addition,  14  other  models  were  considered  having  scores  from  .0  to  66.':' 


Current  Primary  Model 


1-Pearson  Type  6 

solute  Evaluation  of  the  Primary  Model 

Based  on  a  heuristic  evaluation,  we  recommend  being  cautious  about  using  the 
primary  model.  If  you  are  doing  simulation,  then  this  model  may  or  may  not 
provide  an  adequate  representation  for  your  data.  We  strongly  recommend 
further  confirmation  of  the  primary  model.  Press  F3  for  more  information. 

Additional  Information  About  the  Primary  Model 

Result  of  an  Anderson-Darling 

goodness-of-fit  test  at  level  0.1  Not  applicable 
"Error"  in  the  model  mean 

relative  to  the  sample  mean  .00313  =  .07% 


Model  Moment  Comparisons  With  Sample: ^Data^From^A^SVRINT. TXT 
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M  3el 


Mean 


Variance  _  Skewness  __ 

jaaaaaaaaaaaa  aaaaaaaaaaaaaaa  aaaaaaaaaaaaaaa  aaaaaaaaaaaaa 

4.47692  3.87837  1.52038 


Sample  Values 


Kurtosis 
aaaaaaaaaaaaaa  ! 
5 .74908 


1- Pearson  6 

2- Log-logistic 

3- Weibull  (E) 

4- Pearson  5 

5- Inv.  Gaussian 

^6  -  Lognormal _ 

7- Log-Laplace 

8 - Gamma  (E) 

9- Random  Walk 
A-Ext.  Value  B 
B-Johnson  SB 
C-Inv.  Weibull 
D-Expo.  (E) 
E-Gamma 

F-Rand 7  Walk ( E  j 
G-Lognormal  (E) 
H-In.  Gaus.  (E) 
I -Weibull 
J-Exponential 


4.47380 

4.43541 

4.47017 

4 .48849 

4.47692 

4^46870 

4  .’37  252 

4.47692 

4.47692 

4.43230 

4.45366 

4 .82686 

4.47692 

4.47692 

4.47692 

4 .95612 

4 . 47692 

4.48357 

4.47692 


3 . 80285 
1.67977 
3.83547 
4.21873 
3.40973 

3.44911 _ 

5 . 15222 
4 . 36009 
3 .33766 
3 .03947 

Can  Not  Compute 
16 . 1781 
6.77000 
3 . 18301 
7764405 
17 .6025 
14.0194 
3 .99428 
20 . 0428 


1 . 87556 
3 . 09177 
1 . 28600 
2.31523 
1 .23738 

1 . 31857 _ 

10.5353 
1.60503 
1 . 15928 
1 .13955 

Can  Not  Compute 
Does  Not  Exist 
2 . 00000 
.79702 
2 . 41509 
6.60992 
4 . 31710 
.41049 
2 . 00000 


10 . 8166 
60.8163 
5 .18220 
16 .0992 
5 .55183 
6.24283 

Does  Not  Exist 
6.86418 
5 . 19403 
5 .40000 

Can  Not  Compute 

Does  Not  Exist 
9 . 00000 
3 . 95286 
11  79611 
134 . 307 
34 . 0622 
2.91202 
9 . 00000 


Models  Available  For  Sample:  Data  From  A_SVRINT.TXT 
eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee 

M  lei  1:  Pearson  Type  6  Distribution 


Location  Parameter 

0, 

Default 

Scale  Parameter 

1.00000 

Default 

Shape  1  Parameter 

33.2831 

M.L.  Estimate 

Shape  2  Parameter 

8.43956 

M.L.  Estimate 

Model  2:  Log-logistic  Distribution 

Location  Parameter 

0. 

Default 

Scale  Parameter 

4.08519 

M.L.  Estimate 

Shape  Parameter 

4.50918 

M.L.  Estimate 

Model  3:  Weibull  Distribution 

Location  Parameter 

1.87500 

Quantile  Estimate 

Scale  Parameter 

2.82580 

M.L.  Estimate 

Shape  Parameter 

1 .33869 

M.L.  Estimate 

Model  4 :  Pearson  Type  5  Distribution 

Location  Parameter 

0. 

Default 

Scale  Parameter 

25.9233 

M.L.  Estimate 

Shape  Parameter 

6.77550 

M.L.  Estimate 

Model  5:  Inverse  Gaussian  Distribution 

Location  Parameter 

0. 

Default 

Scale  Parameter 

4.47692 

M.L.  Estimate 

Shape  Parameter 

26 . 3159 

M.L.  Estimate 

del  6 :  Lognormal  Distribution 

Location  Parameter 

0. 

Default 

Scale  Parameter 

1.41743 

M.L.  Estimate 

Shape  Parameter 

. 39916 

M.L.  Estimate 

Model  7 :  Log-Laplace  Distribution 

Location  Parameter 

0  . 

Default 

Scale  Parameter 

4.00000 

M.L.  Estimate 

Shape  Parameter 

3.42604 

M.L.  Estimate 

Model  8:  Gamma  Distribution 

Location  Parameter 

1.87500 

Quantile  Estimate 

Scale  Parameter 

1.67572 

M.L.  Estimate 

Shape  Parameter 

1.55272 

M.L.  Estimate 

Model  9:  Random  Walk  Distribution 

Location  Parameter 

0. 

Default 

Scale  Parameter 

.26137 

M.L.  Estimate 

Shape  Parameter 

1.53635 

M.L.  Estimate 

Model  A:  Extreme  Value  Type  B  Distribution 

Location  Parameter 

3.64767 

M.L.  Estimate 

Scale  Parameter 

1.35933 

M.L.  Estimate 

Model  B:  Johnson  SB  Distribution 

Lower  Endpoint  Parameter 

1.08579 

Quantile  Estimate 

Upper  Endpoint  Parameter 

47 .6248 

Quantile  Estimate 

Shape  1  Parameter 

3.81165 

Quantile  Estimate 

Shape  2  Parameter 

1 . 37684 

Quantile  Estimate 

Model-  'C:  Inverted  Weibull  Distribution 


Location  Parameter 

0. 

Default 

Scale  Parameter 

3.40047 

M.L.  Estimate 

Shape  Parameter 

2.71978 

M.L.  Estimate 

M  ^el  D:  Exponential  Distribution 
Location  Parameter 

L 

1 .87500 

Quantile  Estimate 

Scale  Parameter 

2 .60192 

M.L.  Estimate 

Model  E:  Gamma  Distribution 

Location  Parameter 

0  . 

Default 

Scale  Parameter 

.71098 

M.L.  Estimate 

Shape  Parameter 

6.29682 

M.L.  Estimate 

Model  F:  Random  Walk  Distribution 
Location  Parameter 

1.87500 

Quantile  Estimate 

Scale  Parameter 

1 .18021 

M.L.  Estimate 

Shape  Parameter 

.56993 

M.L.  Estimate 

Model  G:  Lognormal  Distribution 

Location  Parameter 

1 . 87500 

Quantile  Estimate 

Scale  Parameter 

.60090 

M.L.  Estimate 

Shape  Parameter 

1.02411 

M.L.  Estimate 

Model  H:  Inverse  Gaussian  Distribution 
Location  Parameter  1.87500 

Quantile  Estimate 

Scale  Parameter 

2 .60192 

M.L.  Estimate 

Shape  Parameter 

1 . 25647 

M.L.  Estimate 

Model  I:  Weibull  Distribution 

Location  Parameter 

0. 

Default 

Scale  Parameter 

5 . 05816 

M.L.  Estimate 

Shape  Parameter 

2 . 38848 

M.L.  Estimate 

Model  J :  Exponential  Distribution 

Location  Parameter 

0. 

Default 

Scale  Parameter 

4.47692 

M.L.  Estimate 

Summary  of  Sample:  Data  From  A_SVRBAG.TXT 
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S-mple  Characteristic 

a  aaaaaaaaaaaaaaaaaaaaa 

Observation  Type 

Number  of  Observations 

Minimum  Observation 

Maximum  Observation 

Mean 

Median 

Variance 

Skewness 


Value 

aaaaaaaaaaaaaaa 
Real  Valued 
54 

2.00000 
4.00000 
2.66667 
3 . 00000 
.33962 
. 18713 


Guided  Selection  Model  Rankings  For  Sample:  Data  From  A  SVRBAG.TXT 
eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee 

F'^ge  of  Random  Variable 

During  the  fitting  process  UniFit  considers  distributions  having  any 
reasonable  range  (not  just  the  specified  range),  provided  they  produce 
values  in  the  specified  range  at  least  99.99%  of  the  time. 

Specified  random  variable  range  At  least  0. 

Relative  Evaluation  of  Candidate  Models 


Models 

aaaaaaaaaaaaaaaaaaaaaa 

1 - Gamma 

2- Random  Walk 

3- Lognormal 

4- Normal 

5- Inverse  Gaussian 


Relative 

Score 

(0-100) 

aaaaaaaa 

69 . 2 

67 . 3 
63.5 

63.5 

61.5 


Random  Variable  Range 

(if  different  from  that  specified) 

aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 


3  3 


Unrestricted 


In  addition. 


9  other  models  were  considered  having  scores  from 


9.6  to  59.6 


Current  Primary  Model 


1 -Gamma 


7  olute  Evaluation  of  the  Primary  Model 

Based  on  a  heuristic  evaluation,  we  do  not  recommend  using  the  primary  mode 
If  you  are  doing  simulation,  then  you  should  use  an  empirical  distribution 
rather  than  the  primary  model  (unless  you  can  show  that  it  is  good). 

Press  F3  for  more  information. 

Additional  Information  About  the  Primary  Model 

Result  of  an  Anderson-Darling 
goodness-of-f it  test  at  level  0.1  Reject 

"Error"  in  the  model  mean 

relative  to  the  sample  mean  0. 


Model  Moment  Comparisons  With  Sample:  Data  From  A  SVRBAG.TXT 
eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee 


M^del 

Mean 

Variance 

Skewness 

Kurtosis 

L  _aaaaaaaaaaaa 

aaaaaaaaaaaaaaa 

aaaaaaaaaaaaaaa 

aaaaaaaaaaaaaaa 

aaaaaaaaaaaaaa 

Sample  Values 

2.66667 

.33962 

.18713 

2 . 24771 

1 -Gamma 

2 .66667 

.33746 

.43569 

3 . 2847 2<~" 

2-Random  Walk 

2.66667 

.  35039 

.  65505 

3.71052  ! 

3-Lognormal 

2 .66872 

.  36042 

.  68626 

3.84890  | 

4-Normal 

2.66667 

.33962 

0. 

3.00000^ 

5-Inv.  Gaussian 

2.66667 

.  35117 

.  66667 

3 .74074 

6-Weibull 

2 .66394 

.  36672 

-.25986 

2 .88660 

7-Pearson  6 

2 .66849 

. 36475 

.  81217 

4 . 27083 

8-Pearson  5 

2.66983 

. 38010 

.  97572 

4 . 87144 

9-Log-logistic 

2.71521 

.  14508 

1  .  36014 

9 .45007 

A-Ext.  Value  B 

2 .67177 

.41190 

1 .13955 

5.40000 

B-Inv.  Weibull 

2 .70853 

.70778 

3 . 48848 

46 . 0233 

C-Pareto  (E) 

2.72139 

1.15584 

79 . 2119 

Does  Not  Exist 

D-Exponential 

2 .66667 

7 . 11111 

2.00000 

9 .00000 

E-Log-Laplace 

3.09329 

.66622 

2. 35109 

24 .5857 

a  ■*  * 

Models  Available  For  Sample:  Data  From  A_SVR3AG.TXT 
eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee 


M^del  1:  Gamma  Distribution 

Location  Parameter 

0. 

Default 

Scale  Parameter 

.12655 

M.L.  Estimate 

Shape  Parameter 

21.0723 

M.L.  Estimate 

Model  2:  Random  Walk  Distribution 

Location  Parameter 

0. 

Default 

Scale  Parameter 

.  39352 

M.L.  Estimate 

Shape  Parameter 

7.96875 

M.L.  Estimate 

Model  3 :  Lognormal  Distribution 

Location  Parameter 

0. 

Default 

Scale  Parameter 

.95691 

M.L.  Estimate 

Shape  Parameter 

.22219 

M.L.  Estimate 

Model  4:  Normal  Distribution 

Location  Parameter 

2.66667 

M.L.  Estimate 

Scale  Parameter 

.58277 

M.L.  Estimate 

Model  5:  Inverse  Gaussian  Distribution 

Location  Parameter 

0. 

Default 

Scale  Parameter 

2.66667 

M.L.  Estimate 

Shape  Parameter 

54.0000 

M.L.  Estimate 

Model  6 :  Weibull  Distribution 

Location  Parameter 

0. 

Default 

Scale  Parameter 

2.89997 

M.L.  Estimate 

Shape  Parameter 

5 . 04168 

M.L.  Estimate 

Model  7 :  Pearson  Type  6  Distribution 

Location  Parameter 

0. 

Default 

Scale  Parameter 

1.00000 

Default 

Shape  1  Parameter 

74.2866 

M.L.  Estimate 

Shape  2  Parameter 

28 .8384 

M.L.  Estimate 

Model  8 :  Pearson  Type  5  Distribution 

Location  Parameter 

0. 

Default 

Scale  Parameter 

52.7365 

M.L.  Estimate 

Shape  Parameter 

20.7528 

M.L.  Estimate 

Model  9:  Log-logistic  Distribution 

Location  Parameter 

0  . 

Default 

Scale  Parameter 

2.63378 

M.L.  Estimate 

Shape  Parameter 

7 . 37275 

M.L.  Estimate 

Model  A:  Extreme  Value  Type  B  Distribution 

Location  Parameter 

.2.38293 

M.L.  Estimate 

Scale  Parameter 

.50040 

M.L.  Estimate 

Model  B:  Inverted  Weibull  Distribution 

Location  Parameter 

0. 

Default 

Scale  Parameter 

2 .33041 

M.L.  Estimate 

Shape  Parameter 

5 . 04440 

M.L.  Estimate 

Model  C:  Pareto  Distribution 

Location  Parameter 

1.99016 

Quantile  Estima 

*  S'ca'le  Parameter 


3 .72167 


M.L.  Estimate 


Model  D:  Exponential  Distribution 

Location  Parameter  0.  Default 

Scale  Parameter  2.66667  M.L.  Estimate 

Model  E:  Log-Laplace  Distribution 
Location  Parameter 
Scale  Parameter 
Shape  Parameter 


0. 

3.00000 

5.75827 


Default 
M.L.  Estimate 
M.L.  Estimate 


4  *  < 

Summary  of  Sample:  Data  From  A_SVRBLD.TXT 
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Sample  Characteristic 

t  .aaaaaaaaaaaaaaaaaaaaa 

Observation  Type 

Number  of  Observations 

Minimum  Observation 

Maximum  Observation 

Mean 

Median 

Variance 

Skewness 


Value 

aaaaaaaaaaaaaaa 
Real  Valued 
53 

8.00000 
39.0000 
21 .1698 
19.0000 
48.4898 
.67436 


Guided  Selection  Model  Rankings  For  Sample:  Data  From  A_SVRBLD.TXT 
eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee 


R^-ige  of  Random  Variable 


During  the  fitting  process  UniFit  considers  distributions  having  any 
reasonable  range  (not  just  the  specified  range),  provided  they  produce 
values  in  the  specified  range  at  least  99.99%  of  the  time. 

Specified  random  variable  range  At  least  0. 

Relative  Evaluation  of  Candidate  Models 


Models 

aaaaaaaaaaaaaaaaaaaaaa 

1- Pearson  Type  5 

2- Extreme  Value  Type  B 

3- Log-logistic  (E) 

4- Inverse  Gaussian 

5- Gamma  (E) 


Relative 
Score 
( 0-100 ) 
aaaaaaaa 

90.8 

85.5 

81.6 

78.9 
78.9 


Random  Variable  Range 

(if  different  from  that  specified) 

aaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaaa 


Unrestricted 
At  least  7.30435 


At  least  7.30435 


In  addition,  15  other  models  were  considered  having  scores  from  .0  to  72.4. 
Current  Primary  Model 


1-Pearson  Type  5 

A  olute  Evaluation  of  the  Primary  Model 

Based  on  a  heuristic  evaluation,  there  is  no  current  evidence  for  not  using 
the  primary  model.  if  you  are  doing  simulation,  then  the  primary  model  will 
probably  provide  a  good  representation  for  your  data.  However,  we  recommend 
further  confirmation  of  the  primary  model.  Press  F3  for  more  information. 

Additional  Information  About  the  Primary  Model 

Result  of  an  Anderson-Darling 

goodness-of-fit  test  at  level  0.1  Do  not  reject 
"Error"  in  the  model  mean 

relative  to  the  sample  mean  -.09635  =  .46% 


A  ^  A 

Model  Moment  Comparisons  With  Sample:  Data  From  A_SVRBLD.TXT 
eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee 


M-'^el 

Mean 

Variance 

Skewness 

Kurtosis 

a  .aaaaaaaaaaaa 

aaaaaaaaaaaaaaa 

aaaaaaaaaaaaaaa 

aaaaaaaaaaaaaaa 

aaaaaaaaaaaaaa 

Sample  Values 

21 . 1698 

48 . 4898 

.67436 

2.67380 

1-Pearson  5 

21 . 2662 

58 . 7141 

1 . 65629 

8 . 8  8  8  7  r 

2-Ext.  Value  B 

21 . 1766 

50.5320 

1 .13955 

5.40000 

3-Log-logis .(E) 

21 .9312 

93.3360 

18 . 5554 

Does  Not  Exist 

4-Inv.  Gaussian 

21 . 1698 

49.2338 

.99434 

4 . 64786 

5 -Gamma  (E) 

21 .1698 

54 . 0641 

1.06060 

4 .68730 

6-Lognormal 

21 . 1984 

50.4161 

1 . 04243 

4.99251 

7 -Random  Walk 

21 . 1698 

48 . 7515 

.  95440 

4 .49704 

8 -Log- logistic 

21 .1600 

21 . 4014 

2 .18730 

21 .4920 

9-Inv.  Weibull 

23 . 1135 

242 . 262 

300.015 

Does  Not  Exist 

A-Gamma 

21 .1698 

45 . 6677 

.63844 

3 .61140 

B-Weibull  (E) 

21 .1669 

47.7828 

.  5  6  253 . 

3'.  12454 

C-Lognormal  (E) 

21 . 7908 

97 .6036 

2.36312 

14 .3367 

D-Log-Lap.  (E) 

21 . 8184 

402 . 009 

Does  Not  Exist 

Does  Not  Exist 

E-Log-Laplace 

20 . 3890 

80 . 5110 

5 .99271 

Does  Not  Exist 

F-Weibull 

21 . 1752 

51 . 2203 

.  09111 

2 .71198 

G-Rand .  Walk(E) 

21 .1698 

96 . 0148 

1 . 82952 

8 . 30719 

H-Pearson  6  (E) 

24 . 0084 

701 . 517 

Does  Not  Exist 

Does  Not  Exist 

I-In.  Gaus.  (E) 

21 . 1698 

110.837 

2 . 27787 

11 . 6478 

J-Expo.  (E) 

21 . 1698 

192.251 

2.00000 

9 . 00000 

K -Exponential 

21 . 1698 

448.161 

2 . 00000 

9 . 0  0  0  0  C 

6  *  v  > 

Models  Available  For  Sample:  Data  From  A_SVRBLD.TXT 
eeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeeee 

M'-'del  1:  Pearson  Type  5  Distribution 


Location  Parameter 

Scale  Parameter 

Shape  Parameter 

0. 

185.070 

9.70256 

Default 

M.L.  Estimate 

M.L.  Estimate 

Model  2:  Extreme  Value  Type  B  Distribution 
Location  Parameter  17.9774 

Scale  Parameter  5.54254 

M.L.  Estimate 

M.L.  Estimate 

Model  3:  Log-logistic  Distribut 
Location  Parameter 

Scale  Parameter 

Shape  Parameter 

ion 

7 . 30435 
12.3696 
3.18633 

Quantile  Estimate 
M.L.  Estimate 

M.L.  Estimate 

Model  4:  Inverse  Gaussian  Distr 
Location  Parameter 

Scale  Parameter 

Shape  Parameter 

ibution 

0. 

21 . 1698 
192.702 

Default 

M.L.  Estimate 

M.L.  Estimate 

Model  5:  Gamma  Distribution 
Location  Parameter 

Scale  Parameter 

Shape  Parameter 

7 . 30435 
3.89919 

3 . 55599 

Quantile  Estimate 
M.L.  Estimate 

M.L.  Estimate 

Model  6 :  Lognormal  Distribution 
Location  Parameter 

Scale  Parameter 

Shape  Parameter 

0. 

3.00076 
. 32609 

Default 

M.L.  Estimate 

M.L.  Estimate 

Model  7 :  Random  Walk  Distribution 

Location  Parameter  0. 

Scale  Parameter  .05243 

Shape  Parameter  .47722 

Default 

M.L.  Estimate 

M.L.  Estimate 

Model  8:  Log-logistic  Distribut 
Location  Parameter 

Scale  Parameter 

Shape  Parameter 

ion 

0. 

19.9696 

5 . 36120 

Default 

M.L.  Estimate 

M.L.  Estimate 

Model  9:  Inverted  Weibull  Distribution 
Location  Parameter  0. 

Scale  Parameter  17.1013 

Shape  Parameter  3.01295 

Default 

M.L.  Estimate 

M.L.  Estimate 

Model  A:  Gamma  Distribution 
Location  Parameter 

Scale  Parameter 

Shape  Parameter 

0  . 

2 . 15721 
9.81351 

Default 

M.L.  Estimate 

M.L.  Estimate 

Model  B:  Weibull  Distribution 
Location  Parameter 

Scale  Parameter 

Shape  Parameter 

7 . 30435 

15 . 6520 

2 . 10771 

Quantile  Estimate 
M.L.  Estimate 

M.L.  Estimate 

Model  C:  Lognormal  Distribution 
Location  Parameter 

7  .  30435 

Quantile  Estimate 

Scale  Parameter 

2.48225 

M.L.  Estimate 

Shape  Parameter 

.61800 

M.L.  Estimate 

Model  D:  Log-Laplace 

Distribution 

Location  Parameter 

7 . 30435 

Quantile  Estimate 

Scale  Parameter 

11 . 6957 

M.L.  Estimate 

Shape  Parameter 

2 . 26931 

M.L.  Estimate 

Model  E:  Log-Laplace 

Distribution 

Location  Parameter 

0 . 

Default 

Scale  Parameter 

19.0000 

M.L.  Estimate 

Shape  Parameter 

3.83127 

M.L.  Estimate 

Model  F:  Weibull  Distribution 

Location  Parameter 

0  . 

Default 

Scale  Parameter 

23.6232 

M.L.  Estimate 

Shape  Parameter 

3 . 25299 

M.L.  Estimate 

Model  G:  Random  Walk 

Distribution 

Location  Parameter 

7 . 30435 

Quantile  Estimate 

Scale  Parameter 

.11370 

M.L.  Estimate 

Shape  Parameter 

.19722 

M.L.  Estimate 

Model  H:  Pearson  Type 

6  Distribution 

Location  Parameter 

7 . 30435 

Quantile  Estimate 

Scale  Parameter 

1 . 00000 

Default 

Shape  1  Parameter 

23.7458 

M.L.  Estimate 

Shape  2  Parameter 

2.42156 

M.L.  Estimate 

Model  I:  Inverse  Gaussian  Distribution 

Location  Parameter 

7 . 30435 

Quantile  Estimate 

Scale  Parameter 

13.8655 

M.L.  Estimate 

Shape  Parameter 

24 . 0501 

M.L.  Estimate 

Model  J:  Exponential 

Distribution 

Location  Parameter 

7 . 30435 

Quantile  Estimate 

Scale  Parameter 

13.8655 

M.L.  Estimate 

Model  K:  Exponential  : 

Distribution 

Location  Parameter 

0 . 

Default 

Scale  Parameter 

21 .1698 

M.L.  Estimate 

Data 


# 

Year 

Month 

Day 

Place 

Forcast 

Drawn 

Defered 

Total  Arrivals 

1 

1994 

1 

■Mil 

20 

14 

2 

16 

2 

1994 

1 

20 

13 

3 

16 

3 

1994 

1 

5 

NNMC 

20 

18 

2 

20 

4 

1994 

1 

6 

VA  Med. 

35 

35 

10 

45 

5 

1994 

1 

7 

uses 

40 

49 

5 

54 

6 

1994 

1 

11 

Bupers 

35 

39 

1 

40 

7 

1994 

1 

12 

NFEC 

40 

33 

5 

38 

8 

1994 

1 

13 

USUHS 

30 

19 

3 

22 

9 

1994 

1 

14 

USCG  HQ 

30 

32 

6 

38 

10 

1994 

1 

19 

NMRI 

20 

8 

0 

8 

11 

1994 

1 

21 

WNY 

35 

15 

1 

16 

12 

1994 

1 

24 

USNA 

25 

•  35 

3 

38 

13 

1994 

1 

25 

ONI 

40 

36 

6 

42 

14 

1994 

1 

26 

OSIA 

45 

47 

8 

55 

15 

1994 

1 

27 

NRC 

35 

32 

3 

35 

16 

1994 

2 

1 

USNA 

40 

55 

18 

73 

17 

2 

2 

ONI 

50 

48 

4 

52 

18 

2 

4 

NSS 

40 

41 

5 

46 

19 

1994 

2 

7 

USNA 

40 

57 

5 

62 

20 

1994 

2 

8 

Bupers 

35 

27 

5 

32 

21 

1994 

2 

10 

Dahlgrin 

45 

53 

2 

55 

22 

1994 

2 

15 

USNA 

40 

30 

6 

36 

23 

2 

16 

Camp  Dav 

50 

57 

3 

60 

24 

2 

17 

Pent 

30 

25 

4 

29 

25 

2 

18 

NRL  (civ) 

65 

63 

9 

72 

26 

2 

22 

NRL  (mil) 

35 

31 

3 

34 

27 

1994 

2 

23 

NIS 

35 

18 

11 

29 

28 

1994 

2 

24 

NSHS 

30 

20 

4 

24 

29 

1994 

2 

25 

Pax  Run 

65 

71 

9 

80 

30 

1994 

2 

28 

BUMED 

20 

17 

1 

18 

31 

1994 

3 

1 

USNA 

30 

26 

6 

32 

32 

1994 

3 

3 

Pent 

20 

29 

4 

33 

33 

1994 

3 

4 

PWBETH 

20 

11 

1 

12 

34 

1994 

3 

7 

Nav  Obs 

30 

20 

2 

22 

35 

1994 
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8 
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16 
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19 

36 

1994 
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40 
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45 

37 

1994 
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10 
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100 
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16 

141 

38 
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11 
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36 
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36 

39 

1994 
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14 

Dior 
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56 

5 

61 

40 

1994 

3 

15 

USNA 

35 

19 
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21 

41 

1994 
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17 

Pent 

30 

20 
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20 

42 

1994 
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18 
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35 

23 
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28 

1994 
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21 
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35 

27 

5 

32 
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1994 
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22 
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35 

10 

45 

45 

1994 

3 

24 

NRC 

30 

28 

6 

34 

46 

1994 

3 

29 

USNA 

35 

30 

0 

30 

47 

1994 

4 

1 

USCG 

25 

21 

0 

21 

48 

1994 

4 

4 

USNA 

25 

13 

1 

14 

49 

1994 

4 

5 

ONI 

35 

26 

5 

31 

50 

1994 

4 

7 

VAMed 

35 

50 

4 

54 
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Data 
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USCG 

Bupers 

ONI 

Dahlgrin 


Forcast 

Drawn 

Defered  Total  Arrivals 

30 

27 

2 

29 

30 

23 

5 

28 

35 

20 

5 

25 

35 

44 

3 

47 

A 


USNA 
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NRL 

USNA 


NNMC 

NIS 


Navy  Band 
Bupers 


NSHS 
Quantico 
Rescue  U 
Rescue  U 
Ft.  Meade 
Pent 
WNY 
NS  Attn 
Camp  Dav 
Quantico 
NRC 

NS  Ficility 
AF 
Pent 

NS  Station 
Quantico 
Dahlgrin 
NS  Station 
Bupers 
Pent 
Pax  Run 
Nav  Ohs 
NRL  (mil) 
NNMC 
Office 
NRL 
NCG 
NNMC 
AIMD 
Quantico 
Quantico 
Philly 
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Data 
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7 

7 

7 

7 

7 
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7 

7 

7 

7 

7 

7 

8 
8 
8 


Place  Forcast  Drawn  Defered  Total  Arrivals 
USCG  30  30  4  34 

SIA _ 40 _ j>4 _ 10 _ 64 

Bupers  35  50  7  57 

VA  Med  40  38  4  42 

WNY  30  26  3  29 

USUHS  30  30  1  31 

USCG _ 30 _ 25 _ 2 _ 27 

Pent  30  20  3  23 


OSIA 

NSA 


NNMC  30  20  7  27 

NRC  30  25  1  26 

NEOS  35  35  0  35 

NSHS _ 30 _ 20 _ 1_ _ 2\_ 

NASP  30  23  2  25 

NIS  30  25  7  32 


NS  Station 

Bupers 

Pent 

NRL 

NRL  (mil) 
NSGA 
Bupers 
Pax  Run 
NC 
NSNA 
PNSY 
NMRC 
USNA 
Navy  Band 
WNY 
USNA 
DIA 
Bupers 
Pent 
NNMC 
ONI 
NRC 
Quantico 
NSA 
USNA 
ONI 
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USNA 
NNMC 


21 

2 

23 

44 

7 

51 

19 

1 

20 

36 

3 

39 

23 

1 

24 

28 

4 

32 

20 

5 

25 

39 

6 

45 

110 

11 

121 

6 

0 

26 

32 

6 

38 

21 

3 

24 

20 

2 

22 

9 

2 

11 

28 
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31 

31 
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35 

21 
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21 

44 

6 
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33 
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13 
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16 

25 
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30 

17 

4 
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24 
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31 

44 

7 

51 

30 

11 

41 
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t 

Drawn 


_ 

1 

U 
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■ 

1 
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21 
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NIS 

NSA  Ann 
NNMC 
USNA 
USNA 
SS  Kenned 
NSA 
ONI 
NRC 
NSWC 
NMRI 
Bupers 
ONI 
Pent 
NRL 
BUMED 
NRL  (mil) 
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37 
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20 
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40 
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25 
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39 

40 

28 

3 

31 
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6 
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34 
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12 
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24 
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14 
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17 
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18 
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Data 


357 

1995 

358 

1995 

IHEE9I 

360 

■SEE 

Place 


USNA 

USCG 


VA  Med 
AIMD 


Bupers 

Dahlgrin 


NRL 

DIA 

USNA 

Pent 


Forcast 

Drawn 

Defered 

Total  Arrivals 

40 

28 

6 

34 

30 

18 

0 

18 

30 

22 

8 

30 

35 

42 

4 

46 

35 

26 

0 

26 

40 

16 

4 

20 

55 

49 

7 

56 

40 

22 

7 

29 

40 

40 

8 

48 

20 

8 

1 

9 

365 

1995 

10 

25 

USCG 

366 

1995 

10 

26 

NR  Comm 

367 

1995 

10 

27 

NNMC 

368 

1995 

10 

30 

NCG 

Page  8 


Data 


# 

Year 

Month 

Day 

Place 

Forcast 

Drawn 

Defered 

Total  Arrivals 

403 

1996 

1 

18 

Pax  Run 

40 

24 

0 

24 

404 

1996 

1 

19 

USUHS 

25 

26 

14 

40 

405 

1996 

1 

22 

Quantico 

40 

25 

2 

27 

406 

1996 

1 

23 

USNA 

40 

79 

16 

95 

407 

1996 

1 

25 

Pent 

20 

25 

7 

32 

408 

1996 

1 

26 

NSA 

30 

30 

13 

43 

409 

1996 

1 

29 

DIA 

30 

28 

7 

35 

410 

1996 

1 

30 

USNA 

40 

76 

6 

82 

411 

1996 

1 

31 

Marine  Bks 

40 

44 

15 

59 

412 

1996 

2 

5 

VAMed 

35 

13 

4 

17 

413 

1996 

2 

6 

USNA 

40 

86 

0 

86 

414 

1996 

2 

8 

Dahlgrin 

40 

24 

2 

26 

415 

1996 

2 

9 

Quantico 

45 

52 

4 

56 

416 

1996 

2 

12 

Bupers 

40 

43 

11 

54 

417 

1996 

2 

15 

Pax  Run 

50 

14 

4 

18 

418 

1996 

2 

16 

GW  NROTC 

40 

26 

9 

35 

419 

1996 

2 

20 

NEOD 

45 

25 

4 

29 

420 

1996 

2 

21 

NRL 

40 

21 

1 

22 

421 

1996 

2 

22 

Pent 

20 

17 

1 

18 

422 

1996 

2 

23 

ONI 

50 

46 

7 

53 

423 

1996 

2 

27 

USNA 

40 

55 

12 

67 

424 

1996 

2 

28 

NRC 

30 

27 

3 

30 

425 

1996 

2 

29 

NRL  (mil) 

35 

14 

4 

18 

# 

Year 

Month 

Day 

Forcast 

Drawn 

Defered 

Total  Arrivals 

15001 

13343 

1859 

14973 

AVE 

36.76716 

32.70343 

4.556373 

36.6985294 

STD 

15.41266 

21.60519 

3.724117 

24.138112 
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o\  WTWM&uc  iai“  "P\  in  <W  t 

vlorXaieet  size:  3500  cells  fUWK.  ^ 

MTB  >  Retrieve  'C:\JENNIFER\DIFFER.MTW'. 

Retrieving  worksheet  from  file: .C:\JENNIFER\DIFFER.MTW 
Worksheet  was  saved  on  5/12/1996 
-  >  nsco  c3  c4 
MTB  >  erase  c4 
MTB  >  let  c4  =  C2-C3 
MTB  >  nsco  c4  cl4 
MTB  >  plot  cl4  c4 

*  ERROR  *  Graph  type  is  not  allowed  for  this  command. 


MTB  >  STest  0.0  '94-95'; 
3UBC>  Alternative  0. 


Vi  fmwce  t^T5 


SIGN  TEST  OF  MEDIAN  =  0.00000  VERSUS  N.E.  0.00000 


N  BELOW  EQUAL  ABOVE  P-VALUE  MEDIAN 
94-95  12  5  0  7  0.7744  31.50 
MTB  >  WTest  0.0  '94-95'; 

3UBC>  Alternative  0. 


TEST  OF  MEDIAN  =  0.000000  VERSUS  MEDIAN  N.E.  0.000000 

N  FOR  WILCOXON  ESTIMATED 

N  TEST  STATISTIC  P-VALUE  MEDIAN 

94-95  12  12  48.0  0.505  19.25 

MTB  >  TTest  0.0  *94-95'; 

3UBC>  Alternative  0. 


J.—3T  OF  MU  =  0.0  VS  MU  N.E.  0.0 


N 

MEAN 

STDEV  SE 

MEAN 

T 

P  VALUE 

94-95 

12 

19.1 

149.5 

43.1 

0.44 

0.67 

MTB  > 

dotplot  cl 

c2 ; 

3UBC> 

same. 

• 

• 

•  •  ... 

•  • 

•  •  •  « 

_+ 

+ 

- .+ - 

- +— 

— Sums 

«  . 

•  .  a  •  .  • 

*  • 

•  mm 

-+ — 

- + - 

— + — - 

—  1994 

240 

360 

480 

600 

720 

840 

4TB  >  Mann-Whitney  95.0  '1994'  '1995'; 

!UBC>  Alternative  0.  g  JggOS 

4ann-Whitney  Confidence  Interval  and  Test 


1994  N  =  12  Median  =  528.5 

i-995  N  =  12  Median  =  520.5 

Point  estimate  for  ETA1-ETA2  is  27.0 

35.4  Percent  C.I.  for  ETA1-ETA2  is  (-48.0,101.0) 

=  164.0 

Lest  of  ETA1  =  ETA2  vs.  ETA1  ~=  ETA2  is  significant  at  0.4357 
Pannot  reject  at  alpha  =0.05 


MTB  >  Mann-Whitney  95.0  '1994'  '1995'; 

SUBC>  Alternative  1. 

M=>nn-Whitney  Confidence  Interval  and  Test 

1994  N  =  12  Median  =  528.5 

1995  N  =  12  Median  =  520.5 

Point  estimate  for  ETA1-ETA2  is  27.0 

95.4  Percent  C.I.  for  ETA1-ETA2  is  (-48.0,101.0) 

W  =  164.0 

Test  of  ETA1  =  ETA2  vs.  ETA1  >  ETA2  is  significant  at  0.2179 

Cannot  reject  at  alpha  =  0.05 

MTB  >  TwoSample  95.0  '1994'  '1995'; 

SUBC>  Alternative  0. 

TWOSAMPLE  T  FOR  1994  VS  1995 

N  MEAN  ST DEV  SE  MEAN 

1994  12  528.3  58.2  17 

1995  12  509  148  43 

95  PCT  Cl  FOR  MU  1994  -  MU  1995:  (  -80,  118) 

TTEST  MU  1994  =  MU  1995  (VS  NE) :  T=  0.41  P=0.68  DF=  14 

.4TB  > 


MTB  *> 

read  'jen94.dat' 

clOl 

cl02 

Entering  data 

from  file 

:  jen94.dat 

12  rows 

read . 

MTB  > 

read  'jen95.dat' 

cl03 

cl04 

Entering  data 

from  file 

:  jen95.dat 

12  rows 

read. 

MTB  > 

print  cl01-cl04 

ROW 

C101 

C102  C103 

C104 

1 

1 

425 

13 

439 

2 

2 

613 

14 

514 

3 

3 

541 

15 

866 

4 

4 

499 

16 

421 

5 

5 

458 

17 

442 

6 

6 

600 

18 

529 

7 

7 

581 

19 

391 v 

8 

8 

475 

20 

547 

9 

9 

543 

21 

598 

10 

10 

574 

22 

527 

11 

11 

515 

23 

588 

12 

12 

516 

24 

249 

MTB  > 

let  cl 

=  cl02 

MTB  > 

let  c2 

=  cl04 

MTB  > 

let  c3 

=  c2  -  cl 

MTB  >  name  cl  '1994'  c2  '1995'  c3  '95  -  94' 
MTB  >  dotplot  c3 


_ h _ + _ + - + - 1- - +-95  -  94 

-240  -120  0  120  240  360 

MTB  >  name  c4  'n  score' 

MTB  >  nsco  c3  c4 

MTB  >  plot  c4  c3 

n  score  - 

* 

1.2+  * 

* 

* 

* 

0.0+  *  * 

_  * 

_  * 

* 


-1.2  + 


* 


_ i _ i - 1 - 1 - 1 - 1-95  —  94 

-240  -120  0  120  240  360 

MTB  >  let  k90  =  3 

MTB  >  execute  'symplote' 

Executing  from  file:  symplote. MTB 


MTB  *>  read  'jen94.dat'  clOl  cl02 
En£ering  data  from  file:  jen94.dat 
12  rows  read. 

MTB  >  read  'jen95.dat'  cl03  cl04 
Entering  data  from  file:  jen95.dat 
12  rows  read. 

MTB  >  print  cl01-cl04 


ROW 

ClOl 

C102 

C103 

C104 

1 

1 

425 

13 

439 

2 

2 

613 

14 

514 

3 

3 

541 

15 

866 

4 

4 

499 

16 

421 

5 

5 

458 

17 

442 

6 

6 

600 

18 

529 

7 

7 

581 

19 

391 

8 

8 

475 

20 

547 

9 

9 

543 

21 

598 

10 

10 

574 

22 

527 

11 

11 

515 

23 

588 

12 

12 

516 

24 

249 

MTB  >  let  cl  =  cl02 
MTB  >  let  c2  =  cl04 
MTB  >  let  c3  =  c2  -  cl 

MTB  >  name  cl  '1994'  c2  '1995'  c3  '95  -  94' 
MTB  >  dotplot  c3 


»  •  •••••  •• 

- + - + - + - + - - - + 

-240  -120  0  120  240 


MTB  >  name  c4  'n  score' 

MTB  >  nsco  c3  c4 

MTB  >  plot  c4  c3 

n  score  - 


1.2  + 


0.0  + 


-1.2  + 


* 


- + - ■ — + - + - --+ - + 

-240  -120  0  120  240 

MTB  >  let  k90  =  3 

MTB  >  execute  'symplote' 

Executing  from  file:  symplotfe .MTB 


-+-95  - 
360 


— +95 
360 


C94 


-  * 
140  + 


70  + 


* 


* 


* 


* 


+ - + - + - - + - + - +- 

0  70  140  210  280  350 


MTB  >  end 
MTB  > 

MTB  >  execute  'skku' 

Executing  from  file:  skku. MTB 
MTB  >  print  k95  k96 

skewness  0.705876 
kurtosis  1.92659 
MTB  >  end 
MTB  >  ttest  0.0  c3 

TEST  OF  MU  =  0.0VSMUN.E.  0.0 

N  MEAN  STDEV  SE  MEAN  T  P  VALUE 

95  -  94  12  -19.1  149.5  43.1  -0.44  0.67 

MTB  >  stest  0.0  c3 

SIGN  TEST  OF  MEDIAN  =  0.00000  VERSUS  N.  E.  ~~0 . 00000 

N  BELOW  EQUAL  ABOVE  P-VALUE  MEDIAN 
95-94  12  7  0  5  0.7744  -31.50 

MTB  >  wtest  0.0  c3 

TEST  OF  MEDIAN  =  0.000000  VERSUS  MEDIAN  N.E.  OTOOOOOO 

N  FOR  WILCOXON  ESTIMATED 

N  TEST  STATISTIC  P-VALUE  MEDIAN 

95  -  94  12  12  30.0  0.505  -19.25 

MTB  >  let  k91  =  0.0 
MTB  >  execute  'johnson' 

Executing  from  file:  johnson. MTB 
MTB  >  print  kl05 

K105  -0.395013 

MTB  >  cdf  kl05  kl; 

SUBC>  t  11. 

MTB  >  let  kl  =  2*kl 
MTB  >  name  kl  'p-value' 

MTB  >  print  kl 

p-value  0.700390 
MTB  >  dotplot  cl  c2 ; 

SUBC>  same. 


C95 


• 

- +  -• 

- + - 

- ,+ - 1994 

- +-L - 1995 

240 

360 

480 

600 

720 

840 

MTB  >  nsco  cl  c4 
MTB  >  plot  c4  cl 

n  score  - 

* 


1.2  + 


0.0  + 


-1.2  + 


* 


+ 

420 


* 


* 


* 


* 


--+- 

455 


— +- 
490 


— +. 
525 


— +- 
560 


— +- 
595 


-1994 


MTB  >  let  k90  =  1 

MTB  >  execute  'symplote' 

Executing  from  file:  symplote. MTB 

C94 

* 


90  + 


* 


60  + 

* 


30  + 


_  *  * 

- + - + - + - + - + - +C95 

15  30  45  60  75  90 

MTB  >  end 
MTB  > 

MTB  >  execute  'skku' 

Executing  from  file:  skku. MTB 
MTB  >  print  k95  k96 

skewness  -0.227973 
kurtosis  -0.764211 
MTB  >  end 

,MTB  >  nsco  c2  c4  • 

MTB  >  plot  c4  c2 


n  score 


1.2  + 


0.0  + 


-1.2  + 


★ 

★ 

★ 

★  ★ 

* 

★ 


* 


★ 


+ - + - - — + - + - + - + 

240  360  480  .600  720  840 


MTB  >  let  k90  =  2 

MTB  >  execute  'symplote' 

Executing  from  file:  symplote. MTB 


300  + 

-  * 

C94 


200  + 


_  * 

100+  * 

_  *  * 


0+  * 


H - 1 - 1 - H - 1 - 1. 

o  70  140  210  280  350 

MTB  >  end 
MTB  > 

MTB  >  execute  ' skku ' 

Executing  from  file:  skku. MTB 
MTB  >  print  k95.  k96 


skewness  0.868077 
kurtosis  2.83897 
MTB  >  end 
MTB  >  twos  cl  c2 

TWO SAMPLE  T  FOR  1994  VS  1995 

N  MEAN  STDEV  SE  MEAN 

1994  12  528.3  58.2  16.8 

1995  12  509  148  42.9 


95  PCT  Cl  FOR  MU  1994  -  MU  1995:  (  -79.66,  117.8) 

TTEST  MU  1994  =  MU  1995  (VS  NE) :  T=  0.41  P=0.68  DF=  14 


1995 


C95 


MTB  >  twos  cl  c2; 

SUBC>  same. 

*  ERROR  *  Subcommand  not  found  in  dictionary. 

*  Subcommand  ignored. 

SUBC>  pooled. 

TWOSAMPLE  T  FOR  1994  VS  1995 

N  MEAN  STDEV  SE  MEAN 

1994  12  528.3  58.2  16.8 

1995  12  509  148  42.9 

95  PCT  Cl  FOR  MU  1994  -  MU  1995:  (  -76.40,  114.6) 

TTEST  MU  1994  =  MU  1995  (VS  NE) :  T=  0.41  P=0.68  DF=  22 

POOLED  STDEV  =  113  - 

MTB  >  mann  cl  c2 

Mann-Whitney  Confidence  Interval  and  Test 

1994  N  =  12  Median  =  528.5 

1995  N  =  12  Median  =  520.5 

Point  estimate  for  ETA1-ETA2  is  27.0 

95.4  pet  c.i.  for  ETA1-ETA2  is  (-48.0,101.0) 

W  =  164.0 

Test  of  ETAl  =  ETA2  vs.  ETA1  n.e.  ETA2  is  significant  at  0.4357 
Cannot  reject  at  alpha  =  0.05 
MTB  >  describe  cl  c2 


N 

MEAN 

MEDIAN 

TRMEAN 

STDEV 

SEME AN 

1994 

12 

528.3 

528.5 

530.2 

58.2 

16.8 

1995 

12 

509.3^ 

520.5 

499.6 

148.5 

42.9 

MIN 

MAX 

Q1 

~qT 

. . 

1994 

425.0 

613.0 

481.0 

579.2 

1995 

249.0 

866.0 

425.5 

577.7 

MTB  >  sort  cl  ell 
MTB  >  sort  c2  cl2 

MTB  >  name  ell  'sort  94'  cl2  'sort  95' 
MTB  >  print  ell  cl2 


ROW 

sort  94 

sort  95 

V 

1 

425 

249 

4  *  ■ 

2 

458 

391 

?  \  .  - 

3 

475 

421 

4 

499 

439 

5 

515 

442 

““ 

6 

516 

514 

- 

7 

541 

527 

8 

543 

529 

9 

574 

547 

10 

581 

588 

11  ' 

600 

598 

12 

613 

866 

MTB  > 

rank  cl 

c21 

MTB  > 

rank  c2 

c22 . 

MTB  > 

name  c21 

'rank  94'  c22  'rank  95' 

MTB->  print  c21  c22 


ROW  rank  94  rank  95 


1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 


1 

12 

7 

4 
2 

11 

10 

3 

8 
9 

5 

6 


4 
6 

12 

3 

5 
8 
2 
9 

11 

7 

10 

1 


MTB  >  save  'jenny' 

Saving  worksheet  in  file:  jenny. MTW 
MTB  >  stop 

***  Mini tab  Release  9.1  ***  Mini tab  Inc. 
Worksheet  size:  5310509  cells 


*** 


V 


c 


JjU7  <■'&  Ms*-  vJ/  /$'  ?'■ 


lU-  .  /'/~  v- 


The  observed  no.  of  runs  =  27 

The  expected rpo.  of \ runs  26.8^62 

24  \0hservati6ns  above  K  ;  .28  beiow\ 

v^\  V/Thi  'test/  isVsigiriifhkcant  aV  0^  9654 
x-'  >veannot>re  j  ecVjn:  al^hd  =  0.0^ 


jorksheet  size:  3500  cells 

TTB  >  Name  c2  =  'ACF1' 

ITB  >  ACF  24  'Sum (All)  '  'ACF1'. 

xv-F  of  Sum  (All) 


-1.0  -0.8  -0.6  -0.4 

+ - + - + - + 

1  -0.066 

2  -0.124 

3  -0.189 

4  0.198 

5  -0.166 

6  -0.068 

7  -0.066 

8  0.035 

9  0.036 
10  -0.008 

11  0.096 

12  0.075 

13  -0.027 

14  0.022 

15  -0.061 

16  0.150 

17  0.075 

18  -0.050 

19  -0.264 

20  -0.009 

1  -0.026 
1  -0.042 

23  0.097 

24  -0.025 


-0.2  0.0  0.2  0.4  0.6  0.8 

— + - + - + - + - + - +_ 

XXX 

xxxx 

xxxxxx 

xxxxxx 

xxxxx 

XXX 

XXX 

XX 

XX 

X 

XXX 

XXX 

XX 

XX 

XXX 

xxxxx 

XXX 

XX 

xxxxxxxx 

X 

XX 

XX 

XXX 

XX 


1.0 

-+ 


TB  > 


Final  Estimates  of  Parameters 


Type  Estimate  st.  Dev.  t-ratio 
AR  1  0.8449  0.1026  8.23 
MA  1  1.0052  0.0572  17.58 
'  istant  40.0585  0.3521  113.78 


i-dn  258.225  2.270 


No.  of  obs.:  52 
Residuals:  SS 

MS 


319593  (backforecasts  excluded) 
6522  DF  =  49 


Modified  Box-Pierce  (Ljung-Box)  chisquare  statistic 

Lag  12  24  36  48 

Chisquare  6.8(DF=10)  17.5(DF=22)  25.7(DF=34)  52.0(DF=46) 


Forecasts  from  period  52 


Period 

Forecast 

53 

283.250 

54 

279.368 

55 

276.088 

56 

273.317 

57 

270.976 

95  Percent  Limits 
Lower  Upper 

124.927  441.573 

119.023  439.712 

114.316  437.860 

110.533  436.100 

107.474  434.477 


Actual 


I f  u  it*'  %  ^  4b  0  ( r  ^ 

('Aftl  "  N/A - 


1 

o 

J 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 
21 
22 

23 

24 


-1.0  -0.8  -0.6 

+ - + - +- 

-0.066 

-0.124 

-0.189 

0.198 

-0.166 

-0.068 

-0.066 

0.035 

0.036 

-0.008 

0.096 

0.075 

-0.027 

0.022 

-0.061 

0.150 

0.075 

-0.050 

-0.264 

-0.009 

-0.026 

-0.042 

0.097 

-0.025 


-0.4  -0.2  0.0  0.2 

— + - + - + - +- 

XXX 

xxxx 

xxxxxx 

xxxxxx 

xxxxx 

XXX 

XXX 

XX 

XX 

X 

XXX 

XXX 

XX 

XX 

XXX 

xxxxx 

XXX 

XX 

xxxxxxxx 

X 

XX 

XX 

XXX 

XX 


0.4  0.6  0.8  1.0 

-+ — - — + - + - + 


P&CF  of  Bi-Mont 


?fl£T IAL  folJDCClie&lhTi ON 


-1.0  -0.8  -0.6  -0.4 
+ - + - + - +- 


-0.066 

z 

-0.129 

3 

-0.211 

4 

0.157 

5 

-0.208 

6 

-0.091 

7 

-0.064 

8 

-0.114 

9 

0.049 

10 

-0.060 

11 

0.094 

12 

0.091 

13 

-0.044 

14 

0.129 

15 

-0.084 

16 

0.213 

17 

0.193 

18 

-0.039 

19 

-0.097 

20 

-0.114 

21 

-0.115 

22 

-0.116 

23 

0.110 

24 

-0.135 

-0.2  0.0  0.2  0.4  0.6  0.8  1.0 

XXX 

xxxx 

xxxxxx 

xxxxx 

xxxxxx 

XXX 

XXX 

xxxx 

XX 

XXX 

XXX 

XXX 

XX 

xxxx. 

XXX 

xxxxxx 

xxxxxx 

XX 

XXX 

xxxx 

xxxx 

xxxx 

xxxx 

xxxx 


MTB  >  retrieve  'jen.mtw' 

•  *  ERROR  *  File  not  found:  jen.mtw 

MTB  >  retrieve  'jenreg.mtw' 

Retrieving  worksheet  from  file:  jenreg.mtw 
Worksheet  was  saved  on  5/12/1996 
MTB  >  rreg  c2  1  cl 


The  regression 

equation  is 

^  _ 
_  -»  ,  v*  .  1^-  • 

/ 

sum  =  531  -  1. 

19  month 

Coefficient 

St.  dev. 

of  coef. 

Predictor 

Rank  Least-sq 

Rank 

Least-sq 

Constant 

530.66  545.32 

39.88 

47.25 

month 

-1.190  -2.122  - 

2.791 

3.307 

•> » *■*  r)  Q  *'“> , J  c4  <■>  .i  **  •  • 

*w  Xs  0  r  f  j  ~ 

Hodges -Lehmann 

estimate  of'tau  =  94.65 

Least 

-squares  S  = 

^  /  A 


r  *  av*  *  * 

,  . 

„  .-o-r  ,  r '  / 

j  r 

{•i /-i  j  ^  v* 

* '  **  i 

vo  J  f'  * 

’  t. .  ~ 


r  / 


••  K  ^  /*)  €*-i 


112.1 


MTB  >  describe  c2 


sum 

N 

24 

MEAN 

518.8 

MEDIAN 

521.5 

TRMEAN 

515.3 

STDEV 

110.7 

SEME AN 
22.6 

sum 

MIN 

249.0 

MAX 

866.0 

Ql 

446.0 

Q3 

579.2 

MTB  >  nsco  c2  c6 
MTB  >  plot  c6  c2 


n  score  - 


1.5  + 


0.0  + 


* 

* 

*  * 

*  2 

*2 
*  * 

*2 

**  * 


★ 


-1.5  + 


★  * 
* 


* 


★ 


— + 
240 


--+ 

360 


- 1 - 1_ 

480  600 


MTB  >  corr  cl  c2 
Correlation  of  month  and  sum 


-0.13  6  -/■ 


o 

6 


MTB  >  stop  - '' 

***  Minitab  Release  9.1  ***  Minitab  Inc.  *** 
Worksheet  size:  5310509  cells 


5 

4 

— + - + - sum 

720  840 


uVlT'*  '■  A'. 


i  f- 


PiL  iMih  "  bUj o  j  AK-  ptg 

MTB  >  ARIMA  10  2  'Bi-Mont.1;  ~ - - ^ 

SUBC>  Constant ; 

SUBC>  Forecast  6  . 


r  '-.imates  at  each  iteration 


*  aration 

SSE 

Parameters 

0 

274153 

0.100 

0.100 

0.100 

233.546 

1 

270133 

0.228 

0.250 

0.121 

200.604 

2 

265353 

0.356 

0.400 

0.137 

167.568 

3 

258902 

0.480 

0.550 

0.148 

135.619 

4 

249802. 

0.589 

0.700 

0.159 

107.506 

5 

239225 

0.538 

0.756 

0.200 

121.469 

6 

235085 

0.442 

0.747 

0.199 

146.993 

7 

234581 

0.389 

0.698 

0.248 

160.853 

8 

234482 

0.345 

0.647 

0.297 

172.461 

9 

234469 

0.341 

0.643 

0.298 

173.643 

10 

234461 

0.344 

0.646 

0.296 

172.852 

Unable  to 

reduce  sum  of 

squares 

any  further 

Final  Estimates  of  Parameters 

Type  . 

Estimate 

St.  Dev. 

t-ratio 

AR  1 

0.3435A 

0.2983 

/"1.15'"n 

v 

MA  1 

0.6461 

0.303.3 

l  2.14 

) 

MA  2 

0.2956/ 

0.245a 

\_1. 21  J 

Constant 

1^7^7152^ 

0.805 

214.84 

Mean 

263.313 

1.226 

t" 


No.  of  obs. 
Residuals: 


48 

SS  = 
MS  = 


234306  (backforecasts  excluded) 
5325  DF  =  44 


Modified  Box-Pierce  (Ljung-Box)  chisquare  statistic 

Lag  12  24  36  48 

Chisquare  _ 2 . 3 ( DF=  9)  11 . 0(DF=21) _ 14.5(DF=33)  /  *  (DF=  *) 

Ol  ■*=■  p-pqXux  4  ,0  2-5)  .OlS^p-'OoJ^CjDS^'  <0  f0  ^  p  i.CS) 

Forecasts  from  period  48  ' 


95  Percent  Limits 


Period  Forecast 

_  Lower 

Upper 

49 

340.688 

197.632 

483.745 

50 

333.900 

184.438 

483.362 

51 

287.563 

127.545 

447.580 

52 

i  271.644  1 

110.426 

432.862 

53 

266.175 

104.816 

427.534 

54 

1  264.296 

102.920 

425.672 

MTB  > 

i  Xt  -- 

3 

\  ^ 

Actual 


* 


>1 / 


v 


MTB  >  ARIMA  202  'Bi-Mont. ' ; 
3UBC>  Constant ; 

3UBC>  Forecast  6  . 

r  u.imates  at  each  iteration 


^  oration 

SSE 

Parameters 

0 

282466 

0.100 

0.100 

0.100 

0.100 

207.597 

1 

260533 

0.028 

-0.012 

0.173 

0.213 

257.339 

2 

253810 

0.066 

0.129 

0.230 

0.363 

210.799 

3 

244027 

0.124 

0.263 

0.320 

0.513 

160.738 

4 

236524 

0.102 

0.245 

0.369 

0.556 

171.506 

5 

235227 

0.062 

0.197 

0.375 

0.551 

195.055 

6 

234962 

0.105 

0.158 

0.423 

0.508 

194.203 

7 

234554 

0.255 

0.063 

0.565 

0.372 

179.643 

8 

234491 

0.300 

0.033 

0.606 

0.332 

175.700 

9 

234491 

0.296 

0.036 

0.602 

0.338 

175.904 

10 

234491 

0.296 

0.036 

0.602 

0.337 

176.087 

Onable  to 

reduce  sum  of  squares 

any  further 

Final  Estimates  of  Parameters 

Type 

Estimate 

St.  Dev. 

t-ratio  / 

AR  1 

0.2956 

1.0501 

0.28  ^ 

AR  2 

0.0356 

0.6691 

K0 . 05! 

MA  1 

0.6020 

1.0043 

0.60  ^ 

MA  2 

0.3375 

0.9862 

0.34  ^ 

Constant 

176.087 

0.865 

203.67 

Mean 

263.303 

1.293 

Mo.  of  obs 

i. :  48 

?osiduals: 

SS  =  234322  (backforecasts 

excluded) 

MS  = 

5449  DF  i 

=  43 

Modified  Box-Pierce  (Ljung-Box)  chisquare  statistic 


Lag 
Chisquare 


(2.3  (] 

Forecasts  from  period  48 


24 

10.9 (DF=20) 


36 

14.5 (DF=32) 


48 

(DF=  *) 


95  Percent  Limits 

Period 

Forecast 

Lower 

Upper 

49 

342.153 

197.437 

486.869 

50 

331.137 

179.780 

482.493 

51 

286.165 

124.505 

447.825 

52 

272.479 

109.779 

435.179 

53 

266.830 

103.959 

429.701 

54 

264.673 

101.777 

427.569 

Actual 


HUU  DATA  Mjo.  JAN-F I'd 


honiHL-i 


MTB  >  ARIMA  0  0  2  0  0  0  24  ' Bi-Mont.' ? 
SUBC>  Constant? 

SUBC>  Forecast  6  . 


timates  at  each  iteration 


.  aration 

i  SSE 

Parameters 

0 

269781 

0.100 

0.100 

259.496 

l 

258997 

0.166 

0.250 

261.581 

2 

253418 

0.235 

0.328 

262.040 

3 

249676 

0.285 

0.388 

262.488 

4 

247391 

0.323 

0.430 

262.834 

5 

246098 

0.348 

0.459 

263.112 

6 

245113 

0.367 

0.483 

263.333 

7 

243901 

0.385 

0.505 

263.544 

8 

242378 

0.402 

0.526 

263.719 

9 

242133 

0.411 

0.529 

263.656 

10 

242131 

0.412 

0.527 

263.633 

11 

242131 

0.411 

0.528 

263.632 

Relative 

change  in  each 

estimate 

less  than  0.00 

Final  Estimates  of  Parameters 

Type 

Estimate 

St.  Dev. 

t-ratio 

MA  1 

0.4lT5\ 

0.1323 

3.11 

MA  2 

0.5276/ 

0.1354 

3.90 

Constant 

Z&T7&32 

0.993 

265.47 

Mean 

263.632 

0.993 

A- 


No.  of  obs.:  48 

Residuals:  SS  =  241928  (backforecasts  excluded) 

MS  =  5376  DF  =  45 


Modified  Box-Pierce  (Ljung-Box)  chisquare  statistic 

Lag  12  24  36  48 

Chisquare  4.5(DF=10)  12.6(DF=22)  14.7(DF=34)  *  (DF=  *) 

Forecasts  from  period  48 

95  Percent  Limits 

Period  Forecast  Lower  Upper  Actual 

49  388.839  245.098  532.580 

50  334.487  179.061  489.912 

51  263.632  90.689  436.576 

52  263.632  90.689  436.576 


MTB  >  read  'jen9495.dat'  cl  c2 
Entering  data  from  file:  jen9495.dat 
24  rows  read1. 

MTB  >  name  cl  'month'  c2  'sum' 

MTB  >  plot  c2  cl 


sum 


800  + 


600  + 


400  + 


*  * 


★  * 


★  * 


★  ★ 


/  ^TL  11  ^ 

ifjj  z  <ruJj£<'cA^  fc- 


+■ 

0.0 


— +- 
5.0 


10.0 


15.0 


20.0 


25.0 


-month 


MTB  >  regress  c2  1  cl  c3  c4; 
SUBC>  tres  c5; 

SUBC>  dw. 

The  regression  equation  is 
sum  =  545  -  2.12  month 


Predictor 

Constant 

month 


Coef 

545.32 

-2.122 


St  dev 
47.25 
3.307 


t-ratio 

11.54 

-0.64 


P 

0.000 

0.528 


.. ..4. 


s  =  112.1  R-sq  =  1.8% 

Analysis  of  Variance 


R-sq(adj)  =  0.0% 


SOURCE 

DF 

SS 

MS 

F  p 

Regression  1 

5179 

5179 

0.41  0.528 

Error 

22 

276665 

12576 

Total 

23 

281844 

v  "  ■ 

Unusual 

Observations 

v  *■*  *  ' 

'  ~ 

s'  A 

Obs . 

month 

sum 

Fit  Stdev.Fit  Residual 

St .Resid 

^ 

15 

15.0 

866.0 

513.5 

24.3  352.5 

3 . 22R 

24 

24.0 

249.0 

494.4 

44.4  -245.4 

-2.38R 

A"'.  * 

R  denotes  an  obs .  with  a  large  st.  resid. 

Durbin-Watson  statistic  =  2.15  '=* 

MTB  >  name  c4  'pred  sum'  c5  'stud  res' 
MTB  >  plot  c5  cl 


stud  res- 


2.5+ 


★  * 


★  * 


0.0  + 


-2.5  + 


MTB  >  nsco  c5  c6 

MTB  >  name  c6  'n  score' 

MTB  >  plot  c6  c5 

n  score  - 


10.0 


15.0 


20.0 


-month 


25.0 


1.5  + 


ft  ffltOsl 


0.0  + 


1 1  js  **— f  ) 


f*  r*  O 


//_j  A  X  ' 


-1.5  + 


-1.5 


MTB  >  lag  c5  c7 

MTB  >  name  c7  'lag  res' 

MTB  >  plot  c7  c5 


—  4  .  0/ 


— + - stud  res 

4.5 


lag  res  ■- 


7  tvhaJ  jPa.  ku  if  yen  '3% 

*  {'f  3fu_^ 


'sr'fttCr&u 


V  d£ur^>t->^  4  NCrdacUf-a 


/KClD  ^ullcajlc^  //V 
f  3 


2.0  + 


o.o+ 


★  *  *  *  * 


★  *  *  * 


-3.0 
N*  =  1 


-stud  res 


MTB  >  plot  c5  c7 


stud  res- 

2.5+ 


of 


_  *  * 

-  **  *  *  *  * 
0.0+  ***  *  *  * 

*  *  2  * 


-2.5  + 


__+. 

-1.0 
N*  =  1 


- +- 

0.0 


- +- 

1.0 


- +- 

2.0 


— +- 
3.0 


— +- 
4.0 


■lag  res 


MTB  >  save  'jenreg' 

Saving  worksheet  in  file:  jenreg. MTW 
MTB  >  Is  jen* 

*  ERROR  *  Name  not  found  in  dictionary. 
MTB  >  stop 
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PROGRAM  INSTRUCTIONS 


/input 

/variables 

/test 

/end 


variables  are  2 . 
format  is  free, 
file  is  'jen9495.dat'. 
names  are  month,  sum. 
kendall . 
spearman. 


PROBLEM  TITLE  IS 

12 -MAY-9 6  17:28:57 


NUMBER  OF  VARIABLES  TO  READ .  2 

NUMBER  OF  VARIABLES  ADDED  BY  TRANSFORMATIONS.  .  0 

TOTAL  NUMBER  OF  VARIABLES  .  2 

CASE  LABELING  VARIABLES  . 

NUMBER  OF  CASES  TO  READ . TO  END 

MISSING  VALUES  CHECKED  BEFORE  OR  AFTER  TRANS.  .  NEITHER 

BLANKS  IN  THE  DATA  ARE  TREATED  AS  . MISSING 

INPUT  FILE.  .  .jen9495.dat 

REWIND  INPUT  UNIT  PRIOR  TO  READING.  .  DATA.  .  .  YES 

NUMBER  OF  INTEGER  WORDS  OF  MEMORY  FOR  STORAGE  .  19998 


VARIABLES  TO  BE  USED 

1  month  2  sum 


DATA  FORMAT :  FREE 


THE  LONGEST  RECORD  MAY  HAVE  UP  TO  80  CHARACTERS. 
USE  ONLY  COMPLETE  CASES 

COMPUTE  KENDALL  RANK  CORRELATION  COEFFICIENT (S) 
COMPUTE  SPEARMAN  RANK  CORRELATION  COEFFICIENT (S) 


NUMBER  OF  CASES  READ. 


VARIABLE 
NO.  NAME 

1  month 

2  sum 
1PAGE  2 


MEAN 


3S 


12.5000 
518.7917 
12 -MAY- 


96 


STANDARD 

DEVIATION 

7.0711 

110.6982 


MINIMUM 

1.0000 

249.0000 

17:28:57 


12.5000 

521.5000 


MAXIMUM  COUNT 

24.0000  24 

866.0000  24 


KENDALL  RANK  CORRELATION  COEFFICIENTS 


7 

* 


month 


1 


sum 


2 


month 

1 

1.0000 

*2  ~  ^  ^  "l  ■-* ;  V 

r  . 

sum 

2 

-0T  0797^) 

1.0000 

— — - •y  /  "T~.  <_  S.  izj  H 

\  XA  .Jr  } 

\  >>«*«/£  *!«■  ***/**>_  - 
I  *  >  .  w  U  «  \a, 

#  f  ; 

SPEARMAN 

RANK 

CORRELATION 

month 

1 

COEFFICIENTS 

sum 

2 

^  JL.  r*  ^  *  '  f  \  ^  ■ 

month 

sum 

1 

2 

1.0000 

;'^0T09'57J 

1.0000 

NUMBER  OF  INTEGER  WORDS  USED  IN  PRECEDING  PROBLEM  652 

CPU  TIME  USED  0.130  SECONDS 
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PROGRAM  INSTRUCTIONS 


/input  variables  are  2. 

format  is  free, 

file  is  'jen94.dat'. 
/variables  names  are  month,  sum. 

/test  kendall. 

spearman. 

/end 

PROBLEM  TITLE  IS 

12 -MAY- 9 6  17:32:23 


NUMBER  OF  VARIABLES  TO  READ .  2 

NUMBER  OF  VARIABLES  ADDED  BY  TRANSFORMATIONS.  .  0 

TOTAL  NUMBER  OF  VARIABLES  .  2 

CASE  LABELING  VARIABLES  . 

NUMBER  OF  CASES  TO  READ . TO  END 

MISSING  VALUES  CHECKED  BEFORE  OR  AFTER  TRANS.  .  NEITHER 

BLANKS  IN  THE  DATA  ARE  TREATED  AS  . MISSING 

INPUT  FILE.  .  .jen94.dat 

REWIND  INPUT  UNIT  PRIOR  TO  READING.  .  DATA.  .  .  YES 

NUMBER  OF  INTEGER  WORDS  OF  MEMORY  FOR  STORAGE  .  19998 


VARIABLES  TO  BE  USED 

1  month  2  sum 

DATA  FORMAT :  FREE 

THE  LONGEST  RECORD  MAY  HAVE  UP  TO  80  CHARACTERS. 

USE  ONLY  COMPLETE  CASES 

COMPUTE  KENDALL  RANK  CORRELATION  COEFFICIENT ( S ) 
COMPUTE  SPEARMAN  RANK  CORRELATION  COEFFICIENT ( S) 

NUMBER  OF  CASES  READ.  .  .  .• . .  12 


VARIABLE 

NO.  NAME 

MEAN 

STANDARD 

DEVIATION 

MINIMUM 

MEDIAN 

MAXIMUM 

COUNT 

1  month 

6.5000 

3.6056 

1.0000 

6.5000 

12.0000 

12 

2  sum 

1PAGE  2  3S 

528.3333 

12 -MAY- 9 6 

58.1524 

425.0000 

17:32:23 

528.5000 

613.0000 

12 

KENDALL  RANK  CORRELATION  COEFFICIENTS 


month 


1 


sum 


2 


~  ~  O  "  i  -  t.  ■>  c  -» 

J  •  ^  -  V 

month  1  1.0000 

sum  2  0.0303  1.0000 

SPEARMAN  RANK  CORRELATION  COEFFICIENTS 


month  sum 

1  2 

month  1  1.0000 

sum  2  0.0979  1.0000 

NUMBER  OF  INTEGER  WORDS  USED  IN  PRECEDING  PROBLEM  568 

CPU  TIME  USED  0.170  SECONDS 
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PROGRAM  INSTRUCTIONS 


/input 

/variables 

/test 

/end 


variables  are  2. 
format  is  free, 
file  is  'jen95.dat'. 
names  are  month,  sum. 
kendall . 
spearman . 


PROBLEM  TITLE  IS 

12 -MAY-9 6  17:34:19 


NUMBER  OF  VARIABLES  TO  READ .  2 

NUMBER  OF  VARIABLES  ADDED  BY  TRANSFORMATIONS.  .  0 

TOTAL  NUMBER  OF  VARIABLES .  2 

CASE  LABELING  VARIABLES  . 

NUMBER  OF  CASES  TO  READ . TO  END 

MISSING  VALUES  CHECKED  BEFORE  OR  AFTER  TRANS.  .  NEITHER 

BLANKS  IN  THE  DATA  ARE  TREATED  AS  . MISSING 

INPUT  FILE.  .  .jen95.dat 

REWIND  INPUT  UNIT  PRIOR  TO  READING.  .  DATA.  .  .  YES 

NUMBER  OF  INTEGER  WORDS  OF  MEMORY  FOR  STORAGE  .  19998 


VARIABLES  TO  BE  USED 

1  month  2  sum 


DATA  FORMAT :  FREE 


THE  LONGEST  RECORD  MAY  HAVE  UP  TO  80  CHARACTERS. 
USE  ONLY  COMPLETE  CASES 

COMPUTE  KENDALL  RANK  CORRELATION  COEFFICIENT (S) 
COMPUTE  SPEARMAN  RANK  CORRELATION  COEFFICIENT ( S) 


NUMBER  OF  CASES  READ 


12 


VARIABLE 
NO.  NAME 

1  month 

2  sum 

1PAGE  2  3S 


MEAN 

18.5000 
509.2501 
12 -MAY-: 


STANDARD 

DEVIATION 

3.6056 

148.4650 


MINIMUM 

13.0000 
249.0000 
:  34 : 19 


MEDIAN 


18.5000 

520.5000 


MAXIMUM  COUNT 

24.0000  12 

866.0000  12 


KENDALL  RANK  CORRELATION  COEFFICIENTS 


month 


1 


sum 


2 


month 

1 

1.0000 

sum 

2 

0.0606 

1.0000 

SPEARMAN 

RANK 

CORRELATION 

COEFFICIENTS 

month 

sum 

1 

2 

month 

1 

1.0000 

sum 

2 

0.0490 

1.0000 

NUMBER  OF  INTEGER  WORDS  USED  IN  PRECEDING  PROBLEM  568 

CPU  TIME  USED  0.140  SECONDS 

1PAGE  3  3S 


BMDP3S  -  NONPARAMETRIC  STATISTICS 

Copyright  1977,  1979,  1981,  1982,  1983,  1985,  1987,  1988,  1990,  1993 

by  BMDP  Statistical  Software,  Inc. 


BMDP  Statistical  Software,  Inc. 
12121  Wilshire  Blvd,  Suite  300 
Los  Angeles,  CA  90025  USA 
Phone  (310)  207-8800 
Fax  (310)  207-8844 


BMDP  Statistical  Software 

Cork  Technology  Park,  Model  Farm  Rd 

Cork,  Ireland 

Phone  +353  21  542722 
Fax  +353  21  542822 


Release:  7.1  (AXP/OpenVMS)  DATE:  12-MAY-96  AT  17:34:19 


PROGRAM  INSTRUCTIONS 


END  OF  INSTRUCTIONS 


PROGRAM  TERMINATED 


I£TB^  >  read  'jenifer.dat' 
Entering  data  from  file: 

24  rows  read. 

MTB  >  name  cl  'month'  c2 
MTB  >  name  c8  'July'  c9 
MTB  >  regress  cl3  12  cl- 
3UBC>  tres  c23; 

SUBC>  dw. 


cl-cl3 

jenifer.dat 

'Jan'  c3  'Feb'  c4  'Mar'  c5  'Apr'  c6  'May'  c7  'June' 
'Aug'  clO  'Sep'  ell  'Oct'  cl2  'Nov'  cl3  'sum' 
cl2  c21  c22 ; 


The  regression  equation  is 


Predictor 

Constant 

month 

Jan 

Feb 

Mar 

Apr 

May 

June 

July 

Aug 

Sep 

Oct 


59  month 

+  32  Jan  + 

165  Feb  + 

307  Mar 

.72  June  + 

96  July  +  122  Aug  + 

- —  -n**  **  -**,<  C-^o  i  *. 

183  Sep 

Coef 

St  dev 

t-ratio 

P 

411.13  ' 

98.87 

4.16 

0.002 

-1.590  l 

3.596 

-0.44 

0.667 

32.0 

112.9 

0.28 

0.782 

165.1 

111.6 

1.48 

0.167 

306.7 

110.5 

2.77 

0.018 

64.8 

109.5 

0.59 

0.566 

56.4 

108.7 

0.52 

0.614 

172.5 

107.9 

1.60 

0.138 

95.5 

107.2 

0.89 

0.392 

122.1 

106.7 

1.15 

0.276 

183.2 

106.2 

1.72 

0.113 

164.8 

105.9 

1.56 

0.148 

167.4 

105.8 

1.58 

0.142 

+  65  Apr  +  56  May 


'  *  fC  c 


% 

\  % 

*  ^ 

\  \ 

1  '  : 
V 


/. 


R-sq  =  56.4% 


R-sq(adj)  =  8.8% 


v 


*,» 


Analysis  of  Variance 


SOURCE 

DF 

SS 

MS 

Regression 

12 

158960 

13247 

Error 

11 

122884 

11171 

Total 

23 

281844 

SOURCE 

DF 

SEQ  SS 

month 

1 

5179 

Jan 

1 

22441 

Feb 

1 

999 

Mar 

1 

64844 

Apr 

1 

5407 

May 

1 

9779 

June 

1 

4000 

July 

1 

1865 

Aug 

1 

109 

Sep 

1 

7664 

Oct 

1 

8680 

Nov 

1 

27994 

Unusual  Observations 

Obs .  month 

sum 

Fit  Stc 

3  3. 

0 

541.0 

713.0 

15  15. 

0 

866.0 

694.0 

R  denotes  an 

obs . 

with  a  large 

st.  resi< 

F 

1.19 


P 

0.392 


77.8 

77.8 


Residual 

-172.0 

172.0 


St .Resid 
-2.40R 
2.40R 


Durbin-Watson  statistic  =2.48 


MTB  >  name  ,c22  'pred  sum'  c23  'stud  res'  c24  'n  score' 
MTB  >  plot  c23  c22 


stud  res- 


♦  w  + 


2.5  + 


0 . 0  + 


*  *  *  * 


*  2  * 


■2.5+ 


-pred  sum 


MTB  >  nsco  c23  c24 
MTB  >  plot  c24  c23 

n  score  - 


1.5  + 


0.0  + 


■1.5  + 


■stud  res 


MTB  >  let  k90  =  23 
MTB  >  execute  'skku' 

Executing  from  file:  skku. MTB 
MTB  >  print  k95  k96 

skewness  0.000000035 
kurtosis  2.99777  - 

MTB'  >  end 

MTB  >  execute  '.qqlog' 

Executing  from  file:  qqlog. MTB 
MTB  >  plot  c92  c90 


2.0  + 


0.0  + 


1  ■  ■'  <  -/"  / 


y?.(  .-u 


2 


*  V  * 


-2.0  + 


★ 


★ 


+  i - 1 - - 1 - f 

-2.4  -1.2  0.0  1.2  2.4 

MTB  >  end 

MTB  >  execute  'qqlap' 

Executing  from  file:  qqlap. MTB 
MTB  >  plot  c92  c90 

C92 


1.6  + 

—  *  ★ 

-  * 

-  *  * 

.  *2 

0.0+  *  2  * 

2* 

-  ** 

-  * 

-  *  * 

-1.6  + 

_  * 


* 


- + - + - + - - - + - + 

-2.4  -1.2  0.0  1.2  2.4 


-C90 


V 


N  V 
rv 


Si, 

V 

\j 

Vs  . 

<  \ 

■  V 

--  O 

v  4 


$ 


-fc 


-C90 


MTB  >  end 

MTB  >  rreg  cl3  12  cl-cl2 
The  regression  equation  is 

sum  =  411  -  1.59  month  +  32.0  Jan  +  165  Feb  +  307  Mar  +  64.8  Apr  +  56.4  May 
+  172  June  +  95.5  July  +  122  Aug  +  183  Sep  +  165  Oct  +  167  Nov 


Coefficient 

St.  dev. 

of  coef. 

Predictor 

Rank 

Least-sq 

Rank 

Least-sq 

Constant 

411.13 

411.13 

93.61 

98.87 

month 

-1.590 

-1.590 

3.405 

3.596 

Jan 

32.0 

32.0 

106.9 

112.9 

Feb 

165.1 

165.1 

105.7 

111.6 

Mar 

306.7 

306.7 

104.7 

110.5 

Apr 

64.8 

64.8 

.  103.7 

109.5 

May 

56.4 

56.4 

102.9 

108.7 

June 

172.5 

172.5 

102.1 

107.9 

July 

95.5 

95.5 

101.5 

107.2 

Aug 

122.1 

122.1 

101.0 

106.7 

Sep 

183.2 

183.2 

100.6 

106.2 

Oct 

164.8 

164.8 

100.3 

105.9 

Nov 

167.4 

167.4 

100.1 

105.8 

Hodges -Lehmann  estimate  of  tau  =  100.1  Least-squares  S  =  105.7 
MTB  >  save  ' j enif er .mtw' 

Saving  worksheet  in  file:  j enif er. mtw 
MTB  >  Is  jen* 

*  ERROR  *  Name  not  found  in  dictionary. 

MTB  >  stop 

***  Minitab  Release  9.1  ***  Minitab  Inc.  *** 

Worksheet  size:  5310509  cells 
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$VSTAT-6b 


REG  500  runs 


Decrease  Number  of  Beds  to  6 

Station  1  Registration  -  Vitals  -  Hemoglobin  ( 500  Replications/4  hours) 
(Infinite  capcity) 


Rep  No. 

L 

W 

Lq 

Wq 

Po 

Max(sys) 

Max(line) 

1 

3.4359 

13.6834 

0 

0 

0.7137 

7 

1 

2 

3.5225 

14.4145 

0 

0 

0.7065 

7 

1 

3 

3.5116 

14.3691 

0 

0 

0.7074 

7 

1 

4 

3.9133 

14.7694 

0 

0 

0.6739 

7 

1 

5 

3.6965 

13.4775 

0 

0 

0.692 

7 

1 

6 

3.4355 

14.4729 

0 

0 

0.7137 

7 

1 

7 

3.4262 

13.0524 

0 

0 

0.7145 

7 

1 

8 

3.4563 

14.0208 

0 

0 

0.712 

6 

1 

9 

3.4616 

14.0557 

0 

0 

0.7115 

8 

1 

10 

3.5865 

13.9964 

0 

0 

0.7011 

7 

1 

11 

3.1864 

13.1344 

0 

0 

0.7345 

7 

1 

488 

3.3624 

13.2171 

0 

0 

0.7198 

6 

1 

489 

3.9036 

14.9705 

0 

0 

0.6747 

8 

1 

490 

3.5461 

13.6878 

0 

0 

0.7045 

7 

1 

491 

3.5362 

14.389 

0 

0 

0.7053 

8 

1 

492 

4.0483 

13.9337 

0 

0 

0.6626 

8 

1 

493 

3.6161 

13.6527 

0 

0 

0.6987 

7 

1 

494 

3.4545 

13.1244 

0 

0 

0.7121 

7 

1 

495 

3.4497 

13.6076 

0 

0 

0.7125 

8 

1 

496 

3.6333 

13.8754 

0 

0 

0.6972 

8 

1 

497 

3.5495 

14.6695 

0 

0 

0.7042 

6 

1 

498 

3.6303 

13.9372 

0 

0 

0.6975 

8 

1 

499 

3.8103 

15.0071 

0 

0 

0.6825 

7 

1 

500 

4.021 

14.4745 

0 

0 

0.6649 

7 

1 

Average 

3.570 

14.100 

0.000 

0.000 

0.702 

7.3 

1.0 

Std  Dev 

0.237 

0.574 

0.000 

0.000 

0.020 

0.7 

0.0 

Max 

4.208 

15.747 

0.000 

0.000 

0.756 

11.0 

1.0 

Min 

2.924 

12.178 

0.000 

0.000 

0.649 

6.0 

1.0 

T(.90) 

1.730 

T(.95) 

2.090 

+/-  (.90) 

0.018 

0.044 

0.000 

0.000 

0.002 

0.1 

0.0 

+/-  (.95) 

0.022 

0.054 

0.000 

0.000 

0.002 

0.1 

0.0 

L 

W 

Lq 

Wq 

Po 

Max(sys) 

Max(line) 

LB  (.95) 

3.548 

14.046 

0.000 

0.000 

0.701 

7.2 

1.0 

LB  (.90) 

3.552 

14.055 

0.000 

0.000 

0.701 

7.2 

1.0 

Avg 

3.570 

14.100 

0.000 

0.000 

0.702 

7.3 

1.0 

UB  (.90) 

3.589 

14.144 

0.000 

0.000 

0.704 

7.3 

1.0 

UB  (.95) 

3.592 

14.153 

0.000 

0.000 

0.704 

7.4 

1.0 

B 

0.179 

0.705 

0.000 

0.000 

0.035 

0.4 

0.1 

n= 

5.273 

1.986 

0.000 

0.000 

0.946 

12.1 

0.0 
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"  \ivSTAT-6b 

i 


Int  500  runs 


Decrease  Number  of  Beds  to  6 

Station  2  -  Interview  (500  Replications/4  hours) 

(2  servers) 


Rep  No. 

L 

W 

Lq 

Wq 

Po 

Max(sys) 

Max(line) 

1 

1.2262 

4.8834 

0.1564 

0.6227 

0.4651 

5 

3 

2 

1.0331 

4.4197 

0.085 

0.3638 

0.526 

4 

2 

3 

1.2423 

5.1592 

0.1475 

0.6128 

0.4526 

5 

3 

4 

1.4219 

5.4409 

0.1839 

0.7039 

0.381 

5 

3 

5 

1.3909 

5.1396 

0.1571 

0.5806 

0.3831 

4 

2 

6 

1.2465 

5.4154 

0.1729 

0.751 

0.4632 

6 

4 

7 

1.3838 

5.4182 

0.2055 

0.8047 

0.4109 

5 

3 

8 

1.2257 

5.2052 

0.1548 

0.6575 

0.4646 

4 

2 

9 

1.3749 

5.6648 

0.217 

0.8942 

0.4211 

5 

3 

10 

1.4391 

5.6163 

0.2677 

1.0449 

0.4143 

5 

3 

11 

1.407 

6.0632 

0.2562 

1.1041 

0.4246 

5 

3 

488 

1.2628 

5.036 

0.198 

0.7898 

0.4676 

5 

3 

489 

1.3447 

5.2945 

0.2225 

0.876 

0.4389 

5 

3 

490 

1.4169 

5.7105 

0.2491 

1.0039 

0.4161 

6 

4 

491 

1.2718 

5.2499 

0.2007 

0.8284 

0.4645 

5 

3 

492 

1.4645 

5.1731 

0.2147 

0.7586 

0.3751 

5 

3 

493 

1.3603 

5.1359 

0.1706 

0.644 

0.4051 

6 

4 

494 

1.6249 

6.1733 

0.3163 

1.2018 

0.3457 

5 

3 

495 

1.4672 

5.9462 

0.3609 

1.4628 

0.4469 

6 

4 

496 

1.3499 

5.2277 

0.2372 

0.9184 

0.4436 

5 

3 

497 

1.1311 

5.0579 

0.1005 

0.4492 

0.4847 

5 

3 

498 

1.7709 

6.7988 

0.5193 

1.9937 

0.3742 

8 

6 

499 

1.1957 

4.7091 

0.1083 

0.4266 

0.4563 

3 

1 

500 

1.5123 

5.589 

0.2608 

0.9639 

0.3743 

5 

3 

Average 

1.322 

5.327 

0.211 

0.843 

0.444 

4.9 

2.9 

Std  Dev 

0.167 

0.516 

0.098 

0.369 

0.044 

0.8 

0.8 

Max 

2.138 

8.403 

0.851 

3.346 

0.564 

8.0 

2.9 

Min 

0.962 

4.310 

0.033 

0.149 

0.312 

3.0 

1.0 

T(.90) 

1.730 

T(.95) 

2.090 

+/-  (.90) 

0.013 

0.040 

0.008 

0.029 

0.003 

0.1 

0.1 

+/-  (.95) 

0.016 

0.048 

0.009 

0.034 

0.004 

0.1 

0.1 

L 

W 

Lq 

Wq 

Po 

Max(sys) 

Max(line) 

LB  (.95) 

1.307 

5.279 

0.202 

0.808 

0.440 

4.8 

2.8 

LB  (.90) 

1.310 

5.287 

0.203 

0.814 

0.441 

4.9 

2.9 

Avg 

1.322 

5.327 

0.211 

0.843 

0.444 

4.9 

2.9 

UB  (.90) 

1.335 

5.367 

0.218 

0.871 

0.448 

5.0 

3.0 

UB  (.95) 

1.338 

5.376 

0.220 

0.877 

0.448 

5.0 

3.0 

B 

0.066 

0.266 

0.011 

0.042 

0.022 

0.2 

0.1 

n= 

18.990 

11.253 

256.697 

229.179 

11.739 

28.8 

81.9 
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WsTAT-6b 


BAG  500  runs 


Decrease  Number  of  Beds  to  6 


Station  3  -  Bag  Table  (500  Replications/4  hours) 
(1  server) 


Rep  No. 

L 

W 

Lq 

Wq 

Po 

Max(sys)  Max(line) 

1 

0.7849 

3.8785 

0.2752 

1.36 

0.4903 

4 

3 

2 

0.636 

3.2652 

0.1225 

0.6288 

0.4865 

3 

2 

3 

0.7932 

3.872 

0.2325 

1.1352 

0.4393 

4 

3 

4 

0.9099 

4.0433 

0.2929 

1.3014 

0.383 

3 

2 

5 

0.8227 

3.7493 

0.2376 

1.0827 

0.4149 

3 

2 

6 

0.5097 

3.081 

0.0745 

0.4505 

0.5648 

2 

1 

7 

0.8152 

3.7067 

0.2086 

0.9486 

0.3934 

3 

2 

8 

0.8024 

3.8262 

0.2248 

1.0718 

0.4224 

4 

3 

9 

0.6675 

3.4003 

0.1464 

0.7457 

0.4789 

3 

2 

10 

0.9655 

4.6425 

0.375 

1.8032 

0.4095 

4 

3 

11 

0.6169 

3.3108 

0.1177 

0.6315 

0.5007 

3 

2 

488 

0.7928 

3.6975 

0.2222 

1.0364 

0.4294 

3 

2 

489 

0.7626 

3.7535 

0.2174 

1.0701 

0.4548 

3 

2 

490 

0.7997 

3.6529 

0.194 

0.8862 

0.3943 

4 

3 

491 

0.7654 

3.8929 

0.2036 

1.0357 

0.4383 

4 

3 

492 

0.9098 

3.7008 

0.2617 

1.0644 

0.3519 

3 

2 

493 

0.7664 

3.7722 

0.2004 

0.9864 

0.434 

3 

2 

494 

0.6845 

3.3572 

0.1396 

0.6845 

0.455 

3 

2 

495 

0.7945 

3.7909 

0.2333 

1.1134 

0.4389 

3 

2 

496 

0.7768 

3.7471 

0.2022 

0.9752 

0.4254 

4 

3 

497 

0.6064 

3.3081 

0.1188 

0.6481 

0.5124 

3 

2 

498 

0.7817 

3.5821 

0.1904 

0.8724 

0.4087 

3 

2 

499 

0.6556 

3.6361 

0.1552 

0.8606 

0.4996 

3 

2 

500 

0.8879 

4.1017 

0.2962 

1.3683 

0.4083 

4 

3 

Average 

0.768 

3.755 

0.217 

1.051 

0.449 

3.4 

2.4 

Std  Dev 

0.123 

0.418 

0.088 

0.385 

0.047 

0.7 

0.7 

Max 

1.304 

5.902 

0.686 

3.236 

0.597 

6.0 

2.4 

Min 

0.481 

2.983 

0.072 

0.425 

0.314 

2.0 

1.0 

T(.90) 

1.730 

T(.95) 

2.090 

+/-  (.90) 

0.010 

0.032 

0.007 

0.030 

0.004 

0.1 

0.1 

+/-  (.95) 

0.012 

0.039 

0.008 

0.036 

0.004 

0.1 

0.1 

L 

W 

Lq 

Wq 

Po 

Max(sys) 

Max(line) 

LB  (.95) 

0.757 

3.716 

0.209 

1.015 

0.444 

3.3 

2.3 

LB  (.90) 

0.759 

3.722 

0.210 

1.021 

0.445 

3.3 

2.3 

Avg 

0.768 

3.755 

0.217 

1.051 

0.449 

3.4 

2.4 

UB  (.90) 

0.778 

3.787 

0.224 

1.081 

0.452 

3.4 

2.4 

UB  (.95) 

0.780 

3.794 

0.225 

1.087 

0.453 

3.4 

2.4 

B 

0.038 

0.188 

0.011 

0.053 

0.022 

0.2 

0.1 

n= 

30.706 

14.823 

197.956 

160.866 

13.288 

44.9 

90.8 

Pagel 


\A/STAT-6b 


BLD  500  runs 


l 


Decrease  Number  of  Beds  to  6 

Station  4  -  Blood  Letting  ( 500  Replications/4  hours) 


(6  servers) 

Rep  No. 

L 

W 

Lq 

Wq 

Po 

Max(sys)  Max(line) 

1 

3 

13.863 

0.0095 

0.0467 

0.534 

7 

1 

2 

3 

13.0498 

0.0021 

0.0108 

0.5767 

7 

1 

3 

3 

13.1797 

0 

0 

0.55 

6 

1 

4 

3 

14.9933 

0.0433 

0.1923 

0.4449 

8 

2 

5 

3 

13.6247 

0 

0 

0.5017 

6 

1 

6 

3 

15.3382 

0 

0 

0.5771 

6 

1 

7 

3 

12.0292 

0 

0 

0.5591 

5 

1 

8 

3 

15.0517 

0.0125 

0.0595 

0.476 

7 

1 

9 

3 

13.6665 

0 

0 

0.5529 

6 

1 

10 

3 

13.3536 

0.0038 

0.0182 

0.5378 

7 

1 

11 

2 

12.116 

0 

0 

0.6237 

6 

1 

488 

3 

12.7265 

0.0031 

0.0146 

0.5457 

7 

1 

489 

3 

14.3097 

0.0067 

0.0327 

0.5165 

7 

1 

490 

3 

13.5938 

0 

0 

0.504 

5 

1 

491 

3 

16.2275 

0.0386 

0.1962 

0.4747 

8 

2 

492 

4 

15.0383 

0.0554 

0.2253 

0.3931 

8 

2 

493 

3 

14.3897 

0.0185 

0.0911 

0.5158 

7 

1 

494 

3 

13.1176 

0.0105 

0.0515 

0.556 

7 

1 

495 

3 

16.2521 

0.1256 

0.5993 

0.4533 

9 

3 

496 

3 

12.2989 

0 

0 

0.5751 

6 

1 

497 

3 

13.6707 

0 

0 

0.5823 

5 

1 

498 

3 

14.2573 

0.0039 

0.0177 

0.4821 

7 

1 

499 

2 

11.4854 

0 

0 

0.6549 

5 

1 

500 

3 

12.736 

0 

0 

0.5405 

6 

1 

Average 

2.888 

13.861 

0.011 

0.054 

0.531 

6.7 

1.2 

Std  Dev 

0.357 

0.972 

0.020 

0.097 

0.049 

0.9 

0.5 

Max 

4.000 

16.753 

0.166 

0.958 

0.668 

9.0 

1.2 

Min 

2.000 

11.400 

0.000 

0.000 

0.367 

5.0 

1.0 

T(.90) 

1.730 

T(.95) 

2.090 

+/-  (.90) 

0.028 

0.075 

0.002 

0.008 

0.004 

0.1 

0.0 

+/-  (.95) 

0.033 

0.091 

0.002 

0.009 

0.005 

0.1 

0.0 

L 

W 

Lq 

Wq 

Po 

Max(sys) 

Max(line) 

LB  (.95) 

2.855 

13.770 

0.009 

0.045 

0.526 

6.6 

1.1 

LB  (.90) 

2.860 

13.785 

0.010 

0.046 

0.527 

6.7 

1.2 

Avg 

2.888 

13.861 

0.011 

0.054 

0.531 

6.7 

1.2 

UB  (.90) 

2.916 

13.936 

0.013 

0.061 

0.534 

6.8 

1.2 

UB  (.95) 

2.921 

13.951 

0.013 

0.063 

0.535 

6.8 

1.2 

B 

0.144 

0.693 

0.001 

0.003 

0.027 

0.3 

0.1 

n= 

18.331 

5.893 

3756.632 

3933.691 

10.399 

20.1 

173.9 
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WSTAT-6b 


TOT  500  runs 


Decrease  Number  of  Beds  to  6 

Total  System  (from  GPSSH) 

(6  servers) 


Rep  No. 

L 

W 

Lmax 

#XACTS 

1 

8 

32.8659 

14 

67 

2 

8 

31.6468 

12 

69 

3 

8 

33.7457 

13 

68 

4 

10 

36.3039 

16 

73 

5 

9 

32.4478 

15 

75 

6 

8 

32.5618 

13 

66 

7 

8 

31.5083 

13 

74 

8 

9 

35.0532 

12 

67 

9 

8 

33.2424 

13 

69 

10 

9 

34.2182 

15 

69 

11 

7 

30.7836 

12 

69 

488 

8 

32.0242 

13 

70 

489 

9 

34.2027 

15 

77 

490 

9 

33.7316 

13 

71 

491 

9 

35.6601 

14 

70 

492 

10 

34.8302 

16 

78 

493 

9 

32.721 

13 

73 

494 

8 

32.06 

14 

71 

495 

9 

35.9641 

14 

75 

496 

8 

31.7336 

12 

72 

497 

8 

32.207 

12 

66 

498 

9 

35.6826 

15 

74 

499 

8 

30.4547 

13 

69 

500 

9 

33.0386 

16 

77 

Average 

8.476 

33.514 

13.946 

70.6 

Std  Dev 

0.712 

1.602 

1.329 

3.6 

T(.90) 

1.730 

T(.95) 

2.1 

+/-  (.90) 

0.055 

0.620 

0.514 

1.4 

+/-  (.95) 

0.067 

0.748 

0.621 

1.7 

L 

W 

Lmax 

#XACTS 

LB  (.95) 

8.409 

32.765 

13.325 

68.9 

LB  (.90) 

8.421 

32.894 

13.432 

69.2 

Avg 

8.476 

33.514 

13.946 

70.6 

UB  (.90) 

8.531 

34.133 

14.460 

71.9 

UB  (.95) 

8.543 

34.262 

14.567 

72.2 

B 

0.424 

1.676 

0.697 

3.5279 

n= 

8.439 

2.734 

10.874 

3.0 
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Station  1  Registration  -  Vitals  -  Hemoglobin  ( 500  Replications/4  hours)  3.57  14.10  0.00  0.00  (Infinite  capcity) 

Station  2  -  Interview  (500  Replications/4  hours)  1.32  5.33  0.21  0.84  (2  servers) 

Station  3  -  Bag  Table  (500  Replications/4  hours)  0.77  3.75  0.22  1.05  (1  server) 

Station  4  -  Blood  Letting  ( 500  Replications/4  hours)  2.89  13.86  0.01  0.05  (6  servers) 
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Wstat-ex 


REG  500  runs 


Exponential  Interarrivals  and  9  beds 

Station  1  Registration  -  Vitals  -  Hemoglobin  ( 500  Replications/4  hours) 
(Infinite  capcity) 


No. 

L 

W 

Lq 

Wq 

Po 

Max(sys) 

Max(line) 

1 

3.4359 

13.6834 

0 

0 

0.7137 

7 

1 

2 

3.5225 

14.4145 

0 

0 

0.7065 

7 

1 

3 

3.5116 

14.3691 

0 

0 

0.7074 

7 

1 

4 

3.9133 

14.7694 

0 

0 

0.6739 

7 

1 

5 

3.6965 

13.4775 

0 

0 

0.692 

7 

1 

6 

3.4355 

14.4729 

0 

0 

0.7137 

7 

1 

7 

3.4262 

13.0524 

0 

0 

0.7145 

7 

1 

8 

3.4563 

14.0208 

0 

0 

0.712 

6 

1 

9 

3.4616 

14.0557 

0 

0 

0.7115 

8 

1 

10 

3.5865 

13.9964 

0 

0 

0.7011 

7 

1 

11 

3.1864 

13.1344 

0 

0 

0.7345 

7 

1 

488 

3.3624 

13.2171 

0 

0 

0.7198 

6 

1 

489 

3.9036 

14.9705 

0 

0 

0.6747 

8 

1 

490 

3.5461 

13.6878 

0 

0 

0.7045 

7 

1 

491 

3.5362 

14.389 

0 

0 

0.7053 

8 

1 

492 

4.0483 

13.9337 

0 

0 

0.6626 

8 

1 

493 

3.6161 

13.6527 

0 

0 

0.6987 

7 

1 

494 

3.4545 

13.1244 

0 

0 

0.7121 

7 

1 

495 

3.4497 

13.6076 

0 

0 

0.7125 

8 

1 

496 

3.6333 

13.8754 

0 

0 

0.6972 

8 

1 

497 

3.5495 

14.6695 

0 

0 

0.7042 

6 

1 

498 

3.6303 

13.9372 

0 

0 

0.6975 

8 

1 

499 

3.8103 

15.0071 

0 

0 

0.6825 

7 

1 

500 

4.021 

14.4745 

0 

0 

0.6649 

7 

1 

Average 
Std  Dev 

3.570 

0.237 

14.100 

0.574 

0.000 

0.000 

0.000 

0.000 

0.702 

0.020 

7.3 

0.7 

1.0 

0.0 

Max 

4.208 

15.747 

0.000 

0.000 

0.756 

11.0 

1.0 

Min 

2.924 

12.178 

0.000 

0.000 

0.649 

6.0 

1.0 

T(.90) 

1.730 

T(.95) 

2.090 

+/-  (.90) 

0.018 

0.044 

0.000 

0.000 

0.002 

0.1 

0.0 

+/-  (.95) 

0.022 

0.054 

0.000 

0.000 

0.002 

0.1 

0.0 

L 

W 

Lq 

Wq 

Po  Max(sys)  Max(Iine)| 

LB  (.95) 

3.548 

14.046 

0.000 

0.000 

0.701 

7.2 

1.0 

LB  (.90) 

3.552 

14.055 

0.000 

0.000 

0.701 

7.2 

1.0 

Avg 

3.570 

14.100 

0.000 

0.000 

0.702 

7.3 

1.0 

UB  (.90) 

3.589 

14.144 

0.000 

0.000 

0.704 

7.3 

1.0 

UB  (.95) 

3.592 

14.153 

0.000 

0.000 

0.704 

7.4 

1.0 

B 

0.179 

0.705 

0.000 

0.000 

0.035 

0.4 

0.1 

n= 

5.273 

1.986 

0.000 

0.000 

0.946 

12.1 

0.0 
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Wstat-ex 


Int  500  runs 


Exponential  Interarrivals  and  9  beds 
Station  2  -  Interview  (500  Replications/4  hours) 
(2  servers) 


Rep  No. 

L 

W 

Lq 

Wq 

Po 

Max(sys) 

Max(line) 

1 

1.2262 

4.8834 

0.1564 

0.6227 

0.4651 

5 

3 

2 

1.0331 

4.4197 

0.085 

0.3638 

0.526 

4 

2 

3 

1.2423 

5.1592 

0.1475 

0.6128 

0.4526 

5 

3 

4 

1.4219 

5.4409 

0.1839 

0.7039 

0.381 

5 

3 

5 

1.3909 

5.1396 

0.1571 

0.5806 

0.3831 

4 

2 

6 

1.2465 

5.4154 

0.1729 

0.751 

0.4632 

6 

4 

7 

1.3838 

5.4182 

0.2055 

0.8047 

0.4109 

5 

3 

8 

1.2257 

5.2052 

0.1548 

0.6575 

0.4646 

4 

2 

9 

1.3749 

5.6648 

0.217 

0.8942 

0.4211 

5 

3 

10 

1.4391 

5.6163 

0.2677 

1.0449 

0.4143 

5 

3 

11 

1.407 

6.0632 

0.2562 

1.1041 

0.4246 

5 

3 

488 

1.2628 

5.036 

0.198 

0.7898 

0.4676 

5 

3 

489 

1.3447 

5.2945 

0.2225 

0.876 

0.4389 

5 

3 

490 

1.4169 

5.7105 

0.2491 

1.0039 

0.4161 

6 

4 

491 

1.2718 

5.2499 

0.2007 

0.8284 

0.4645 

5 

3 

492 

1.4645 

5.1731 

0.2147 

0.7586 

0.3751 

5 

3 

493 

1.3603 

5.1359 

0.1706 

0.644 

0.4051 

6 

4 

494 

1.6249 

6.1733 

0.3163 

1.2018 

0.3457 

5 

3 

495 

1.4672 

5.9462 

0.3609 

1.4628 

0.4469 

6 

4 

496 

1.3499 

5.2277 

0.2372 

0.9184 

0.4436 

5 

3 

497 

1.1311 

5.0579 

0.1005 

0.4492 

0.4847 

5 

3 

498 

1.7709 

6.7988 

0.5193 

1.9937 

0.3742 

8 

6 

499 

1.1957 

4.7091 

0.1083 

0.4266 

0.4563 

3 

1 

500 

1.5123 

5.589 

0.2608 

0.9639 

0.3743 

5 

3 

Average 

1.322 

5.327 

0.211 

0.843 

0.444 

4.9 

2.9 

Std  Dev 

0.167 

0.516 

0.098 

0.369 

0.044 

0.8 

0.8 

Max 

2.138 

8.403 

0.851 

3.346 

0.564 

8.0 

2.9 

Min 

0.962 

4.310 

0.033 

0.149 

0.312 

3.0 

1.0 

T(.90) 

1.730 

T(.95) 

2.090 

+/-  (.90) 

0.013 

0.040 

0.008 

0.029 

0.003 

0.1 

0.1 

+/-  (.95) 

0.016 

0.048 

0.009 

0.034 

0.004 

0.1 

0.1 

L 

W 

Lq 

Wq 

Po 

Max(sys) 

Max(line) 

LB  (.95) 

1.307 

5.279 

0.202 

0.808 

0.440 

4.8 

2.8 

LB  (.90) 

1.310 

5.287 

0.203 

0.814 

0.441 

4.9 

2.9 

Avg 

1.322 

5.327 

0.211 

0.843 

0.444 

4.9 

2.9 

UB  (.90) 

1.335 

5.367 

0.218 

0.871 

0.448 

5.0 

3.0 

UB  (.95) 

1.338 

5.376 

0.220 

0.877 

0.448 

5.0 

3.0 

B 

0.066 

0.266 

0.011 

0.042 

0.022 

0.2 

0.1 

n= 

18.990 

11.253 

256.697 

229.179 

11.739 

28.8 

81.9 

Page  1 


'Wstat-ex 


BAG  500  runs 


Exponential  Interarrivals  and  9  beds 
Station  3  -  Bag  Table  (500  Replications/4  hours) 
(1  server) 


Rep  No. 

L 

W 

Lq 

Wq 

Po 

Max(sys)  Max(line) 

1 

0.7849 

3.8785 

0.2752 

1.36 

0.4903 

4 

3 

2 

0.636 

3.2652 

0.1225 

0.6288 

0.4865 

3 

2 

3 

0.7932 

3.872 

0.2325 

1.1352 

0.4393 

4 

3 

4 

0.9099 

4.0433 

0.2929 

1.3014 

0.383 

3 

2 

5 

0.8227 

3.7493 

0.2376 

1.0827 

0.4149 

3 

2 

6 

0.5097 

3.081 

0.0745 

0.4505 

0.5648 

2 

1 

7 

0.8152 

3.7067 

0.2086 

0.9486 

0.3934 

3 

2 

8 

0.8024 

3.8262 

0.2248 

1.0718 

0.4224 

4 

3 

9 

0.6675 

3.4003 

0.1464 

0.7457 

0.4789 

3 

2 

10 

0.9655 

4.6425 

0.375 

1.8032 

0.4095 

4 

3 

11 

0.6169 

3.3108 

0.1177 

0.6315 

0.5007 

3 

2 

488 

0.7928 

3.6975 

0.2222 

1.0364 

0.4294 

3 

2 

489 

0.7626 

3.7535 

0.2174 

1.0701 

0.4548 

3 

2 

490 

0.7997 

3.6529 

0.194 

0.8862 

0.3943 

4 

3 

491 

0.7654 

3.8929 

0.2036 

1.0357 

0.4383 

4 

3 

492 

0.9098 

3.7008 

0.2617 

1.0644 

0.3519 

3 

2 

493 

0.7664 

3.7722 

0.2004 

0.9864 

0.434 

3 

2 

494 

0.6845 

3.3572 

0.1396 

0.6845 

0.455 

3 

2 

495 

0.7945 

3.7909 

0.2333 

1.1134 

0.4389 

3 

2 

496 

0.7768 

3.7471 

0.2022 

0.9752 

0.4254 

4 

3 

497 

0.6064 

3.3081 

0.1188 

0.6481 

0.5124 

3 

2 

498 

0.7817 

3.5821 

0.1904 

0.8724 

0.4087 

3 

2 

499 

0.6556 

3.6361 

0.1552 

0.8606 

0.4996 

3 

2 

500 

0.8879 

4.1017 

0.2962 

1.3683 

0.4083 

4 

3 

Average 

0.768 

3.755 

0.217 

1.051 

0.449 

3.4 

2.4 

Std  Dev 

0.123 

0.418 

0.088 

0.385 

0.047 

0.7 

0.7 

Max 

1.304 

5.902 

0.686 

3.236 

0.597 

6.0 

2.4 

Min 

0.481 

2.983 

0.072 

0.425 

0.314 

2.0 

1.0 

T(.90) 

1.730 

T(.95) 

2.090 

+/-  (.90) 

0.010 

0.032 

0.007 

0.030 

0.004 

0.1 

0.1 

+/-  (.95) 

0.012 

0.039 

0.008 

0.036 

0.004 

0.1 

0.1 

L 

W 

Lq 

Wq 

Po 

Max(sys) 

Max(line) 

LB  (.95) 

0.757 

3.716 

0.209 

1.015 

0.444 

3.3 

2.3 

LB  (.90) 

0.759 

3.722 

0.210 

1.021 

0.445 

3.3 

2.3 

Avg 

0.768 

3.755 

0.217 

1.051 

0.449 

3.4 

2.4 

UB  (.90) 

0.778 

3.787 

0.224 

1.081 

0.452 

3.4 

2.4 

UB  (.95) 

0.780 

3.794 

0.225 

1.087 

0.453 

3.4 

2.4 

B 

0.038 

0.188 

0.011 

0.053 

0.022 

0.2 

0.1 

n= 

30.706 

14.823 

197.956 

160.866 

13.288 

44.9 

90.8 
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Wstat-ex 


BLD  500  runs 


Exponential  Interarrivals  and  9  beds 
Station  4  -  Blood  Letting  ( 500  Replications/4  hours) 


(6  servers) 

Rep  No. 

L 

W 

Lq 

Wq 

Po 

Max(sys)  Max(line) 

1 

3 

13.863 

0.0095 

0.0467 

0.534 

7 

1 

2 

3 

13.0498 

0.0021 

0.0108 

0.5767 

7 

1 

3 

3 

13.1797 

0 

0 

0.55 

6 

1 

4 

3 

14.9933 

0.0433 

0.1923 

0.4449 

8 

2 

5 

3 

13.6247 

0 

0 

0.5017 

6 

1 

6 

3 

15.3382 

0 

0 

0.5771 

6 

1 

7 

3 

12.0292 

0 

0 

0.5591 

5 

1 

8 

3 

15.0517 

0.0125 

0.0595 

0.476 

7 

1 

9 

3 

13.6665 

0 

0 

0.5529 

6 

1 

10 

3 

13.3536 

0.0038 

0.0182 

0.5378 

7 

1 

11 

2 

12.116 

0 

0 

0.6237 

6 

1 

488 

3 

12.7265 

0.0031 

0.0146 

0.5457 

7 

1 

489 

3 

14.3097 

0.0067 

0.0327 

0.5165 

7 

1 

490 

3 

13.5938 

0 

0 

0.504 

5 

1 

491 

3 

16.2275 

0.0386 

0.1962 

0.4747 

8 

2 

492 

4 

15.0383 

0.0554 

0.2253 

0.3931 

8 

2 

493 

3 

14.3897 

0.0185 

0.0911 

0.5158 

7 

1 

494 

3 

13.1176 

0.0105 

0.0515 

0.556 

7 

1 

495 

3 

16.2521 

0.1256 

0.5993 

0.4533 

9 

3 

496 

3 

12.2989 

0 

0 

0.5751 

6 

1 

497 

3 

13.6707 

0 

0 

0.5823 

5 

1 

498 

3 

14.2573 

0.0039 

0.0177 

0.4821 

7 

1 

499 

2 

11.4854 

0 

0 

0.6549 

5 

1 

500 

3 

12.736 

0 

0 

0.5405 

6 

1 

Average 

2.888 

13.861 

0.011 

0.054 

0.531 

6.7 

1.2 

Std  Dev 

0.357 

0.972 

0.020 

0.097 

0.049 

0.9 

0.5 

Max 

4.000 

16.753 

0.166 

0.958 

0.668 

9.0 

1.2 

Min 

2.000 

11.400 

0.000 

0.000 

0.367 

5.0 

1.0 

T(.90) 

1.730 

T(.95) 

2.090 

+/-  (.90) 

0.028 

0.075 

0.002 

0.008 

0.004 

0.1 

0.0 

+/-  (.95) 

0.033 

0.091 

0.002 

0.009 

0.005 

0.1 

0.0 

L 

W 

Lq 

Wq 

Po 

Max(sys) 

Max(line) 

LB  (.95) 

2.855 

13.770 

0.009 

0.045 

0.526 

6.6 

1.1 

LB  (.90) 

2.860 

13.785 

0.010 

0.046 

0.527 

6.7 

1.2 

Avg 

2.888 

13.861 

0.011 

0.054 

0.531 

6.7 

1.2 

UB  (.90) 

2.916 

13.936 

0.013 

0.061 

0.534 

6.8 

1.2 

UB  (.95) 

2.921 

13.951 

0.013 

0.063 

0.535 

6.8 

1.2 

B 

0.144 

0.693 

0.001 

0.003 

0.027 

0.3 

0.1 

n= 

18.331 

5.893 

3756.632 

3933.691 

10.399 

20.1 

173.9 
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Wstat-ex 


TOT  500  runs 


Exponential  Interarrivals  and  9  beds 

Total  System  (from  GPSSH) 

(6  servers) 


Rep  No. 

L 

W 

Lmax 

#XACTS 

1 

8 

32.8659 

14 

67 

2 

8 

31.6468 

12 

69 

3 

8 

33.7457 

13 

68 

4 

10 

36.3039 

16 

73 

5 

9 

32.4478 

15 

75 

6 

8 

32.5618 

13 

66 

7 

8 

31.5083 

13 

74 

8 

9 

35.0532 

12 

67 

9 

8 

33.2424 

13 

69 

10 

9 

34.2182 

15 

69 

11 

7 

30.7836 

12 

69 

488 

8 

32.0242 

13 

70 

489 

9 

34.2027 

15 

77 

490 

9 

33.7316 

13 

71 

491 

9 

35.6601 

14 

70 

492 

10 

34.8302 

16 

78 

493 

9 

32.721 

13 

73 

494 

8 

32.06 

14 

71 

495 

9 

35.9641 

14 

75 

496 

8 

31.7336 

12 

72 

497 

8 

32.207 

12 

66 

498 

9 

35.6826 

15 

74 

499 

8 

30.4547 

13 

69 

500 

9 

33.0386 

16 

77 

Average 

8.476 

33.514 

13.946 

70.6 

Std  Dev 

0.712 

1.602 

1.329 

3.6 

T(.90) 

1.730 

T(.95) 

2.1 

+/-  (.90) 

0.055 

0.620 

0.514 

1.4 

+/-  (.95) 

0.067 

0.748 

0.621 

1.7 

L 

W 

Lmax 

#XACTS 

LB  (.95) 

8.409 

32.765 

13.325 

68.9 

LB  (.90) 

8.421 

32.894 

13.432 

69.2 

Avg 

8.476 

33.514 

13.946 

70.6 

UB  (.90) 

8.531 

34.133 

14.460 

71.9 

UB  (.95) 

8.543 

34.262 

14.567 

72.2 

B 

0.424 

1.676 

0.697 

3.5279 

n= 

8.439 

2.734 

10.874 

3.0 
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Station  1  Registration  -  Vitals  -  Hemoglobin  ( 500  Replications/4  hours)  3.57  14.10  0.00  0.00  (Infinite  capcity) 

Station  2  -  Interview  (500  Replications/4  hours)  1.32  5.33  0.21  0.84  (2  servers) 

Station  3  -  Bag  Table  (500  Replications/4  hours)  0.77  3.75  0.22  1.05  (1  server) 

Station  4  -  Blood  Letting  ( 500  Replications/4  hours)  2.89  13.86  0.01  0.05  (6  servers) 
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Sums  500  runs 


(Infinite  capcity) 
(2  servers) 

(1  server) 

(9  servers) 
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Wstatbas 


REG  500  runs 


Base  Case 

Station  1  Registration  -  Vitals  -  Hemoglobin  ( 500  Replications/4  hours) 
(Infinite  capcity) 


Rep  No. 

L 

W 

Lq 

Wq 

Po 

Max(sys)  Max(line) 

1 

3.3421 

13.7609 

0 

0 

0.7215 

7 

1 

2 

3.3344 

14.7102 

0 

0 

0.7221 

7 

1 

3 

3.5628 

14.3258 

0 

0 

0.7031 

8 

1 

4 

3.4712 

13.96 

0 

0 

0.7107 

7 

1 

5 

3.2547 

14.0702 

0 

0 

0.7288 

6 

1 

6 

3.2259 

14.3327 

0 

0 

0.7312 

7 

1 

7 

3.0904 

12.9921 

0 

0 

0.7425 

8 

1 

8 

3.2716 

14.1147 

0 

0 

0.7274 

7 

1 

9 

3.5991 

14.3755 

0 

0 

0.7001 

8 

1 

10 

3.3052 

13.4354 

0 

0 

0.7246 

8 

1 

11 

3.0486 

13.1302 

0 

0 

0.746 

7 

1 

488 

3.6787 

14.9604 

0 

0 

0.6934 

7 

1 

489 

3.3049 

13.5241 

0 

0 

0.7246 

7 

1 

490 

3.6422 

14.7991 

0 

0 

0.6965 

7 

1 

491 

3.2679 

14.034 

0 

0 

0.7277 

7 

1 

492 

2.7077 

12.6561 

0 

0 

0.7744 

6 

1 

493 

3.1866 

13.5184 

0 

0 

0.7345 

7 

1 

494 

3.0829 

13.446 

0 

0 

0.7431 

7 

1 

495 

3.2605 

13.9349 

0 

0 

0.7283 

7 

1 

496 

3.7407 

14.8558 

0 

0 

0.6883 

8 

1 

497 

2.9255 

13.9929 

0 

0 

0.7562 

7 

1 

498 

3.6583 

14.6385 

0 

0 

0.6951 

8 

1 

499 

3.6192 

14.8387 

0 

0 

0.6984 

7 

1 

500 

3.4097 

13.8527 

0 

0 

0.7159 

7 

1 

Average 

3.350 

14.097 

0.000 

0.000 

0.721 

7.3 

1.0 

Std  Dev 

0.221 

0.563 

0.000 

0.000 

0.018 

0.7 

0.0 

Max 

3.994 

16.068 

0.000 

0.000 

0.774 

11.0 

1.0 

Min 

2.708 

12.467 

0.000 

0.000 

0.667 

6.0 

1.0 

T(-90) 

1.730 

T(.95) 

2.090 

+/-  (.90) 

0.017 

0.044 

0.000 

0.000 

0.001 

0.1 

0.0 

+/-  (.95) 

0.021 

0.053 

0.000 

0.000 

0.002 

0.1 

0.0 

L 

W 

Lq 

Wq 

Po 

Max(sys) 

Max(line) 

LB  (.95) 

3.329 

14.045 

0.000 

0.000 

0.719 

7.3 

1.0 

LB  (.90) 

3.333 

14.054 

0.000 

0.000 

0.719 

7.3 

1.0 

Avg 

3.350 

14.097 

0.000 

0.000 

0.721 

7.3 

1.0 

UB  (.90) 

3.367 

14.141 

0.000 

0.000 

0.722 

7.4 

1.0 

UB  (.95) 

3.370 

14.150 

0.000 

0.000 

0.723 

7.4 

1.0 

B 

0.167 

0.705 

0.000 

0.000 

0.036 

0.4 

0.1 

n= 

5.224 

1.910 

0.000 

0.000 

0.783 

11.3 

0.0 
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Wstatbas 


Int  500  runs 


Base  Case 

Station  2  -  Interview  (500  Replications/4  hours) 
(2  servers) 


Rep  No. 

L 

w 

Lq 

Wq 

Po 

Max(sys)  Max(line) 

1 

3.3301 

13.7116 

2.3013 

9.4753 

0.2847 

13 

11 

2 

2.8494 

13.1334 

1.9441 

8.9607 

0.3783 

11 

10 

3 

2.863 

11.9852 

1.7677 

7.3998 

0.2942 

13 

12 

4 

2.8846 

11.7598 

1.7351 

7.0736 

0.296 

14 

12 

5 

2.1586 

9.4733 

1.1114 

4.8776 

0.3435 

8 

8 

6 

2.5775 

11.6339 

1.5817 

7.139 

0.3366 

10 

8 

7 

2.6825 

11.2772 

1.5657 

6.5821 

0.2905 

8 

8 

8 

3.11 

13.612 

2.0692 

9.0566 

0.2755 

12 

12 

9 

3.5202 

14.2528 

2.3148 

9.3723 

0.2324 

12 

10 

10 

3.4434 

13.9974 

2.3332 

9.4845 

0.2816 

12 

10 

11 

2.8477 

12.4429 

1.7707 

7.7371 

0.2999 

10 

10 

488 

2.4309 

10.1681 

1.2772 

5.3423 

0.2745 

10 

10 

489 

3.3154 

13.9339 

2.2683 

9.5332 

0.3192 

14 

12 

490 

3.1272 

12.7068 

2.0412 

8.2939 

0.2969 

11 

11 

491 

3.5826 

15.8251 

2.5801 

11.397 

0.3196 

15 

14 

492 

2.7992 

13.4924 

1.76 

8.4836 

0.3112 

13 

11 

493 

2.0815 

8.9641 

1.0356 

4.4599 

0.329 

8 

7 

494 

2.4054 

10.6501 

1.3916 

6.1615 

0.3544 

10 

9 

495 

2.2964 

9.8144 

1.233 

5.2696 

0.3171 

9 

8 

496 

2.9728 

12.1338 

1.834 

7.4856 

0.275 

13 

12 

497 

2.3126 

11.0613 

1.4488 

6.9297 

0.4054 

13 

11 

498 

2.8818 

11.6913 

1.747 

7.0877 

0.2702 

10 

9 

499 

3.6293 

15.7821 

2.5431 

11.0586 

0.25 

11 

10 

500 

2.8222 

11.6299 

1.6757 

6.9054 

0.2459 

10 

9 

Average 

2.812 

12.077 

1.768 

7.591 

0.322 

11.0 

10.0 

Std  Dev 

0.474 

1.905 

0.440 

1.820 

0.042 

1.6 

1.5 

Max 

Min 

5.148 

1.852 

20.789 

7.772 

3.943 

0.833 

16.089 

3.567 

0.446 

0.205 

16.0 

7.0 

10.0 

5.0 

T(.90) 

1.730 

T(.95) 

2.090 

+/-  (.90) 

0.037 

0.147 

0.034 

0.141 

0.003 

0.1 

0.1 

+/-  (.95) 

0.044 

0.178 

0.041 

0.170 

0.004 

0.2 

0.1 

L 

W 

Lq 

Wq 

Po  Max(sys)  Max(line)| 

LB  (.95) 

2.768 

11.899 

1.727 

7.421 

0.318 

10.9 

9.8 

LB  (.90) 

2.775 

11.930 

1.734 

7.451 

0.319 

10.9 

9.9 

Avg 

2.812 

12.077 

1.768 

7.591 

0.322 

11.0 

10.0 

UB  (.90) 

2.849 

12.224 

1.802 

7.732 

0.325 

11.2 

10.1 

UB  (.95) 

2.857 

12.255 

1.809 

7.762 

0.326 

11.2 

10.1 

B 

0.141 

0.604 

0.088 

0.380 

0.016 

0.6 

0.5 

n= 

34.069 

29.777 

74.307 

68.782 

19.919 

25.3 

27.9 
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Wstatbas 


BAG  500  runs 


Base  Case 

Station  3  -  Bag  Table  (500  Replications/4  hours) 
(1  server) 


Rep  No. 

L 

w 

Lq 

Wq 

Po 

Max(sys)  Max(iine) 

1 

0.7679 

4.2953 

0.3125 

1.7481 

0.5446 

4 

3 

2 

0.681 

3.824 

0.2018 

1.1331 

0.5208 

4 

3 

3 

1.235 

6.6214 

0.7049 

3.7793 

0.4699 

7 

6 

4 

0.911 

4.6745 

0.3465 

1.7779 

0.4355 

4 

3 

5 

0.8248 

4.049 

0.2724 

1.3371 

0.4476 

4 

3 

6 

0.867 

4.5655 

0.3465 

1.8247 

0.4795 

5 

4 

7 

1.0201 

5.1759 

0.439 

2.2276 

0.4189 

4 

3 

8 

0.9182 

4.8649 

0.395 

2.093 

0.4768 

5 

4 

9 

0.9086 

4.4026 

0.3741 

1.8125 

0.4654 

4 

3 

10 

1.1295 

5.3931 

0.5378 

2.5677 

0.4083 

4 

3 

11 

0.7034 

4.0781 

0.2291 

1.3281 

0.5257 

3 

2 

488 

0.832 

4.1291 

0.2787 

1.3833 

0.4467 

4 

3 

489 

0.8046 

3.9718 

0.2869 

1.4163 

0.4823 

4 

3 

490 

0.9488 

4.6767 

0.4033 

1.9882 

0.4546 

4 

3 

491 

0.7874 

4.1977 

0.2829 

1.508 

0.4955 

4 

3 

492 

0.7183 

4.1809 

0.2483 

1.445 

0.53 

4 

3 

493 

0.8781 

4.303 

0.3292 

1.6134 

0.4512 

5 

4 

494 

0.7008 

3.7925 

0.2217 

1.1999 

0.5209 

4 

3 

495 

0.8495 

4.3315 

0.2921 

1.4894 

0.4426 

3 

2 

496 

0.808 

3.9574 

0.2601 

1.2741 

0.4522 

3 

2 

497 

0.7665 

4.6392 

0.3147 

1.9045 

0.5481 

5 

4 

498 

0.9937 

4.8376 

0.4185 

2.0376 

0.4249 

4 

3 

499 

0.8742 

4.7338 

0.3585 

1.9414 

0.4843 

4 

3 

500 

0.6648 

3.4865 

0.1898 

0.9956 

0.5251 

3 

2 

Average 

0.864 

4.458 

0.342 

1.754 

0.478 

4.1 

3.1 

Std  Dev 

0.183 

0.767 

0.155 

0.734 

0.042 

0.9 

0.9 

Max 

1.522 

7.730 

0.965 

5.075 

0.604 

8.0 

3.1 

Min 

0.530 

3.170 

0.112 

0.580 

0.369 

3.0 

2.0 

T(-90) 

1.730 

T(.95) 

2.090 

+/-  (.90) 

0.014 

0.059 

0.012 

0.057 

0.003 

0.1 

0.1 

+/-  (.95) 

0.017 

0.072 

0.015 

0.069 

0.004 

0.1 

0.1 

L 

W 

Lq 

Wq 

Po 

Max(sys) 

Max(line) 

LB  (.95) 

0.847 

4.386 

0.328 

1.686 

0.474 

4.0 

3.0 

LB  (.90) 

0.850 

4.398 

0.330 

1.698 

0.475 

4.1 

3.1 

Avg 

0.864 

4.458 

0.342 

1.754 

0.478 

4.1 

3.1 

UB  (.90) 

0.878 

4.517 

0.354 

1.811 

0.482 

4.2 

3.2 

UB  (.95) 

0.881 

4.529 

0.357 

1.823 

0.482 

4.2 

3.2 

B 

0.043 

0.223 

0.017 

0.088 

0.024 

0.2 

0.2 

n= 

53.666 

35.474 

245.995 

209.662 

9.294 

59.4 

103.5 
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Wstatbas 


BLO  500  runs 


Base  Case 


Station  4  - 

Blood  Letting  ( 500  Replications/4  hours) 

(9  servers) 

Rep  No. 

L 

w 

Lq 

Wq 

Po 

Max(sys)  Max(line) 

1 

2 

13.9007 

0 

0 

-0.2761 

8 

1 

2 

2 

12.9271 

0 

0 

-0.2558 

6 

1 

3 

2 

13.2711 

0 

0 

-0.275 

6 

1 

4 

3 

14.3214 

0 

0 

-0.3101 

7 

1 

5 

3 

14.2888 

0 

0 

-0.3234 

7 

1 

6 

3 

14.7458 

0 

0 

-0.3111 

8 

1 

7 

2 

12.1404 

0 

0 

-0.2659 

5 

1 

8 

3 

15.0505 

0 

0 

-0.3156 

8 

1 

9 

3 

13.3733 

0 

0 

-0.3067 

8 

1 

10 

3 

13.2941 

0 

0 

-0.3094 

7 

1 

11 

2 

12.5645 

0 

0 

-0.2408 

6 

1 

488 

3 

13.4618 

0 

0 

-0.3014 

7 

1 

489 

3 

14.4144 

0 

0 

-0.3244 

7 

1 

490 

3 

14.8326 

0 

0 

-0.3343 

9 

1 

491 

2 

12.2994 

0 

0 

-0.2563 

8 

1 

492 

2 

14.3364 

0 

0 

-0.2737 

7 

1 

493 

3 

14.2133 

0 

0 

-0.3223 

7 

1 

494 

2 

10.9595 

0 

0 

-0.225 

5 

1 

495 

3 

13.1014 

0 

0 

-0.2855 

6 

1 

496 

3 

13.5562 

0 

0 

-0.3075 

7 

1 

497 

2 

13.4661 

0 

0 

-0.2472 

6 

1 

498 

3 

13.3013 

0 

0 

-0.3036 

6 

1 

499 

2 

11.9878 

0 

0 

-0.246 

6 

1 

500 

3 

14.3283 

0 

0 

-0.3036 

8 

1 

Average 

2.742 

13.808 

0.000 

0.001 

•0.296 

7.2 

1.0 

Std  Dev 

0.438 

0.929 

0.001 

0.006 

0.028 

0.9 

0.1 

Max 

3.000 

16.709 

0.017 

0.090 

-0.199 

11.0 

1.0 

Min 

2.000 

10.892 

0.000 

0.000 

-0.379 

5.0 

1.0 

T(.90) 

1.730 

T(.95) 

2.090 

+/-  (.90) 

0.034 

0.072 

0.000 

0.000 

0.002 

0.1 

0.0 

+/-  (.95) 

0.041 

0.087 

0.000 

0.001 

0.003 

0.1 

0.0 

L 

W 

Lq 

Wq 

Po 

Max(sys) 

Max(line) 

LB  (.95) 

2.701 

13.721 

0.000 

0.000 

-0.299 

7.1 

1.0 

LB  (.90) 

2.708 

13.736 

0.000 

0.000 

-0.298 

7.1 

1.0 

Avg 

2.742 

13.808 

0.000 

0.001 

-0.296 

7.2 

1.0 

UB  (.90) 

2.776 

13.880 

0.000 

0.001 

-0.294 

7.3 

1.0 

UB  (.95) 

2.783 

13.895 

0.000 

0.001 

-0.293 

7.3 

1.0 

B 

0.137 

0.690 

0.000 

0.000 

-0.015 

0.4 

0.1 

n= 

30.543 

5.417  mtimrn  ######## 

10.780 

20.2 

4.7 
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4 

Wstatbas 


TOT  500  runs 


Base  Case 

Total  System  (from  GPSSH) 
(9  servers) 


Rep  No. 

L 

W 

Lmax 

#XACTS 

1 

8 

32.8283 

14 

67 

2 

8 

31.6382 

12 

69 

3 

8 

33.7457 

13 

68 

4 

10 

36.1406 

16 

73 

5 

9 

32.4478 

15 

75 

6 

8 

32.5618 

13 

66 

7 

8 

31.5083 

13 

74 

8 

9 

35.0026 

12 

67 

9 

8 

33.2424 

13 

69 

10 

9 

34.2035 

15 

69 

11 

7 

30.7836 

12 

69 

488 

8 

32.0119 

13 

70 

489 

9 

34.1772 

15 

77 

490 

9 

33.7316 

13 

71 

491 

9 

35.5031 

14 

70 

492 

10 

34.6396 

16 

78 

493 

9 

32.6511 

13 

73 

494 

8 

32.0201 

14 

71 

495 

9 

35.4686 

14 

75 

496 

8 

31.7336 

12 

72 

497 

8 

32.207 

12 

66 

498 

9 

35.6678 

15 

74 

499 

8 

30.4547 

13 

69 

500 

9 

33.0386 

16 

77 

Average 

8.464 

33.470 

13.942 

70.6 

Std  Dev 

0.708 

1.580 

1.329 

3.6 

T(.90) 

1.730 

T(.95) 

2.1 

+/-  (.90) 

0.055 

0.611 

0.514 

1.4 

+/-  (.95) 

0.066 

0.739 

0.621 

1.7 

L 

W 

Lmax 

#XACTS 

LB  (.95) 

8.398 

32.731 

13.321 

68.9 

LB  (.90) 

8.409 

32.858 

13.428 

69.2 

Avg 

8.464 

33.470 

13.942 

70.6 

UB  (.90) 

8.519 

34.081 

14.456 

71.9 

UB  (.95) 

8.530 

34.208 

14.563 

72.2 

B 

0.423 

1.673 

0.697 

3.5279 

n= 

8.384 

2.669 

10.878 

3.0 
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Sums  500  runs 


Station  1  Registration  -  V 

L 

W 

Lq 

Wq 

Po  Max(sys)  Max(line) 

3.35 

14.10 

0.00 

0.00 

Station  2  -  Interview  (50 

2.81 

12.08 

1.77 

7.59 

Station  3  -  Bag  Table  (5 

0.86 

4.46 

0.34 

1.75 

Station  4  -  Blood  Letting 

2.74 

13.81 

0.00 

0.00 

9.77 

44.44 

2.11 

9.35 

Page  1 1 


Sums  500  runs 


(Infinite  capcity) 
(2  servers) 

(1  server) 

(9  servers) 


Page  2 


V 


base48hrs 


REG  500  runs 


Two  Interviewers  at  all  Times 

Station  1  Registration  -  Vitals  -  Hemoglobin  ( 500  Replications/48  hours) 
(Infinite  capcity) 


Rep  No. 

L 

W 

Lq 

Wq 

Po 

Max(sys)  Max(line) 

1 

4.0343 

14.6855 

0.0002 

0.0006 

0.6638 

13 

1 

2 

4.2876 

15.0588 

0 

0 

0.6427 

12 

1 

3 

3.9541 

14.0628 

0 

0 

0.6705 

11 

1 

4 

3.8724 

13.467 

0 

0 

0.6773 

11 

1 

5 

3.6955 

13.7784 

0 

0 

0.692 

10 

1 

6 

3.8842 

13.7764 

0 

0 

0.6763 

12 

1 

7 

3.7278 

13.7715 

0 

0 

0.6893 

11 

1 

8 

4.1851 

14.5569 

0 

0 

0.6512 

12 

1 

9 

4.0297 

14.5527 

0 

0 

0.6642 

12 

1 

10 

3.6475 

13.4331 

0 

0 

0.696 

12 

1 

11 

4.143 

14.2216 

0 

0 

0.6547 

12 

1 

488 

4.0303 

14.4689 

0 

0 

0.6641 

12 

1 

489 

4.7263 

15.0573 

0.0032 

0.0103 

0.6064 

15 

3 

490 

4.1159 

14.3162 

0.0065 

0.0225 

0.6575 

14 

2 

491 

3.7305 

13.6895 

0.0008 

0.003 

0.6892 

13 

1 

492 

3.9733 

14.6158 

0 

0 

0.6689 

11 

1 

493 

4.0941 

14.5137 

0 

0 

0.6588 

12 

1 

494 

4.0704 

14.0225 

0.0043 

0.0149 

0.6612 

15 

3 

495 

3.8456 

13.8788 

0 

0 

0.6795 

12 

1 

496 

4.2672 

14.4932 

0 

0 

0.6444 

12 

1 

497 

4.1393 

13.9222 

0.0001 

0.0004 

0.6551 

13 

1 

498 

4.1096 

14.2597 

0 

0 

0.6575 

11 

1 

499 

3.7109 

13.3991 

0 

0 

0.6908 

11 

1 

500 

4.4978 

14.9695 

0 

0 

0.6252 

12 

1 

Average 

4.027 

14.127 

0.001 

0.002 

0.664 

12.1 

1.2 

Std  Dev 

0.198 

0.510 

0.002 

0.007 

0.017 

1.3 

0.6 

Max 

4.726 

16.318 

0.025 

0.088 

0.712 

17.0 

1.2 

Min 

3.462 

12.666 

0.000 

0.000 

0.606 

10.0 

1.0 

T(.90) 

1.730 

T(.95) 

2.090 

+/-  (.90) 

0.015 

0.039 

0.000 

0.001 

0.001 

0.1 

0.0 

+/-  (.95) 

0.019 

0.048 

0.000 

0.001 

0.002 

0.1 

0.1 

L 

W 

Lq 

Wq 

Po 

Max(sys) 

Max(line) 

LB  (.95) 

4.008 

14.080 

0.000 

0.002 

0.663 

12.0 

1.2 

LB  (.90) 

4.012 

14.088 

0.000 

0.002 

0.663 

12.0 

1.2 

Avg 

4.027 

14.127 

0.001 

0.002 

0.664 

12.1 

1.2 

UB  (.90) 

4.042 

14.167 

0.001 

0.003 

0.666 

12.2 

1.3 

UB  (.95) 

4.045 

14.175 

0.001 

0.003 

0.666 

12.2 

1.3 

B 

0.201 

0.706 

0.000 

0.000 

0.033 

0.6 

0.1 

n= 

2.908 

1.562 

mmtmm 

0.740 

14.1 

289.9 
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base48hrs 


Int  500  runs 


Two  Interviewers  at  all  Times 


Station  2  -  Interview  (500  Replications/48  hours) 
(2  servers) 


Rep  No. 

L 

W 

Lq 

Wq 

Po 

Max(sys) 

Max(line) 

1 

1.5042 

6.4821 

0.4016 

1.7307 

0.4487 

9 

7 

2 

1.4161 

5.8936 

0.4067 

1.6927 

0.4953 

12 

10 

3 

1.3999 

5.8015 

0.3752 

1.5548 

0.4876 

10 

8 

4 

1.751 

7.1472 

0.6228 

2.5423 

0.4359 

10 

8 

5 

1.5383 

6.6743 

0.494 

2.1434 

0.4779 

14 

12 

6 

1.4546 

6.0018 

0.3509 

1.4476 

0.4481 

8 

6 

7 

1.4866 

6.3517 

0.4153 

1.7746 

0.4644 

9 

7 

8 

1.597 

6.4418 

0.4599 

1.855 

0.4314 

11 

9 

9 

1.4808 

6.3003 

0.4486 

1.9085 

0.4839 

12 

10 

10 

1.299 

5.6512 

0.3255 

1.4161 

0.5133 

10 

8 

11 

1.6934 

6.9769 

0.5848 

2.4093 

0.4457 

13 

11 

488 

1.4441 

6.1658 

0.3402 

1.4524 

0.4481 

9 

7 

489 

1.6915 

6.4098 

0.4971 

1.8838 

0.4028 

9 

7 

490 

1.448 

6.0792 

0.3953 

1.6597 

0.4736 

12 

10 

491 

1.3925 

6.1241 

0.4038 

1.7757 

0.5056 

10 

8 

492 

1.2844 

5.6869 

0.279 

1.2354 

0.4973 

7 

5 

493 

1.6235 

6.6099 

0.5126 

2.0868 

0.4445 

9 

7 

494 

1.5592 

6.4893 

0.5277 

2.196 

0.4842 

13 

11 

495 

1.4436 

6.2519 

0.4072 

1.7637 

0.4818 

13 

11 

496 

1.5766 

6.2782 

0.4833 

1.9244 

0.4533 

11 

9 

497 

1.589 

6.2365 

0.4848 

1.9026 

0.4479 

11 

9 

498 

1.8011 

7.3683 

0.6796 

2.7801 

0.4392 

12 

10 

499 

1.42 

6.0497 

0.3529 

1.5036 

0.4665 

8 

6 

500 

1.4503 

5.7498 

0.3871 

1.5345 

0.4684 

8 

6 

Average 

1.528 

6.318 

0.448 

1.847 

0.460 

10.0 

8.0 

Std  Dev 

0.180 

0.638 

0.133 

0.516 

0.029 

2.2 

2.2 

Max 

2.377 

9.310 

1.096 

4.175 

0.544 

20.0 

8.0 

Min 

1.122 

4.770 

0.210 

0.893 

0.350 

6.0 

4.0 

T(.90) 

1.730 

T(.95) 

2.090 

+/-  (.90) 

0.014 

0.049 

0.010 

0.040 

0.002 

0.2 

0.2 

+/-  (.95) 

0.017 

0.060 

0.012 

0.048 

0.003 

0.2 

0.2 

L 

W 

Lq 

Wq 

Po 

Max(sys) 

Max(line) 

LB  (.95) 

1.511 

6.258 

0.435 

1.798 

0.457 

9.8 

7.8 

LB  (.90) 

1.514 

6.268 

0.437 

1.807 

0.458 

9.9 

7.9 

Avg 

1.528 

6.318 

0.448 

1.847 

0.460 

10.0 

8.0 

UB  (.90) 

1.541 

6.367 

0.458 

1.887 

0.462 

10.2 

8.2 

UB  (.95) 

1.544 

6.377 

0.460 

1.895 

0.463 

10.2 

8.2 

B 

0.076 

0.316 

0.022 

0.092 

0.023 

0.5 

0.4 

n= 

16.578 

12.197 

105.670 

93.596 

4.695 

55.2 

86.0 
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base48hrs  BAG  500  runs 


Two  Interviewers  at  all  Times 

Station  3  -  Bag  Table  (500  Replications/4  hours) 

(1  server) 


Rep  No. 

L 

W 

Lq 

Wq 

Po 

Max(sys)  Max(line) 

1 

2.362 

10.4284 

1.6761 

7.3998 

0.314 

14 

13 

2 

2.6168 

11.0994 

1.9128 

8.1131 

0.296 

15 

14 

3 

2.2927 

9.6965 

1.5778 

6.6731 

0.2851 

12 

11 

4 

3.9098 

16.2812 

3.1332 

13.0474 

0.2234 

15 

14 

5 

2.0213 

8.9429 

1.33 

5.8846 

0.3088 

13 

12 

6 

2.9795 

12.5819 

2.2454 

9.4822 

0.266 

14 

13 

7 

2.1948 

9.6338 

1.5172 

6.6597 

0.3224 

12 

11 

8 

2.8905 

11.8416 

2.1309 

8.7297 

0.2404 

14 

13 

9 

1.7733 

7.7252 

1.1219 

4.8872 

0.3485 

9 

8 

10 

2.4907 

10.9179 

1.7913 

7.8524 

0.3007 

16 

15 

11 

2.1636 

9.1099 

1.4623 

6.157 

0.2987 

11 

10 

488 

1.8567 

8.1567 

1.1905 

5.2302 

0.3338 

11 

10 

489 

3.686 

14.2684 

2.8829 

11.1594 

0.1968 

18 

17 

490 

2.4482 

10.508 

1.7038 

7.3131 

0.2556 

12 

11 

491 

1.9335 

8.701 

1.2911 

5.8103 

0.3577 

14 

13 

492 

1.792 

8.0941 

1.1733 

5.2995 

0.3813 

12 

11 

493 

2.2634 

9.3458 

1.568 

6.4746 

0.3047 

17 

16 

494 

2.0579 

8.7286 

1.367 

5.7982 

0.3091 

12 

11 

495 

2.1623 

9.5659 

1.4493 

6.4116 

0.287 

11 

10 

496 

2.0511 

8.4958 

1.3401 

5.5508 

0.289 

11 

10 

497 

3.623 

14.4949 

2.8477 

11.3931 

0.2247 

16 

15 

498 

2.3807 

9.9083 

1.691 

7.0377 

0.3103 

13 

12 

499 

2.167 

9.3569 

1.4797 

6.3889 

0.3126 

16 

15 

500 

2.731 

11.0243 

1.9716 

7.9587 

0.2406 

11 

10 

Average 

2.423 

10.217 

1.719 

7.238 

0.296 

14.2 

13.2 

Std  Dev 

0.583 

2.245 

0.553 

2.176 

0.038 

3.7 

3.7 

Max 

5.127 

20.766 

4.316 

17.481 

0.401 

32.0 

13.2 

Min 

1.459 

6.403 

0.829 

3.639 

0.175 

8.0 

7.0 

T(.90) 

1.730 

T(.95) 

2.090 

+/-  (.90) 

0.045 

0.174 

0.043 

0.168 

0.003 

0.3 

0.3 

+/-  (.95) 

0.054 

0.210 

0.052 

0.203 

0.004 

0.3 

0.3 

L 

W 

Lq 

Wq 

Po 

Max(sys) 

Max(line) 

LB  (.95) 

2.369 

10.008 

1.667 

7.034 

0.292 

13.8 

12.8 

LB  (.90) 

2.378 

10.044 

1.676 

7.069 

0.293 

13.9 

12.9 

Avg 

2.423 

10.217 

1.719 

7.238 

0.296 

14.2 

13.2 

UB  (.90) 

2.468 

10.391 

1.762 

7.406 

0.299 

14.5 

13.5 

UB  (.95) 

2.478 

10.427 

1.771 

7.441 

0.299 

14.5 

13.5 

B 

0.121 

0.511 

0.086 

0.362 

0.015 

0.7 

0.7 

n= 

69.188 

57.776 

123.977 

108.234 

20.104 

80.6 

93.3 
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4  * 

,  base48hrs 


BLD  500  runs 


Two  Interviewers  at  all  Times 

Station  4  -  Blood  Letting  ( 500  Replications/48  hours) 

(9  servers) 


Rep  No. 

L 

W 

Lq 

Wq 

Po 

Max(sys)  Max(line) 

1 

4.7751 

21.2446 

0.0111 

0.0493 

0.4707 

12 

3 

2 

5.005 

21.2289 

0.1179 

0.4999 

0.457 

16 

7 

3 

5.0877 

21.5174 

0.1103 

0.4664 

0.4469 

15 

6 

4 

5.5234 

23.034 

0.126 

0.5253 

0.4003 

14 

5 

5 

4.7615 

21.0667 

0.0662 

0.2927 

0.4783 

15 

6 

6 

4.8217 

20.603 

0.082 

0.3503 

0.4734 

15 

6 

7 

4.7656 

20.9497 

0.0747 

0.3283 

0.4788 

14 

5 

8 

4.6899 

19.2957 

0.0449 

0.1846 

0.4839 

14 

5 

9 

5.0448 

21.9765 

0.0936 

0.4077 

0.4499 

14 

5 

10 

4.8437 

21.2649 

0.039 

0.1713 

0.4662 

12 

3 

11 

5.0458 

21.2453 

0.0946 

0.3982 

0.4499 

15 

6 

488 

5.0361 

22.1585 

0.067 

0.295 

0.4479 

14 

5 

489 

5.4469 

21.0848 

0.2094 

0.8104 

0.418 

16 

7 

490 

5.0064 

21.5843 

0.086 

0.3708 

0.4533 

14 

5 

491 

4.7835 

21.5938 

0.0373 

0.1682 

0.4726 

13 

4 

492 

4.5466 

20.5678 

0.038 

0.1718 

0.499 

14 

5 

493 

5.2003 

21.4728 

0.1714 

0.7075 

0.4412 

21 

12 

494 

5.3934 

22.8762 

0.0779 

0.3305 

0.4094 

13 

4 

495 

4.6068 

20.3802 

0.0524 

0.2316 

0.494 

15 

6 

496 

5.1395 

21.4416 

0.0819 

0.3416 

0.438 

14 

5 

497 

5.6658 

22.7307 

0.3657 

1.4672 

0.4111 

21 

12 

498 

5.2464 

22.1548 

0.1891 

0.7984 

0.4381 

17 

8 

499 

4.8792 

21.2259 

0.0945 

0.4111 

0.4684 

14 

5 

500 

5.6501 

22.9686 

0.2155 

0.8759 

0.3961 

16 

7 

Average 

5.048 

21.426 

0.100 

0.420 

0.450 

14.6 

5.6 

Std  Dev 

0.305 

0.937 

0.068 

0.279 

0.029 

2.1 

2.1 

Max 

5.950 

24.715 

0.490 

1.980 

0.532 

23.0 

5.6 

Min 

4.247 

18.504 

0.004 

0.017 

0.366 

10.0 

1.0 

T(.90) 

1.730 

T(.95) 

2.090 

+/-  (.90) 

0.024 

0.072 

0.005 

0.022 

0.002 

0.2 

0.2 

+/-  (.95) 

0.028 

0.088 

0.006 

0.026 

0.003 

0.2 

0.2 

L 

W 

Lq 

Wq 

Po 

Max(sys) 

Max(line) 

LB  (.95) 

5.020 

21.338 

0.093 

0.394 

0.447 

14.4 

5.4 

LB  (.90) 

5.025 

21.353 

0.094 

0.398 

0.448 

14.5 

5.5 

Avg 

5.048 

21.426 

0.100 

0.420 

0.450 

14.6 

5.6 

UB  (.90) 

5.072 

21.498 

0.105 

0.441 

0.452 

14.8 

5.8 

UB  (.95) 

5.077 

21.513 

0.106 

0.446 

0.453 

14.8 

5.8 

B 

0.252 

1.071 

0.005 

0.021 

0.023 

0.7 

0.3 

n= 

4.357 

2.288 

557.372 

529.732 

5.062 

25.6 

172.2 
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base48hrs 


TOT  500  runs 


Two  Interviewers  at  all  Times 

Total  System  (from  GPSSH) 

(9  servers) 


Rep  No. 

L 

W 

Lmax 

#XACTS 

1 

8 

32.8283 

14 

67 

2 

8 

31.6382 

12 

69 

3 

8 

33.7457 

13 

68 

4 

10 

36.1406 

16 

73 

5 

9 

32.4478 

15 

75 

6 

8 

32.5618 

13 

66 

7 

8 

31.5083 

13 

74 

8 

9 

35.0026 

12 

67 

9 

8 

33.2424 

13 

69 

10 

9 

34.2035 

15 

69 

11 

7 

30.7836 

12 

69 

488 

8 

32.0119 

13 

70 

489 

9 

34.1772 

15 

77 

490 

9 

33.7316 

13 

71 

491 

9 

35.5031 

14 

70 

492 

10 

34.6396 

16 

78 

493 

9 

32.6511 

13 

73 

494 

8 

32.0201 

14 

71 

495 

9 

35.4686 

14 

75 

496 

8 

31.7336 

12 

72 

497 

8 

32.207 

12 

66 

498 

9 

35.6678 

15 

74 

499 

8 

30.4547 

13 

69 

500 

9 

33.0386 

16 

77 

Average 

8.464 

33.470 

13.942 

70.6 

Std  Dev 

0.708 

1.580 

1.329 

3.6 

T(.90) 

1.730 

T(.95) 

2.1 

+/-  (.90) 

0.055 

0.611 

0.514 

1.4 

+/-  (.95) 

0.066 

0.739 

0.621 

1.7 

L 

W 

Lmax 

#XACTS 

LB  (.95) 

8.398 

32.731 

13.321 

68.9 

LB  (.90) 

8.409 

32.858 

13.428 

69.2 

Avg 

8.464 

33.470 

13.942 

70.6 

UB  (.90) 

8.519 

34.081 

14.456 

71.9 

UB  (.95) 

8.530 

34.208 

14.563 

72.2 

B 

0.423 

1.673 

0.697 

3.5279 

n= 

8.384 

2.669 

10.878 

3.0 
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I 


L 

W 

Station  1  Registration  -  V 

4.03 

14.13 

Station  2  -  Interview  (50 

1.53 

6.32 

Station  3  -  Bag  Table  (5 

2.42 

10.22 

Station  4  -  Blood  Letting 

5.05 

21.43 

13.03 

52.09 

Wq  Po  Max(sys)  Max(line) 

0.00  0.00 

0.45  1.85 

1.72  7.24 

0.10  0.42 

2.27  9.51 


Sums  500  runs 
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